Search NASA⌕ Search

SEARCH · Search NASA

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

Balancing solvation: stabilizing lithium metal batteries via optimized cosolvents for ionic-liquid electrolytes

In this study, we examined three cosolvents with distinct solvation capabilities for ionic-liquid electrolytes based on 1-methyl-1-propyl pyrrolidinium bis(fluorosulfonyl)imide (Py13FSI). We demonstrate that 1,1,1-trifluoro-2-(2-(2-(2,2,2-trifluoroethoxy)ethoxy)ethoxy)ethane (FDG) notably enhances the cycle life of Py13FSI-based electrolytes, outperforming 1,1,2,2-tetrafluoroethyl 2,2,3,3-tetrafluoropropylether (TTE) and diglyme (DG). Electrochemical and surface analyses showed that this improvement could be attributed to the formation of a favorable cathode interphase, promoting efficient Li+ transport with reduced overpotential. Spectroscopic techniques (FTIR, Raman, and NMR spectroscopy) and molecular dynamics simulations revealed that cosolvents with varying solvation abilities can influence the solvation structures in Py13FSI-based electrolytes. The mild solvating strength and lithium stability of FDG are key contributors to its effectiveness. Conversely, DG, a strong solvating solvent, destabilized the Py13FSI-DG electrolyte at the lithium metal anode, while TTE, a non-solvating solvent, failed to enhance lithium transport or form a stable cathode interphase. Our findings highlight that balanced solvation exerted by the cosolvents is critical for forming a stable electrolyte–cathode interface, potentially through FSI decomposition. This study offers valuable insights into the development of durable ionic-liquid electrolytes, emphasizing the importance of selecting cosolvents with optimal solvation properties.

