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.
Solvation of Molecular Excited-States: (Non)-Equilibrium Transition Energies in Donor-Acceptor
Explore the source record for details and available documents.
First-Principles-Based Computational Study on Nucleation and Growth Mechanisms of U on Mo(110) Surface Solvated in an Eutectic LiCl-KCl Molten Salt
In nuclear engineering and industry electrochemical methods have been intensively employed in pyro-processing technology to recycle used nuclear fuels for substantial alleviation of the concern with the safe treatment of high-level radioactive nuclear wastes. Pyro-processing is based on sequential electrochemical reactions at high temperature (typically above 700 K at ambient pressure) to enhance the efficiency for decoupling nuclear fuels from unfavorable fission fragments. Using a molten salt in electrorefining process, at which actinides can be selectively deposited in solid cathode from ionic states, gives several advantages over an aqueous solution: higher thermal conductivity (it is liquid phase at high temperature, molten salt plays a role as medium at high temperature over 773 K in electrochemical reaction, lower kinetic barrier in ionic diffusion mobility and better stability against undesired side reactions with radioactive materials. Therefore, the volume of high-level nuclear wastes to be, otherwise, disposed into secured sites is dramatically reduced. Most of all it conforms to nuclear non-proliferation treaty. Electrochemical reactions in a molten salt medium are, however, much more complicated than in aqueous solution rendering the accurate understanding of thermodynamic and kinetic properties formidably challenging. It is represented by extensive researches in experiments and computations to identify underlying mechanisms on how U and other fission fragments behave in the electrorefining system. Regardless of the researches there is a general agreement that the initial stage of the nucleation and growth process of the U deposit is critical for determining the properties of final product and the electrorefining efficiency. In spite of the focused studies it still remains unclear. Theoretical and computational researches have been often carried out with various empirical parameter for interatomic potential energy, which are unknown variables until numerical fits to experimental measurements or just assumption is made. Thus, empirical simulations sometimes predicted outcomes are inconsistent depending on materials of interest with experimental measurements. This may be largely because the accuracy of an empirical interatomic potential energy models critically depends on the input parameters and thus, can be limited in transferability of the results. Furthermore, it may be devoid of fundamental physics or chemistry. This situation requires special caution when an electrorefining system is studied to recycle spent nuclear fuels because of the extreme difficulty in reliable acquisition of thermodynamic and kinetic data from direct experimental measurements. In this presentation, we extensively utilized first principles density functional theory (DFT) calculations and ab-initio molecular dynamic (AIMD) simulations to figure out underlying mechanisms of the initial stage of nucleation and growth behaviors of U on Mo(110) surface in a eutectic LiCl-KCl molten salt at T = 773 K. To rigorously investigate the effect of the molten salt electrolyte we compared the results with those calculated under vacuum condition.
Unusual dynamic properties of water near the ice-binding plane of hyperactive antifreeze protein
The dynamical properties of solvation water of hyperactive antifreeze protein from Choristoneura fumiferana (CfAFP) are analyzed and discussed in context of its antifreeze activity. The protein comprises of three well-defined planes and one of them binds to the surface of ice. The dynamical properties of solvation water around each of these planes were analyzed separately; the results are compared with the dynamical properties of solvation water of ice around its two crystallographic planes: basal and prism. Three main conclusions are inferred from our investigations. The first one is that the solvation shell of CfAFP does not seem to be particularly far-ranged, at least not beyond what is usually observed for proteins that do not interact with ice. Therefore, it does not appear to us that the antifreeze activity is enhanced by a long-ranged retardation of water mobility. Also the correlation between the collective mobility of water and the collective mobility of protein atoms highly resembles the one measured for the protein that does not interact with ice. Our second conclusion is that the dynamical properties of solvation water of CfAFP are non-uniform. The dynamics of solvation water of ice-binding plane is, in some respects, different from the dynamics of solvation water of the two remaining planes. The feature that distinguishes the dynamics of solvation water of the three planes is the activation energy of diffusion process. The third conclusion is that—from the three analyzed solvation shells of CfAFP—the dynamical properties of solvation water of the ice-binding plane resemble the most the properties of solvation water of ice; note, however, that these properties still clearly differ from the dynamic properties of solvation water of ice.
Water Lone Pair Delocalization in Classical and Quantum Descriptions of the Hydration of Model Ions
Understanding the nature of ionic hydration at a fundamental level has eluded scientists despite intense interest for nearly a century. In particular, the microscopic origins of the asymmetry of ion solvation thermodynamics with respect to the sign of the ionic charge remains a mystery. Here, we determine the response of accurate quantum mechanical water models to strong nanoscale solvation forces arising from excluded volumes and ionic electrostatic fields. This is compared to the predictions of two important limiting classes of classical models of water with fixed point changes, differing in their treatment of “lone pair” electrons. Using the quantum water model as our standard of accuracy, we find that a single fixed classical treatment of lone pair electrons cannot accurately describe solvation of both apolar and cationic solutes, emphasizing the need for a more flexible description of local electronic effects in solvation processes. However, we explicitly show that all water models studied respond to weak long-ranged electrostatic perturbations in a manner that follows macroscopic dielectric continuum models, as would be expected. Here, we emphasize the importance of these findings in the context of realistic ion models, using density functional theory and empirical models, and discuss the implications of our results for quantitatively accurate reduced descriptions of solvation in dielectric media.
Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning
Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.
Liquid-Phase Effects on Adsorption Processes in Heterogeneous Catalysis
Aqueous solvation free energies of adsorption have recently been measured for phenol adsorption on Pt(111). Endergonic solvent effects of ~1 eV suggest solvents dramatically influence a metal catalyst’s activity with significant implications for the catalyst design. However, measurements are indirect and involve adsorption isotherm models, which potentially reduces the reliability of the extracted energy values. Computational, implicit solvation models predict exergonic solvation effects for phenol adsorption, failing to agree with measurements even qualitatively. In this study, an explicit, hybrid quantum mechanical/molecular mechanical approach for computing solvation free energies of adsorption is developed, solvation free energies of phenol adsorption are computed, and experimental data for solvation free energies of phenol adsorption are reanalyzed using multiple adsorption isotherm models. Explicit solvation calculations predict an endergonic solvation free energy for phenol adsorption that agrees well with measurements to within the experimental and force field uncertainties. Computed adsorption free energies of solvation of carbon monoxide, ethylene glycol, benzene, and phenol over the (111) facet of Pt and Cu suggest that liquid water destabilizes all adsorbed species, with the largest impact on the largest adsorbates.
Localized High Concentration Electrolytes for High Voltage Lithium-Metal Batteries: Correlation between the Electrolyte Composition and Its Reductive/Oxidative Stability
We demonstrate a first-principles screening methodology as an effective tool to explore electrolyte formulations for the new generation of high energy density rechargeable batteries. We study the liquid structure and electronic properties in dilute electrolytes, high concentration electrolytes (HCE), and localized high concentration electrolytes (LHCE), with focus on electrolyte formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and bis(2,2,2-trifluoroethyl) ether (BTFE) as a diluent. We describe the solvation complexes in the dilute electrolyte and explore structural changes triggered by the increase in lithium salt concentration for HCEs and the diluent effects in LHCEs. In HCE formulations, there is a 4-fold coordination environment of lithium-ions as in the dilute electrolyte, but the number of lithium-ion interactions with O atoms from FSI- anions dominates. In these solutions, the ability of the FSI- anions to interact with multiple lithium-ions allows complex 3D network formation and influences the reductive/oxidative behavior of the electrolyte. Interestingly, in LHCEs, the BTFE diluent molecules do not change the 3D solution structure when diluting the HCE formulation from 5.49 to 3.83 M. However, there is a composition threshold where the structural and electronic behavior may change. We show that diluting the HCE electrolyte with BTFE down to 1.77 M breaks the three-dimensional solution structure into an island-like solvation complex. We relate these structural changes to the electronic properties of the electrolytes finding a causal relationship between the reductive/oxidative behavior and the lithium-oxygen interaction mechanisms in the solvated complexes. The coordination with lithium-ions lowers the electrolyte LUMO and HOMO levels: the higher is the number of interactions with lithium-ions, the more likely the solvent molecule, FSI- anion, or diluent molecule is to be reduced and the less likely it is to become oxidized. The evolution of the solvated ion structure in HCE and LHCE suggests a close connection to a corresponding change in the lithium-ion transport mechanisms for these electrolytes.
Coordination-Dependent Chemical Reactivity of TFSI Anions at a Mg Metal Interface
Charge transfer across the electrode–electrolyte interface is a highly complex and convoluted process involving diverse solvated species with varying structures and compositions. Despite recent advances in in situ and operando interfacial analysis, molecular specific reactivity of solvated species is inaccessible due to a lack of precise control over the interfacial constituents and/or an unclear understanding of their spectroscopic fingerprints. However, such molecular-specific understanding is critical to the rational design of energy-efficient solid–electrolyte interphase layers. We have employed ion soft landing, a versatile and highly controlled method, to prepare well-defined interfaces assembled with selected ions, either as solvated species or as bare ions, with distinguishing molecular precision. Equipped with precise control over interfacial composition, we employed in situ multimodal spectroscopic characterization to unravel the molecular specific reactivity of Mg solvated species comprising (i.e., bis(trifluoromethanesulfonyl)imide, TFSI – ) anions and solvent molecules (i.e., dimethoxyethane, DME/G1) on a Mg metal surface relevant to multivalent Mg batteries. In situ multimodal spectroscopic characterization revealed higher reactivity of the undercoordinated solvated species [Mg-TFSI-G1] + compared to the fully coordinated [Mg-TFSI-(G1) 2 ] + species or even the bare TFSI – . These results were corroborated by the computed reaction pathways and energy barriers for decomposition of the TFSI – within Mg solvated species relative to bare TFSI – . Finally, we evaluated the TFSI reactivity under electrochemical conditions using Mg(TFSI) 2 –DME-based phase-separated electrolytes representing different solvated constituents. Based on our multimodal study, we report a detailed understanding of TFSI – decomposition processes as part of coordinated solvated species at a Mg-metal anode that will aid the rational design of improved sustainable electrochemical energy technologies.
Computational Investigation of Aqueous Phase Effects on the Dehydrogenation and Dehydroxylation of Polyols over Pt(111)
Prediction of solvation effects on the kinetics of elementary reactions occurring at metal–water interfaces is of high importance for the rational design of catalysts for the biomass and electrocatalysis communities. A lack of knowledge of the reliability of various computational solvation schemes for processes on metal surfaces is currently a limiting factor. Using a multilevel quantum mechanical/molecular mechanical (QM/MM) description of the potential energy surface, we determined characteristic time and length scales for typical free-energy perturbation (FEP) calculations of bond cleavages in ethylene glycol, a sugar surrogate molecule, over Pt(111). Our approach is based on our explicit solvation model for metal surfaces and the repetition of FEP calculations to estimate confidence intervals. Results indicate that aqueous phase effects on the free energies of elementary processes can be determined with 95% confidence intervals from limited configuration space sampling and the fixed charge approximation used in the QM/MM-FEP methodology of smaller 0.1 eV. Next, we computed the initial O–H, C–H, and C–OH bond cleavages in ethylene glycol over Pt(111) in liquid water utilizing two different metal–water interaction potentials. Our calculations predict that aqueous phase effects are small (<0.1 eV) for the C–H bond cleavage and the activation barrier of the C–OH bond cleavage. In contrast, solvation effects are large (>0.35 eV) for the O–H bond cleavage and the reaction free energy of the C–OH bond scission. While the choice of a different Pt–water force field can lead to differences in predicted solvation effects of up to 0.2 eV, the differences are usually smaller (<0.1 eV), and the trends are always the same. In contrast, implicit solvation methods appear to currently not be able to reliably describe solvation effects originating from hydrogen bonding for metal surfaces even qualitatively.
Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water
Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.
Assessing the impacts of dynamic soft-templates innate to switchable ionic liquids on nanoparticulate green rust crystalline structures
Explore the source record for details and available documents.
Rapid water dynamics structures the OH-stretching spectra of solitary water in ionic liquids and dipolar solvents
In a recent study [J. Phys. Chem. B 126, 4584–4598 (2022)], we have used infrared spectroscopy to investigate the solvation and dynamics of solitary water in ionic liquids and dipolar solvents. Complex shapes observed for water OH-stretching bands, common to all high-polarity solvents, were assigned to water in several solvation states. Here, in the present study, classical molecular dynamics simulations of a single water molecule in four ionic liquids and three dipolar solvents were used to test and refine this interpretation. Consistent with past assignments, simulations show solitary water usually donates two hydrogen bonds to distinct solvent molecules. Such symmetrically solvated water produces the primary pair of peaks identified in the OH spectra of water in nearly all solvents. We had further proposed that additional features flanking this main peak are due to asymmetric solvation states, states in which only one OH group makes a hydrogen bond to solvent. Such states were found in significant concentrations in all of the systems simulated. Simulations of the OH stretching spectra using a semiclassical description and the vibrational map developed by Auer and Skinner [J. Chem. Phys. 128, 224511–224512 (2008)] provided semi-quantitative agreement with experiment. Analysis of species-specific spectra confirmed assignment of the additional features in the experimental spectra to asymmetrically solvated water. The simulations also showed that rapid water motions cause a marked motional narrowing compared with the inhomogeneous limit. This narrowing is largely responsible for making the additional features due to minority solvation states manifest in the spectra.
Solvent-Induced Shift of Spectral Lines in Polar–Polarizable Solvents
Solvent-induced shift of optical transition lines is traditionally described by the Lippert- McRae equation given in terms of the Onsager theory for dipole solvation. It splits the overall shift into the equilibrium solvation by induced dipoles and the reaction field by the permanent dipoles in equilibrium with the chromophore in the ground state. Here we have reconsidered this classical problem from the perspective of microscopic solvation theories. A microscopic solvation functional is derived and continuum solvation is consistently introduced by taking the limit of zero wavevector in the reciprocal-space solvation susceptibility functions. We show that the phenomenological expression for the reaction field of permanent dipoles in the Lippert-McRae equation is not consistent with the microscopic theory. The main deficiency of the Lippert- McRae equation equation is the use of additivity of the response by permanent and induced dipoles of the liquid. An alternative closed-form equation for the spectral shift is derived. Its continuum limit allows a new, non-additive functionality for the solvent-induced shift in terms of the high-frequency and static dielectric constants. Finally, the main qualitative outcome of the theory is a significantly weaker dependence of the spectral shift on the polarizability of the solvent than predicted by the Lippert-McRae formula.
Microsolvation Structures of Protonated Glycine and L -Alanine
The IR predissociation spectra of microsolvated glycine and l-alanine, GlyH + (H 2 O) n and AlaH + (H 2 O) n , n = 1–6, are presented. The assignments of the solvation structures are aided by H 2 O/D 2 O substitution, IR-IR double resonance spectroscopy, and computational efforts. The analysis reveals the water–amino acid as well as the water–water interactions, and the subtle effects of the methyl side chain in L-alanine on the solvation motif are also highlighted. The bare amino acids exhibit an intramolecular hydrogen bond between the protonated amine and carboxyl terminals. In the n = 1–2 clusters, the water molecules preferentially solvate the protonated amine group, and we observed differences in the relative isomer stabilities in the two amino acids due to electron donation from the methyl weakening the intramolecular hydrogen bond. The structures in the n = 3 clusters show a further preference for solvation of the carboxyl group in l-alanine. For n = 4–6 clusters, the solvation structure of the two amino acids is remarkably similar, with one dominant isomer present in each cluster size. The first solvation shell is completed at n = 4, evidenced by a lack of free NH and OH stretches on the amino acid, as well as the first observation of H 2 O–H 2 O interactions in the spectra of n = 5. Finally, we note that calculations at the density functional theory (DFT) level show excellent agreement with the experiment for the smaller clusters. Furthermore, when water–water interactions compete with water–amino acid interactions in the larger clusters, DFT results show greater disagreement with experiment when compared to MP2 results.
Ligand-channel-enabled ultrafast Li-ion conduction
Li-ion batteries (LIBs) for electric vehicles and aviation demand high energy density, fast charging and a wide operating temperature range, which are virtually impossible because they require electrolytes to simultaneously have high ionic conductivity, low solvation energy and low melting point and form an anion-derived inorganic interphase. We report guidelines for designing such electrolytes by using small-sized solvents with low solvation energy. The tiny solvent in the secondary solvation sheath pulls out the Li + in the primary solvation sheath to form a fast ion-conduction ligand channel to enhance Li + transport, while the small-sized solvent with low solvation energy also allows the anion to enter the first Li + solvation shell to form an inorganic-rich interphase. The electrolyte-design concept is demonstrated by using fluoroacetonitrile (FAN) solvent. The electrolyte of 1.3 M lithium bis(fluorosulfonyl)imide (LiFSI) in FAN exhibits ultrahigh ionic conductivity of 40.3 mS cm -1 at 25 °C and 11.9 mS cm -1 even at -70 °C, thus enabling 4.5-V graphite||LiNi 0.8 Mn 0.1 Co 0.1 O 2 pouch cells (1.2 Ah, 2.85 mAh cm -2 ) to achieve high reversibility (0.62 Ah) when the cells are charged and discharged even at -65 °C. The electrolyte with small-sized solvents enables LIBs to simultaneously achieve high energy density, fast charging and a wide operating temperature range, which is unattainable for the current electrolyte design but is highly desired for extreme LIBs. This mechanism is generalizable and can be expanded to other metal-ion battery electrolytes.
Exploring the impact of ions on oxygen K-edge X-ray absorption spectroscopy in NaCl solution using the GW-Bethe-Salpeter-equation approach
X-ray absorption spectroscopy (XAS) is a powerful experimental tool to probe the local structure in materials with the core hole excitations. Here, the oxygen K-edge XAS spectra of the NaCl solution and pure water are computed by using a recently developed GW-Bethe-Salpeter equation approach, based on configurations modeled by path-integral molecular dynamics with the deep-learning technique. The neural network is trained on ab initio data obtained with strongly constrained and appropriately normed density functional theory. The observed changes in the XAS features of the NaCl solution, compared to those of pure water, are in good agreement between experimental and theoretical results. We provided detailed explanations for these spectral changes that occur when NaCl is solvated in pure water. Specifically, the presence of solvating ion pairs leads to localization of electron-hole excitons. As a result, our theoretical XAS results support the theory that the effects of the solvating ions on the H-bond network are mainly confined within the first hydration shell of ions, however beyond the shell the arrangement of water molecules remains to be comparable to that observed in pure water.
Phenolic Polymer Interactions with Water and Ethylene Glycol Solvents
Interactions between pre-cured phenolic polymer chains and a solvent have a significant impact on the structure and properties of the final post-cured phenolic resin. Developing an understanding of the nature of these interactions is important and will aid in the selection of the proper solvent that will lead to the desired final product. Here, we investigate the role of the phenolic chain structure and the solvent type on the overall solvation performance of the system through ab initio techniques and molecular dynamics simulations. Two types of solvents are considered: ethylene glycol (EGL) and H2O. Three phenolic chain structures are considered, including two novolac-type chains with either an ortho-ortho (OON) or an ortho-para (OPN) backbone network and a resole-type (RES) chain with an ortho-ortho network. Each system is characterized through a structural analysis of the solvation shell and the hydrogen-bonding environment as well as through a quantification of the solvation free energy along with partitioned interaction energies between specific molecular species. The combination of simulations and the analyses indicate that EGL provides a higher solvation free energy than H2O due to more energetically favorable hydrophilic interactions as well as favorable hydrophobic interactions between CH element groups. In addition, the phenolic chain structure significantly affects the solvation performance, with OON having limited intermolecular hydrogen-bond formations, while OPN and RES interact more favorably with the solvent molecules. The results suggest that a resole-type phenolic chain with an ortho-para network should have the best solvation performance in EGL, H2O, and other similar solvents.