Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOLVATION”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

IR spectral assignments for the hydrated excess proton in liquid water

The local environmental sensitivity of infrared (IR) spectroscopy to a hydrogen-bonding structure makes it a powerful tool for investigating the structure and dynamics of excess protons in water. Although of significant interest, the line broadening that results from the ultrafast evolution of different solvated proton-water structures makes the assignment of liquid-phase IR spectra a challenging task. In this paper, we apply a normal mode analysis using density functional theory of thousands of proton-water clusters taken from reactive molecular dynamics trajectories of the latest generation multistate empirical valence bond proton model (MS-EVB 3.2). These calculations are used to obtain a vibrational density of states and IR spectral density, which are decomposed on the basis of solvated proton structure and the frequency dependent mode character. Decompositions are presented on the basis of the proton sharing parameter δ, often used to distinguish Eigen and Zundel species, the stretch and bend character of the modes, the mode delocalization, and the vibrational mode symmetry. We find there is a wide distribution of vibrational frequencies spanning 1200-3000 cm -1 for every local proton configuration, with the region 2000-2600 cm -1 being mostly governed by the distorted Eigen-like configuration. Finally, we find a continuous red shift of the special-pair O∙∙∙H + ∙∙∙O stretching frequency, and an increase in the flanking water bending intensity with decreasing δ. Also, we find that the flanking water stretch mode of the Zundel-like species is strongly mixed with the flanking water bend, and the special pair proton oscillation band is strongly coupled with the bend modes of the central H 5 O 2 + moiety.

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

On combining the conductor-like screening model and optimally tuned range-separated hybrid density functionals

Range-separated hybrid functionals whose range-separation parameter γ has been nonempirically tuned to a particular molecule have been shown to yield frontier orbital energies and other properties in very good agreement with experiments. Yet, many cases, such as organic optoelectronic devices, require the description of molecules embedded in an environment. This can be done by combining the γ-tuning procedure with polarizable continuum models in general and the very versatile conductor-like screening model in particular. There are at least two different ways of performing this combination. The partially vertical γ-tuning employs equilibrium solvation throughout. The strictly vertical γ-tuning, on the other hand, employs nonequilibrium solvation to obtain ionization energies. Here, we compare ground-state and excited-state properties of several different molecules relevant to organic optoelectronics that were obtained using both of the two different tuning procedures. While there are significant differences in the ground-state properties, we see virtually no difference in the excited-state properties. Given these conclusions, we find that both tuning procedures have to be used in conjunction for the correct description of both ground-state and excited-state properties.

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

Molecular insights into ionic liquid/aqueous interface of phosphonium based phase-separable ionic liquids

Phosphonium based phase-separable ionic liquids (PSILs) are promising green solvents for dissolution of cellulose and lignin, a necessary step for conversion of biomass to fuels and chemicals. The knowledge of interfacial behavior of ionic liquid/solvent systems is critical for designing efficient separation processes. Here, molecular dynamics simulations are carried out for aqueous interface of tetraalkylphosphonium ionic liquids with chloride and acetate as anions to investigate IL miscibility with water. The transition zone from miscible to immiscible behavior was observed for alkyl chain lengths of 6 to 8. Emulsion phase was observed for [P8888]+ ion and multiple IL/water interface was observed for [P12121212]+ phosphonium cation. IL/water interface is observed to be enriched with solvated anions with phosphorous atom oriented towards the aqueous phase. The potential of mean force calculations suggest that as the alkyl chain of the phosphonium cations increase, so does the driving force for aggregation of cations in the aqueous phase. Due to the difference in the molecular shape and size, anions are better solvated by the water molecules, and these ions diffuse into the aqueous phase. Furthermore, the analysis of selected pair interactions provides insights into the nature of intermolecular forces and the role of the alkyl side chains on the interfacial properties.

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

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

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

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
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

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2025

Key Factors in Semi-Generic Coarse-Grained Modeling of Solid Polymer Electrolytes

