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

Environmental photosensitizers can exhibit enhanced actinic absorption in microhydrated clusters compared to solution

Brown carbon chromophores at environmental air-water interfaces often act as photosensitizers that absorb sunlight and subsequently transfer energy to nearby molecules, initiating a wide variety of chemical reactions. Despite their importance to understanding daytime chemistry at these airwater interfaces, little is known about the role of the solvation environment on the photophysical properties of these photosensitizers. In this work, we present a joint experimental-theoretical study of the vibrational and photophysical properties of microhydrated protonated and deprotonated 4-benzoylbenzoic acid (4-BBA), a key model system for environmental photosensitizers. We find that for protonated 4-BBAH + ∙(H 2 O) 0-1 , representing photosensitizers in very acidic conditions, a single bright state dominates the UV-Vis spectrum between 280 and 400 nm. Comparing the experimental UV-Vis spectra and quantum chemistry-predicted spectra of 4-BBA + ∙(H 2 O) 0-2 , we find that the degree of microhydration has little effect on the UV-Vis spectra or the orbitals of the dominant feature. For deprotonated 4-BBA ‒ , representing photosensitizers in basic conditions, quantum chemistry calculations predict that the UV-Vis spectra are ~3x weaker in intensity than the brightest 4-BBAH + ∙(H 2 O) 0-1 features and were not observed experimentally. Quantum chemistry calculations predict one dominant UV-Vis feature is present in the spectra of 4-BBA ‒ ∙(H 2 O) 0-2 , which exhibit minor shifts with degree of microhydration. While 4-BBA in bulk solution over a range of pH values has relatively weak absorption within the solar actinic region, we show that microhydrated 4-BBA has bright transitions within the actinic region. This indicates that the complex structure of environmental air-water interfaces can shift the absorption maximum of photosensitizers into the actinic region for enhanced absorption of sunlight and subsequent enhancement of photosensitizer-driven reactions.

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

Influence of fixed charge concentration and water uptake on ion sorption in AMPS/PEGDA membranes

Water content and charge density strongly influence ion sorption in charged membranes. However, isolating the individual effect of each parameter on ion sorption can be difficult due to experimental challenges associated with independently controlling these two variables during membrane synthesis. This frustrates development and validation of fundamental models and structure/property relationships. In this study, three series of uncharged and negatively charged hydrogel membranes were prepared using poly(ethylene glycol) diacrylate (PEGDA) as a cross-linker and 2-acrylamido-2-methylpropanesulfonic acid (AMPS) as an ionic monomer. Here, water content and charge density were controlled by varying AMPS content and PEGDA chain length. Ion sorption in these materials was measured using aqueous NaCl solutions at concentrations ranging from 0.01 to 1.0 M. At constant charge density, salt sorption increased as water content increased. At fixed water content, salt sorption decreased as charge density increased. Salt sorption in membranes with the highest water content and charge density was predicted using the Donnan/Manning model, which accounts for solution and membrane phase thermodynamic non-idealities. However, this approach could not predict sorption in weakly charged and uncharged membranes, presumably due to additional membrane phase thermodynamic non-idealities associated with ion solvation and uncharged polymer segments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ca-dimers, solvent layering, and dominant electrochemically active species in Ca(BH4)2 in THF

Abstract Divalent ions (Mg, Ca, and Zn) are being considered as competitive, safe, and earth-abundant alternatives to Li-ion electrochemistry, but present challenges for stable cycling due to undesirable interfacial phenomena. We explore the formation of electroactive species in the electrolyte Ca(BH 4 ) 2 ∣THF using molecular dynamics coupled with a continuum model of bulk and interfacial speciation. Free-energy analysis and unsupervised learning indicate a majority population of neutral Ca dimers and monomers with diverse molecular conformations and an order of magnitude lower concentration of the primary electroactive charged species – the monocation, $${\rm{CaBH}}_{4}^{+}$$ CaBH 4 + – produced via disproportionation of neutral complexes. Dense layering of THF molecules within ~1 nm of the electrode surface strongly modulates local electrolyte species populations. A dramatic increase in monocation population in this interfacial zone is induced at negative bias. We see no evidence for electrochemical activity of fully-solvated Ca 2+ . The consequences for performance are discussed in light of this molecular-scale insight.

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

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine‐Residue Low Complexity Peptide, as Probed by 2 H Solid‐State NMR

Abstract Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid‐state NMR to assess rotameric inter‐conversions and other dynamic modes of the methionine in the context of a nine‐residue random‐coil peptide (RC9) with the low‐complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine‐containing systems, characterized previously. Namely, the M35 of hydrated amyloid‐β 1–40 in the three‐fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)‐methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

Vugmeyster, Liliya↗
DOE OSTI2020

Gating of Hydrophobic Nanopores with Large Anions