Li, Xinlin [Argonne National Laboratory (ANL), Arg↗
DOE OSTI2026

Deciphering the Dynamic Balance Between Solvation Strength and Polysulfides Reaction Heterogeneity in Practical Lithium‐Sulfur Batteries

Achieving stable interfacial chemistry in lithium–sulfur batteries under practical conditions remains a key barrier to commercialization. Here, we demonstrate that interfacial dynamics can be effectively regulated by coupling solvation-power control with intrinsic heterogeneity of sulfur redox chemistry through the introduction of a weakly solvating fluorinated cosolvent, LIB 1200ET (1200ET). Compared with conventional fluorinated ethers, 1200ET efficiently shifts Li + solvation environment toward a more non-coordinated configuration at low volume fractions, enabling substantial solvation modulation without significantly impairing sulfur redox kinetics. This solvation transition weakens Li + –solvent interactions while strengthening Li + –anion and Li + –lithium polysulfide (LPS) coordination, suppressing LPS solubility and promoting reconstruction of solid–electrolyte interphase (SEI). Regulated LPS chemistry, together with 1200ET, leads to formation of a S 4+ -rich, LiF-reinforced SEI with enhanced ionic conductivity and mechanical robustness. Spatially resolved sulfur K-edge X-ray absorption spectroscopy on pouch cells reveals pronounced current-density-dependent chemical heterogeneity, distinguishing kinetically dominated and solvation-controlled regions. Under practical conditions (3.7 mg cm −2 sulfur loading, E/S = 6 µL mg −1 ), a single-layer pouch cell delivers 527 mAh g −1 over 200 cycles at C/3, while an Ah-level multilayer pouch cell achieves an energy density of 358 Wh kg −1 . These results establish non-coordinating cosolvent-driven solvation engineering as a scalable strategy for practical Li–S batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DFT-based investigation of solvation of Nd(III)/Yb(III) cations in room-temperature ionic-liquid

We report a DFT study on the solvation of bare rare-earth (RE) cations [Nd(III) and Yb(III)] in the room-temperature ionic-liquid (RTIL) based on 1,3-dimethyl-imidazolium + , PF 6 - , [MMI][PF 6 ]. In the present article, we aim to investigate whether the chosen RTIL can preferentially solvate any of the RE cations under study. First, the RTIL binding to the RE cation is studied in the gas phase, and then explored with continuum solvation model employing RTIL parameters. First, we considered the formation of the first solvation shell, where the PF 6 - anions directly bond with the Ln(III) cation. We then further surrounded the first solvation shell with [MMI]+ cations to form the second solvation shell. The stability of PF 6 - coordinated complexes is studied by calculating the binding energy values. The number of PF 6 - anions for both Ln(III) cations starts with three ligands and goes until the maximum binding energy in the RTIL medium is reached. The coordination number and binding modes of PF 6 - depend on the size of the RE cation. Comparing the binding energy values between the RE cations for the first solvation shell complexes, it is noticed that Yb-complexes have relatively greater binding energy than Nd-complexes. We have further performed a Boltzmann population analysis to find out the most dominant conformer among each category of PF 6 - coordinated complexes. A study using Natural Population Analysis shows that after complexation with PF 6 - the charge on the Yb-centers lower to a greater extent than Nd-centers, implying stronger bonding between Yb and PF 6 - , which is in harmony with the binding energy analysis. Overall, the computational investigation provides an in-depth understanding of the complexation phenomena of RE cations in RTIL, which can be useful for the practical separation of RE cations using RTIL.

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

Machine Learning-Guided Approach for Studying Solvation Environments

Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, in this work, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1–, and 2–, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.

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

Correlating Solvation Free Energy to Electrolyte Properties for Lithium Metal Batteries

The electrolyte plays a critical role in lithium metal batteries. In particular, ion solvation profoundly impacts key electrolyte properties and battery performance. Here, in this study, we systematically investigate solvation-property relationships in a series of electrolytes with different solvent-diluent ratios. Through potentiometric techniques that measure the relative solvation free energies of electrolytes, we find that weaker solvation correlates with larger ion clusters, lower ionic conductivity and diffusion coefficient, and superior electrochemical stability. Weaker solvation leads to the formation of a small number of Li clusters with large hydrodynamic radii, which lowers the Li + diffusivity and ionic conductivity of the electrolyte. Concurrently, weaker solvation leads to improved electrochemical stability at both the cathode and anode interfaces. Understanding these solvation-property relationships and trade-offs is important to designing electrolytes for optimized lithium metal battery performance.

25 ENERGY STORAGE↗
DOE OSTI2020

Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches

Determining the influence of the solvent on electrochemical reaction energetics is a central challenge in our understanding of electrochemical interfaces. To date, it is unclear how well existing methods predict solvation energies at solid/liquid interfaces, since they cannot be assessed experimentally. Ab initio molecular dynamics (AIMD) simulations present a physically highly accurate, but also a very costly approach. In this work, we employ extensive AIMD simulations to benchmark solvation at charge-neutral metal/water interfaces against commonly applied continuum solvent models. We consider a variety of adsorbates including *CO, *CHO, *COH, *OCCHO, *OH, and *OOH on Cu, Au, and Pt facets solvated by water. The surfaces and adsorbates considered are relevant, among other reactions, to electrochemical CO 2 reduction and the oxygen redox reactions. We determine directional hydrogen bonds and steric water competition to be critical for a correct description of solvation at the metal/water interfaces. As a consequence, we find that the most frequently applied continuum solvation methods, which do not yet capture these properties, do not presently provide more accurate energetics over simulations in vacuum. We find most of the computed benchmark solvation energies to linearly scale with hydrogen bonding or competitive water adsorption, which strongly differ across surfaces. Furthermore, we determine solvation energies of adsorbates to be non-transferable between metal surfaces, in contrast to standard practice.

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

Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide

Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.

Chemical compounds↗
DOE OSTI2026

Isolating solvent–solute hydrogen bonding interactions via 2D IR solvation shell spectroscopy

The solvation shell around a solute is a fundamental feature of liquid-phase solutions, determining the behavior and properties of both the solute and the overall solution. Direct experimental measurements of the solvation shell properties are challenging due to the strong signals generated from the bulk solvent, which overwhelm the small contribution of the solvation shell. Here, we use ultrafast two dimensional infrared (2D IR) spectroscopy and intermolecular cross-peaks to isolate the IR absorption spectrum of methanol molecules in the solvation shell surrounding the solute N-methylacetamide. We demonstrate that the intermolecular coupling between the solvent and solute vibrations is indirectly mediated by a low-frequency hydrogen-bonding mode, suggesting an important mechanism for anharmonic coupling induced by hydrogen bonds. From the relative frequency shifts and cross-peak anisotropy, we find that methanol molecules surrounding N-methylacetamide form stronger and distinctly oriented hydrogen bonds than those in the bulk solvent. Here, we also compare these results with the solvent spectra of the solute N,N-dimethylacetamide to investigate how solute structural changes alter the solvation shell and the contribution of N–H hydrogen bond donation. Our results are supported by molecular dynamics simulations, which provide detailed insights into the hydrogen-bonding distributions. Through these results, we demonstrate 2D solvation shell spectroscopy to be a valuable method for investigating solvation structures and dynamics without interference from the bulk solvent.

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 OSTI2015

Solvation effects on the band edge positions of photocatalysts from first principles

The band edge positions of photocatalysts relative to the redox potentials of water play an important role in determining the efficiency of photoelectrochemical cells. These band positions depend on the structure of the solid–liquid interface, but direct ab initio molecular dynamics calculations of these interfaces, while expected to be accurate, are too computationally demanding for high-throughput materials screening. Thus rapid theoretical screening of new photocatalyst materials requires simplified continuum solvation models that are suitable for treating solid–liquid interfaces. In this paper, we evaluate the accuracy of the recently developed CANDLE and SaLSA continuum solvation models for predicting solvation effects on the band positions of several well-studied surfaces [Si(111), TiO 2 (110), IrO 2 (110) and WO 3 (001)] in water. We find that the solvation effects vary considerably, ranging from <0.5 eV for hydrophobic surfaces, 0.5–1 eV for many hydrophilic oxide surfaces, to ~2 eV for oxygen-deficient surfaces. The solvation model predictions are in excellent agreement (within ~0.1 eV) with ab initio molecular dynamics results where available, and in good agreement (within ~0.2–0.3 eV) with experimental measurements. We also predict the energetics for surface oxygen vacancies and their effect on the band positions of the hydrated WO 3 (001) surface, leading to an explanation for why the solvation shift observed experimentally is substantially larger than predicted for the ideal surface. Finally, based on the correlation between solvation shift and the type of surface and solvent, we suggest approaches to engineer the band positions of surfaces in aqueous and non-aqueous solutions.

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

Impact of solvation on the electronic resonances in uracil

Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.

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

Metadynamics investigation of lanthanide solvation free energy landscapes and insights into separations energetics

Lanthanide ion solvation chemistry in nonaqueous phases is key to understanding and developing effective separation processes for these critical materials. Due to the complexity and inherent disorder of the solution phase, a comprehensive picture of the solvated metal ion is often difficult to generate solely from conventional spectroscopic approaches and electronic structure calculations, particularly in the extractant phase. In this work, we use classical molecular dynamics (MD) simulation with an advanced sampling technique, metadynamics, supplemented by experimental spectroscopy and speciation analysis, to measure lanthanide solvation free energy landscapes. We define coordination-based collective variables to probe the entire range of solvation configurations in the organic phase of lanthanum (La), europium (Eu), and lutetium (Lu) nitrate salts bound with a commonly used extractant, N,N′-dimethyl, N,N′-dioctylhexylethoxymalonamide (DMDOHEMA). The known lanthanide extraction trend of La ≈ Eu > Lu is readily explained by the measured free energy surfaces, which show consistent DMDOHEMA coordination from La to Eu, followed by loss of DMDOHEMA coordination from Eu to Lu. These simulations suggest how ligand crowding at the metal center can control selectivity, in this case resulting in the opposite extraction trend as observed with other conventional extractants, where the enthalpic contribution from increasing lanthanide charge density across the series dominates the extraction energetics. We also find that the presence of inner-sphere water, verified by time-resolved fluorescence, diversifies the accessible solvation structures. As a result, understanding solvation requires consideration of an entire thermodynamic ensemble, rather than the single dominant lowest-energy structure, as is often considered out of necessity in interpretation of spectroscopic data or in electronic structure-based ligand design approaches. In general, we demonstrate how metadynamics uniquely enables investigation of complex, multidimensional solvation energetic landscapes, and how it can explain selectivity trends where extraction is controlled by more complex mechanisms than simple charge density-based selectivity.

Wang, Xiaoyu↗
DOE OSTI2021

Solvation-protection-enabled high-voltage electrolyte for lithium metal batteries

To facilitate the practical application of lithium metal batteries (LMBs), stable interfaces between the electrolyte and the lithium metal must be achieved. Herein, we introduce a solvation protection strategy for designing a functional electrolyte for high-voltage LMBs. Fluoroethylene carbonate (FEC) was introduced as a solvation protection solvent for the difluoroethylene carbonate (DFEC)/trifluoroethyl methyl carbonate (FEMC) electrolyte system to enable the cycling of lithium metal anode. The addition of FEC alters the structures of lithium complexes in solution because of its relatively high solvating power. Through the precise control of the solvation number (> 1) of fluorinated cyclic carbonate (i.e., FEC:DFEC > critical ratio), lithium complexes with Li + solvated solely by FEMC, which decompose on the lithium surface to form detrimental by-products, can be effectively eliminated. Here, the new ternary FEC/DFEC/FEMC system not only maintains the beneficial effect of DFEC in forming a robust solid-electrolyte interphase on the lithium anode, but also confers outstanding anodic stability provided by FEMC, while eliminating detrimental FEMC decomposition through the solvation protection effect of FEC. Clearly, this ternary system outperforms the FEC/FEMC and DFEC/FEMC binary systems in facilitating the stable cycling of LMBs.

25 ENERGY STORAGE↗
DOE OSTI2021

Investigating the Variation in Solvation Interactions of Choline Chloride-Based Deep Eutectic Solvents Formed Using Different Hydrogen Bond Donors

Choline chloride-based deep eutectic solvents (DESs) have gained widespread recognition as green solvents in catalysis and separation science as their designer properties permit solvation properties to be modulated by the choice of the hydrogen bond donor (HBD) as well as its molar ratio. Despite being one of the most popular classes of DESs, very little is currently understood regarding how their solvation characteristics differ among themselves as well as from other classes of DESs. Previous studies have shown that the catalytic activity, extraction efficiency, and solubilities of natural compounds can be influenced by the type of choline salt and HBD, but no study has systematically related these parameters with its solvation properties. In this study, inverse gas chromatography was employed for the first time to study the individual solvation characteristics of a diverse range of choline chloride and acetate-based DESs composed of amides, diols, organic acids, and carbohydrates. Furthermore, hydrogen bond acidity was found to be a dominant interaction for all DESs, especially those composed of organic acid HBDs. DESs composed of isomers of butane diol and hexane diol as HBDs exhibited higher hydrogen bond basicity and dispersive-type interactions compared to those composed of urea, acetamide, and organic acids with the position of the hydroxyl functional group influencing their solvation interactions. Choline acetate-based DESs exhibited lower dipolarity and stronger hydrogen bond basicity and dispersive-type interactions compared to choline chloride DESs. The solvation models developed in this study were used to interpret and explain DES behavior in previously reported studies involving catalysis, organic synthesis, and extractions demonstrating that the measured solute–solvent interactions can serve to predict the performance of choline chloride-based DESs when used in various applications.

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

Influence of Van der Waals Interactions on the Solvation Energies of Adsorbates at Pt-Based Electrocatalysts

Solvation can significantly modify the adsorption energy of species at surfaces, thereby influencing the performance of electrocatalysts and liquid-phase catalysts. Thus, it is important to understand adsorbate solvation at the nanoscale. Here we evaluate the effect of van der Waals (vdW) interactions described by different approaches on the solvation energy of *OH adsorbed on near-surface alloys (NSAs) of Pt. Our results show that the studied functionals can be divided into two groups, each with rather similar average *OH solvation energies: (1) PBE and PW91; and (2) vdW functionals, RPBE, PBE-D3 and RPBE-D3. On average, *OH solvation energies are less negative by ~0.14 eV in group (2) compared to (1), and the values for a given alloy can be extrapolated from one functional to another within the same group. Depending on the desired level of accuracy, these concrete observations and our tabulated values can be used to rapidly incorporate solvation into models for electrocatalysis and liquid-phase catalysis.

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

Role of solvation structure in the shuttling of the hydrated excess proton

The classic Marcus electron transfer reaction model demonstrated that a barrierless electron transfer reaction can occur when both the reactant and product have almost similar solvation environment. In our recently developed proton model, we have incorporated the pre-solvation concept and showed that it indeed facilitates the proton diffusion in aqueous solution. Here in this work, we further quantify the degree of pre-solvation using different structural parameters, e.g., tetrahedral order parameter, average numbers of hydrogen bonds. All the above said parameters exhibit a very strong correlation with the proton share parameter. The more Zundel-like configurations have almost identical solvation environment for both the water molecules and support the pre-solvation concept. However, in the case of Eigen-like configurations, the central hydronium and “special pair” water have distinctly different solvation structures.

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

Entropic Penalty Switches Li + Solvation Site Formation and Transport Mechanisms in Mixed Polarity Copolymer Electrolytes

Emerging solid polymer electrolyte (SPE) designs for efficient Li-ion (Li + ) conduction have relied on polarity and mobility contrast to improve conductivity. To further develop this concept, we employ simulations to examine Li + solvation and transport in poly(oligo ethylene methacrylate) (POEM) and its copolymers with poly(glycerol carbonate methacrylate) (PGCMA). We find that Li + is solvated by ether oxygens instead of the highly polar PGCMA, due to lower entropic penalties. The presence of PGCMA promotes single-chain solvation, thereby suppressing interchain Li + hopping. The conductivity difference between random copolymer PGCMA-r-POEM and block copolymer PGCMA-b-POEM is explained in terms of a hybrid solvation site mechanism. With diffuse microscopic interfaces between domains, PGCMA near the POEM contributes to Li + transport by forming hybrid solvation sites. The formation of such sites is hindered when PGCMA is locally concentrated. These findings help explain how thermodynamic driving forces govern Li + solvation and transport in mixed SPEs.

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

Understanding the influence of capillary waves on solvation at the liquid-vapor interface

This work investigates the question if surface capillary waves (CWs) affect interfacial solvation thermodynamic properties that determine the propensity of small molecules toward the liquid-vapor interface. We focus on (1) the evaluation of these properties from molecular simulations in a practical manner and (2) understanding them from the perspective of theories in solvation thermodynamics, especially solvent reorganization effects. Concerning the former objective, we propose a computational method that exploits the relationship between an external field acting on the liquid-vapor interface and the magnitude of CWs. The system considered contains the solvent, an externally applied field (f) and the solute molecule fixed at a particular location. The magnitude of f is selected to induce changes in CWs. The difference between the solvation free energies computed in the presence and in the absence of f is then shown to quantify the contribution of CWs to interfacial solvation. We describe the implementation of this method in the canonical ensemble by using a Lennard-Jones solvent and a non-ionic solute. Results are shown for three types of solutes that differ in the nature of short-ranged repulsive (hard-core) interactions. Overall, we observe that CWs have a negligible or very small effect on the interfacial solvation free energy of a solute molecule fixed near the liquid-vapor interface for the above systems. We also explain how the effects of pinning or dampening of CWs caused by a fixed solute are effectively compensated and do not contribute to the solvation free energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