Given the long length and time scales of interest and strong interactions present in solid polymer electrolytes, coarse-grained molecular models can be useful in understanding the molecular basis of their structure and transport properties. Rather than representing atoms individually, coarse-grained models use beads representing groups of atoms, making the system simpler and more efficient. Generic coarse-grained models, that are not built based on matching a particular atomistic system, can provide insight into the important physical considerations that apply across different chemical systems, though it can be unclear how to best match to experimental systems and capture relevant experimental behaviors with as simple of a model as possible. Here, we build on generic bead-spring models, but use stiff angle potentials, set different bead properties for different polymer types, and include additional ion parameters, creating a semi-generic coarse-grained model that can map more specifically to polymers and copolymers with different chain architectures, component glass transition temperatures, and ion solvation behaviors. We set most parameters with basic homopolymer data such as glass transition temperatures, Kuhn length, density, and dielectric constant, with further adjustment based on data from the polymer electrolyte system. We specifically model polystyrene-block-poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) with lithium triflate salt, considering the POEM to be made of a poly(methylmethacrylate) backbone with poly(ethylene oxide) side chains. Solvation of ions is accounted for by additional polymer-ion interactions of the form -S/r4, plus additional lithium-polymer Lennard-Jones interactions. We find these potentials have different effects and discuss strategies for setting these parameters.

Zhang, Yuanhao↗
DOE OSTI2025

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2021

An investigation of aqueous ammonium nitrate aerosols with soft X-ray spectroscopy

Aqueous aerosols are important in atmospheric chemistry, drug delivery, and human health. X-ray photoelectron spectroscopy (XPS), a surface-sensitive technique and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy which informs on the bulk, is deployed to study the solvation of ammonium nitrate (NH 4 NO 3 ). Aerosolised solutions of NH 4 NO 3 were introduced into a photoelectron spectrometer via an aerodynamic lens, and interrogated with soft X-ray photons, the resulting electrons were imaged via velocity map imaging. Density functional theory calculations were performed to compare with the measured binding energies of NH 4 + and NO 3 - solvated in water. The results reveal that the nitrate anion has a slight propensity for the surface, while both ammonium and nitrate ions are present in equal measure in the bulk.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

HCOO - aq Degradation in Droplets by OH aq in an Atmospheric Pressure Glow Discharge

Plasmas in contact with liquids can degrade organic molecules in solution, as reactive oxygen and nitrogen species produced in the plasma solvate into the liquid. Immersing small droplets (tens of microns in diameter) in the plasma can more rapidly activate the liquid compared to treating a large volume of liquid with a smaller surface-to-volume ratio. The interactions between a radio frequency glow discharge sustained in He/H 2 O and a water droplet containing formate (HCOO - aq ) immersed in and flowing through the plasma were modeled using a 0-dimensional global plasma chemistry model to investigate these activation processes. HCOO - aq interacts with OH aq , produced from solvation of OH from the gas phase. HCOO - aq concentrations were benchmarked with previously reported experimental measurements. The diameter of the droplet, initial HCOO - aq concentration, and gas flow rate affect only the HCOO - aq concentration and OH aq density, leaving the OH density in the gas phase unaffected. Power deposition and gas mixture (e.g., percentage of H 2 O) change both the gas and liquid phase chemistry. A general trend was observed – during the first portion of droplet exposure to the plasma, OH aq primarily consumes HCOO - aq . However, O2 - aq , a byproduct of HCOO - aq consumption, consumes OH aq once O2 - aq reaches a critically large density. Using HCOO - aq as a surrogate for OH aq sensitive contaminants, combinations of residence time, droplet diameter, water vapor density, and power will determine the optimum remediation strategy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Propagation of positive discharges in an air bubble having an embedded water droplet

Low temperature plasmas generated inside gas bubbles immersed in water is an effective method of rapidly transferring plasma generated reactive species to the water for applications in biomedicine, agriculture and environment. Reactive species are generally produced in the gas phase plasma and then solvate into the liquid. The large surface-to-volume ratio (SVR) of the bubble accelerates this process. In generating bubbles in water, aerosols and droplets are also contained within the bubble. These droplets also have a large SVR and so can be rapidly plasma activated. However, the presence of the droplets can also impact the propagation of the plasma in the bubble. In this paper, results are discussed from computational and experimental investigations of the formation and evolution of discharges in an air bubble immersed in water with an embedded water droplet. The computations were performed with a two-dimensional plasma hydrodynamics model. Experiments were performed with a quasi-2D bubble apparatus. In bubbles having a droplet, a plasma filament typically bridges from the powered electrode to the droplet, and then from the droplet to the bubble surface. A surface-hugging streamer also occurs on the inner bubble surface and on the surface of the droplet. Both surface streamers result in part from surface charge accumulation and can dominate the formation of reactive species that transport into the droplet. Increasing droplet conductivity suppresses propagation of the surface discharge and leads to a lower density of aqueous reactive species. Increasing conductivity of the surrounding water does not change the overall structure of the discharge but does slightly elevate the discharge intensity. Finally, the size and shape of the embedded droplet can significantly affect the formation and propagation of the streamer.