In this report, understanding ion transport in nanoporous materials is critical to a wide variety of energy and environmental technologies, ranging from ion-selective membranes, drug delivery, and biosensing, to ion batteries and supercapacitors. While nanoscale transport is often described by continuum models that rely on a point charge description for ions and a homogeneous dielectric medium for the solvent, here, we show that transport of aqueous solutions at a hydrophobic interface can be highly dependent on the size and hydration strength of the solvated ions. Specifically, measurements of ion current through single silicon nitride nanopores that contain a hydrophobic-hydrophilic junction show that transport properties are dependent not only on applied voltage but also the type of anion. We find that in Cl - -contained solutions the nanopores only conducted ionic current above a negative voltage threshold. On the other hand, introduction of large polarizable anions, such as Br - and I - , facilitated the pore wetting, making the pore conductive at all examined voltages. Molecular dynamics simulations revealed that the large anions, Br - and I - , have a weaker solvation shell compared to that of Cl - , and consequently were prone to migrate from the aqueous solution to the hydrophobic surface, leading to the anion accumulation responsible for pore wetting. The results are essential for designing nanoporous systems that are selective to ions of the same charge, for realization of ion-induced wetting in hydrophobic pores, as well as for fundamental understanding on the role of ion hydration shell on the properties of solid/liquid interfaces.

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

Two-bond 13 C– 13 C spin-coupling constants in saccharides: dependencies on exocyclic hydroxyl group conformation

Seven doubly 13 C-labeled isotopomers of methyl β-D-glucopyranoside, methyl β-D-xylopyranoside, methyl β-D-galactopyranoside, methyl β-D-galactopyranosyl-(1→4)-β-D-glucopyranoside and methyl β-D-galactopyranosyl-(1→4)-β-D-xylopyranoside were prepared, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C NMR spectroscopy to determine experimentally the dependence of 2 J C1,C3 values in aldopyranosyl rings on the C1–C2–O2–H torsion angle, θ 2 , involving the C2 carbon of the C1–C2–C3 coupling pathway. Using X-ray crystal structures to determine θ 2 in crystalline samples and by selecting compounds that exhibit a relatively wide range of θ 2 values in the crystalline state, 2 J C1,C3 values measured in crystalline samples were plotted against θ 2 and the resulting plot compared to that obtained from density functional theory (DFT) calculations. For θ 2 values ranging from ~90° to ~240°, very good agreement was observed between the experimental and theoretical plots, providing strong validation of DFT-calculated spin-coupling dependencies on exocyclic C–O bond conformation involving the central carbon of geminal C–C–C coupling pathways. we report that these findings provide new experimental evidence supporting the use of 2 J CCC values as non-conventional spin-coupling constraints in MA'AT conformational modeling of saccharides in solution, and the use of NMR spin-couplings not involving coupled hydroxyl hydrogens as indirect probes of C–O bond conformation. Solvomorphism was observed in crystalline βGal-(1→4)-βGlcOCH 3 wherein the previously-reported methanol solvate form was found to spontaneously convert to a monohydrate upon air-drying, leading to small but discernible conformational changes in, and a new crystalline form of, this disaccharide.

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

Luminescent Bimetallic Two-Coordinate Gold(I) Complexes Utilizing Janus Carbenes

A series of bimetallic carbene-metal-amide (cMa) complexes have been prepared with bridging biscarbene ligands to serve as a model for the design of luminescent materials with large oscillator strengths and small energy differences between the singlet and triplet states (ΔE ST ). The complexes have a general structure (R 2 N)Au(:carbene—carbene:)Au(NR 2 ). The bimetallic complexes show solvation-dependent absorption and emission that is analyzed in detail. It is found that the molar absorptivity of the bimetallic complexes is correlated with the energy barrier to rotation of the metal–ligand bond. The bimetallic cMa complexes also exhibit short emission lifetimes (τ = 200–300 ns) with high photoluminescence efficiencies (Φ PL > 95%). The radiative rates of bimetallic cMa complexes are 3–4 times faster than that of the corresponding monometallic complexes. Analysis of temperature-dependent luminescence data indicates that the lifetime for the singlet state (τS 1 ) of bimetallic cMa complexes is near 12 ns with a ΔE ST of 40–50 meV. Here the presented compounds provide a general design for cMa complexes to achieve small values for ΔE ST while retaining high radiative rates. Solution-processed organic light-emitting devices (OLEDs) made using two of the complexes as luminescent dopants show high efficiency and low roll-off at high luminance.

14 SOLAR ENERGY↗
DOE OSTI2022

Coarse-grained molecular dynamics integrated with convolutional neural network for comparing shapes of temperature sensitive bottlebrushes

Quantification of shape changes in nature-inspired soft material architectures of stimuli-sensitive polymers is critical for controlling their properties but is challenging due to their softness and flexibility. Here, we have computationally designed uniquely shaped bottlebrushes of a thermosensitive polymer, poly(N-isopropylacrylamide) (PNIPAM), by controlling the length of side chains along the backbone. Coarse-grained molecular dynamics simulations of solvated bottlebrushes were performed below and above the lower critical solution temperature of PNIPAM. Conventional analyses (free volume, asphericity, etc.) show that lengths of side chains and their immediate environments dictate the compactness and bending in these architectures. We further developed 100 unique convolutional neural network models that captured molecular-level features and generated a statistically significant quantification of the similarity between different shapes. Thus, our study provides insights into the shapes of complex architectures as well as a general method to analyze them. The shapes presented here may inspire the synthesis of new bottlebrushes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Correction to Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