2D fluid modelling↗
DOE OSTI2018

Ion Correlation Effects in Salt-Doped Block Copolymers

Here, we apply classical density functional theory to study how salt changes the microphase morphology of diblock copolymers. Polymers are freely jointed and one monomer type favorably interacts with ions, to account for the selective solvation that arises from different dielectric constants of the microphases. By including correlations from liquid state theory of an unbound reference fluid, the theory can treat chain behavior, microphase separation, ion correlations, and preferential solvation, at same coarse grained level. We show good agreement with molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Importance of nuclear quantum effects on the hydration of chloride ion

Path-integral ab initio molecular dynamics (PI-AIMD) calculations have been employed to probe the nature of chloride ion solvation in aqueous solution. Nuclear quantum effects (NQEs) are shown to weaken hydrogen bonding between the chloride anion and the solvation shell of water molecules. As a consequence, the disruptive effect of the anion on the solvent water structure is significantly reduced compared to what is found in the absence of NQEs. Additionally, the chloride hydration structure obtained from PI-AIMD agrees well with information extracted from neutron scattering data. In particular, the observed satellite peak in the hydrogen-chloride-hydrogen triple angular distribution serves as a clear signature of NQEs. The present results suggest that NQEs are likely to play a crucial role in determining the structure of saline solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anti-Coulomb ion-ion interactions: A theoretical and computational study

The free energy of ion solvation can be decomposed into enthalpic and entropic contributions. This helps one to understand the connection between the dielectric properties and the underlying forces. We present a simple linear-response model of screened charge interactions that provides an alternative understanding of solvation barriers. Moreover, it explains the “anti-Coulomb” interactions (attraction between like-charged ions and repulsion between opposite-charged ions) observed in both simulations and experiments. We show that this is a universal behavior associated to the nonlocal response function of any dielectric or metallic systems. Published by the American Physical Society 2024

Wills, Alec (ORCID:0000000280184692)↗
DOE OSTI2017

Solvent exchange in a metal–organic framework single crystal monitored by dynamic in situ X-ray diffraction

Understanding the processes by which porous solid-state materials adsorb and release guest molecules would represent a significant step towards developing rational design principles for functional porous materials. To elucidate the process of liquid exchange in these materials, dynamic in situ X-ray diffraction techniques have been developed which utilize liquid-phase chemical stimuli. Using these time-resolved diffraction techniques, the ethanol solvation process in a flexible metal–organic framework [Co(AIP)(bpy) 0.5 (H 2 O)]·2H 2 O was examined. The measurements provide important insight into the nature of the chemical transformation in this system including the presence of a previously unreported neat ethanol solvate structure.

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

Probing liquid/liquid interfaces at and away from equilibrium using vibrational sum frequency generation

The self-assembly and transport of species at oil/aqueous interfaces represent key mechanistic steps in processes spanning from liquid extraction to soft-matter neuromorphic electrical devices. While technologically and fundamentally important, there is a notable lack of chemical understanding regarding these buried interfaces due to challenges in differentiating species located at the surface from the nearby bulk phases. In this work, we will describe some of the challenges in probing chemistry at these liquid/liquid interfaces and detail the approaches we have developed using vibrational sum frequency generation to enable studies of these complex interfaces at and away from equilibrium. From these spectroscopic developments, new insight into solvation, ion pairing, aggregation, complexation, and transport will be discussed in the context of solvent extraction chemistry. Notably, we show how the presence of both oil and aqueous phases provides an environment to drive fundamentally different phenomena from what is observed at model air/aqueous interfaces. This difference arises from key differences in solvation and bulk populations that feedback onto the surface to drive emergent and often functional molecular assemblies.

Doughty, Benjamin L↗