This errata is to correct the acknowledgment in the original published article. G.B. (partial support), V.-A.G., M.-T.N., and R.R. were supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. Computer resources were provided by Research Computing at Pacific Northwest National Laboratory and the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. PNNL is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830. G.B. acknowledges the EPSRC for partial support in PhD funding. The authors acknowledge the use of the IRIDIS High Performance Computing Facility (IRIDIS 5) and associated support services at the University of Southampton in the completion of this work. We are grateful to the UK Materials and Molecular Modeling Hub (Thomas HPC) for computational resources, which is partially funded by EPSRC (EP/P020194/1).

continuum solvation, metal/water interfaces, ab in↗
DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Selectivity in gas–liquid interactions: Molecular beam scattering of CD 4 and ND 3 from an aqueous flat liquid jet

The dynamics of polar and nonpolar molecules colliding with an aqueous surface are characterized by scattering molecular beams of deuterated methane and ammonia, CD 4 and ND 3 (E i = 28.9 and 30.3 kJ mol −1 , respectively), from a flat liquid jet of cold salty water (8 m LiBr, 230 K). Translational energy distributions of scattered species collected as a function of collision geometry probe both impulsive scattering (IS) and thermal desorption (TD) mechanisms. Here, we find that CD 4 scattering is dominated by IS and exhibits a super-specular angular distribution. The fraction of TD scattering events is notably smaller for cold salty water than for dodecane, consistent with a higher free energy of solvation for CD 4 in the water jet. In contrast, no scattering signal is seen for ND 3 from the water jet, a result attributed to the high solubility and efficient protonation of ND 3 in liquid water. The IS channel for CD 4 was analyzed using a soft-sphere model, yielding a higher internal energy (E int ) and lower effective surface mass (m eff ) than was seen for Ne/water; the higher value of E int is attributed to rotational excitation of the scattered CD 4 . These findings demonstrate that the outcomes of a gas–liquid collision—scattering trajectory, surface adherence, and energy transfer—are directed at the molecular level by both the gaseous scatterer and liquid surface.

Foreman, Madison M. [University of California, Ber↗
DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Theoretical Model for Computing Freezing Point Depression of Lithium-Ion Battery Electrolytes

Reliable prediction of freezing point depression in liquid electrolytes will accelerate the development of improved Li-ion batteries which can operate in low temperature environments. In this work we establish a computational methodology to calculate activity coefficients and liquidus lines for battery-relevant liquid electrolytes. Electronic structure methods are used in conjuction with classical molecular dynamics simulations and theoretical expressions for Born solvation energy, ion-atmosphere effects from Debye-Hückel theory and solvent entropic effects. The framework uses no a priori knowledge beyond neat solvent properties and the concentration of salt. LiPF 6 in propylene carbonate (PC), LiPF 6 in dimethyl carbonate (DMC) and LiClO 4 in DMC are investigated up to 1 molal with accuracy better than 3 °C when compared to experimental freezing point measurements. Finally, we find that the difference in freezing point depression between the propylene carbonate-based electrolyte and the dimethyl carbonate electrolytes originates from the difference in the solvent dielectric constant.

25 ENERGY STORAGE↗
DOE OSTI2020

Predicting Small Molecule Transfer Free Energies by Combining Molecular Dynamics Simulations and Deep Learning

Accurately predicting small molecule partitioning and hydrophobicity is critical in the drug discovery process. There are many heterogeneous chemical environments within a cell and entire human body. For example, drugs must be able to cross the hydrophobic cellular membrane to reach their intracellular targets, and hydrophobicity is an important driving force for drug–protein binding. Atomistic molecular dynamics (MD) simulations are routinely used to calculate free energies of small molecules binding to proteins, crossing lipid membranes, and solvation but are computationally expensive. Machine learning (ML) and empirical methods are also used throughout drug discovery but rely on experimental data, limiting the domain of applicability. We present atomistic MD simulations calculating 15,000 small molecule free energies of transfer from water to cyclohexane. This large data set is used to train ML models that predict the free energies of transfer. We show that a spatial graph neural network model achieves the highest accuracy, followed closely by a 3D-convolutional neural network, and shallow learning based on the chemical fingerprint is significantly less accurate. A mean absolute error of ~4 kJ/mol compared to the MD calculations was achieved for our best ML model. We also show that including data from the MD simulation improves the predictions, tests the transferability of each model to a diverse set of molecules, and show multitask learning improves the predictions. This work provides insight into the hydrophobicity of small molecules and ML cheminformatics modeling, and our data set will be useful for designing and testing future ML cheminformatics methods.

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

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

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

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