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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Confined Lithium–Sulfur Reactions in Narrow-Diameter Carbon Nanotubes Reveal Enhanced Electrochemical Reactivity

Here, we demonstrate an unusual electrochemical reaction of sulfur with lithium upon encapsulation in narrow-diameter (subnanometer) single-walled carbon nanotubes (SWNTs). Our study provides mechanistic insight on the synergistic effects of sulfur confinement and Li + ion solvation properties that culminate in a new mechanism of these sub-nanoscale-enabled reactions (which cannot be solely attributed to the lithiation–delithiation of conventional sulfur). Two types of SWNTs with distinct diameters, produced by electric arc (EA-SWNTs, average diameter 1.55 nm) or high-pressure carbon monoxide (HiPco-SWNTs, average diameter 1.0 nm), are investigated with two comparable electrolyte systems based on tetraethylene glycol dimethyl ether (TEGDME) and 1,4,7,10,13-pentaoxacyclopentadecane (15-crown-5). Electrochemical analyses indicate that a conventional solution-phase Li–S reaction occurs in EA-SWNTs, which can be attributed to the smaller solvated [Li(TEGDME)] + and [Li(15-crown-5)] + ions within the EA-SWNT diameter. In stark contrast, the Li–S confined in narrower diameter HiPco-SWNTs exhibits unusual electrochemical behavior that can be attributed to a solid-state reaction enabled by the smaller HiPco-SWNT diameter compared to the size of solvated Li + ions. Our results of the electrochemical analyses are corroborated and supported with various spectroscopic analyses including operando Raman, X-ray photoelectron spectroscopy, and first-principles calculations from density functional theory. Taken together, our findings demonstrate that the controlled solid-state lithiation–delithiation of sulfur and an enhanced electrochemical reactivity can be achieved by sub-nanoscale encapsulation and one-dimensional confinement in narrow-diameter SWNTs.

25 ENERGY STORAGE↗
DOE OSTI2025

Revealing Solvent-Assisted Li+ Transport in the Solid Electrolyte Interphase operando

The performance of energy-dense lithium metal batteries is critically influenced by the properties of the solid electrolyte interphase (SEI). Yet, progress in understanding this layer has been limited by the lack of accurate operando characterization because the SEI evolves dynamically during cycling. Here, we apply dynamic electrochemical impedance spectroscopy (dEIS) to resolve the real-time evolution of the SEI on lithium metal in ether-based electrolytes with varying degrees of fluorination. We find that faster stabilization of the compact SEI resistance correlates with improved passivation and higher Coulombic efficiency. Unexpectedly, compact SEI resistance correlates directly with Li+ solvation energy, revealing that weaker Li+ solvation increases not only bulk but also interphase resistance. These findings challenge the conventional view of the SEI as a purely solid-phase conductor and instead support a solvent-assisted Li+ transport mechanism within the compact SEI. This framework emphasizes the need to balance SEI ionic conductivity with the Li+ solvation environment to maximize lithium metal battery performance.

Sacci, Robert [ORNL] (ORCID:0000000200735221)↗
DOE OSTI2016

Interplay of Ion–Water and Water–Water Interactions within the Hydration Shells of Nitrate and Carbonate Directly Probed with 2D IR Spectroscopy

The long-range influence of ions in solution on the water hydrogen-bond (H-bond) network remains a topic of vigorous debate. Recent spectroscopic and theoretical studies have, for the most part, reached the consensus that weakly coordinating ions only affect water molecules in the first hydration shell. In this work, we apply ultrafast broadband two-dimensional infrared (2D IR) spectroscopy to aqueous nitrate and carbonate in neat H 2 O to study the solvation structure and dynamics of ions on opposite ends of the Hofmeister series. By exciting both the water OH stretches and ion stretches and probing the associated cross-peaks between them, we are afforded a comprehensive view into the complex nature of ion hydration. We show in aqueous nitrate that weak ion–water H-bonding leads to water–water interactions in the ion solvation shells dominating the dynamics. In contrast, the carbonate CO stretches show significant mixing with the water OH stretches due to strong ion–water H-bonding such that the water and ion modes are intimately correlated. Further, the excitonic nature of vibrations in neat H 2 O, which spans multiple water molecules, is an important factor in describing ion hydration. Finally, we attribute these complex dynamics to the likely presence of intermediate-range effects influenced by waters beyond the first solvation shell.

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

Towards a unified picture of the water self-ions at the air-water interface: a density functional theory perspective

The propensities of the water self ions, H3O+ and OH- , for the air-water interface has implications for interfacial acid-base chemistry. Despite numerous experimental and computational studies, no consensus has been reached on the question of whether or not H3O+ and/or OH- prefer to be at the water surface or in the bulk. Here we report a molecular dynamics simulation study of the bulk vs. interfacial behavior of H3O+ and OH- that employs forces derived from density functional theory with a generalized gradient approximation exchangecorrelation functional (specifically, BLYP) and empirical dispersion corrections. We computed the potential of mean force (PMF) for H3O+ as a function of the position of the ion in a 215-molecule water slab. The PMF is flat, suggesting that H3O+ has equal propensity for the air-water interface and the bulk. We compare the PMF for H3O+ to our previously computed PMF for OH- adsorption, which contains a shallow minimum at the interface, and we explore how differences in solvation of each ion at the interface vs. the bulk are connected with interfacial propensity. We find that the solvation shell of H3O+ is only slightly dependent on its position in the water slab, while OH- partially desolvates as it approaches the interface, and we examine how this difference in solvation behavior is manifested in the electronic structure and chemistry of the two ions. DJT was supported by National Science Foundation grant CHE-0909227. CJM was supported by the U.S. Department of Energy‘s (DOE) Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated for the Department of Energy by Battelle. The potential of mean force required resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC05-00OR22725. The remaining simulations and analysis used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. at at Lawrence Berkeley National Laboratory. MDB is grateful for the support of the Linus Pauling Distinguished Postdoctoral Fellowship Program at PNNL.

Baer, Marcel D.↗
DOE OSTI2014

Solvent friction effects propagate over the entire protein molecule through low-frequency collective modes

Protein solvation dynamics has been investigated using atom-dependent Langevin friction coefficients derived directly from molecular dynamics (MD) simulations. To determine the effect of solvation on the atomic friction coefficients, solution and vacuum MD simulations were performed for lysozyme and staphylococcal nuclease and analyzed by Langevin mode analysis. The coefficients thus derived are roughly correlated with the atomic solvent-accessible surface area (ASA), as expected from the fact that friction occurs as the result of collisions with solvent molecules. However, a considerable number of atoms with higher friction coefficients are found inside the core region. Hence, the influence of solvent friction propagates into the protein core. The internal coefficients have large contributions from the low-frequency modes, yielding a simple picture of the surface-to-core long-range damping via solvation governed by collective low-frequency modes. To make use of these findings in implicit-solvent modeling, we compare the all-atom friction results with those obtained using Langevin dynamics (LD) with two empirical representations: the constant-friction and the ASA-dependent (Pastor Karplus) friction models. The constant-friction model overestimates the core and underestimates the surface damping whereas the ASA-dependent friction model, which damps protein atoms only on the solvent-accessible surface, reproduces well the friction coefficients for both the surface and core regions observed in the explicit-solvent MD simulations. Furthermore, in LD simulation, the solvent friction coefficients should be imposed only on the protein surface.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Understanding the apparent fractional charge of protons in the aqueous electrochemical double layer

A detailed atomic-scale description of the electrochemical interface is essential to the understanding of electrochemical energy transformations. In this work, we investigate the charge of solvated protons at the Pt(111) | H 2 O and Al(111) | H 2 O interfaces. Using semi-local density functional theory as well as hybrid functionals and embedded correlated wavefunction methods as higher-level benchmarks, we show that the effective charge of a solvated proton in the electrochemical double layer or outer Helmholtz plane at all levels of theory is fractional, when the solvated proton and solvent band edges are aligned correctly with the Fermi level of the metal ( E F ). In conclusion, the observed fractional charge in the absence of frontier band mis-alignment arises from a significant overlap between the proton and the electron density from the metal surface, and results in an energetic difference between protons in bulk solution and those in the outer Helmholtz plane.

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

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗
DOE OSTI2013

An ab initio approach to understanding the specific ion effect

Although there is a consensus that large, polarizable anions will adsorb to the air-water interface, the precise interactions that give rise to surface enhancement are still being debated. Previously, we have demonstrated that there is a significant dependence on the choice of molecular interaction potential for iodide adsorption at the air-water interface. Specifically, density functional theory (DFT) based interaction potentials lead to significantly less adsorption than empirical interaction potentials that include polarization. We have also demonstrated that DFT based interaction potentials can accurately capture the structure of the first solvation shell of both simple and polyatomic anions as compared to multi-edge X-ray absorption fine structure (XAFS) experiments. We continue to push DFT by examining the surface propensity and structure of polyatomic anions that exhibit both chaotropic and kosmotropic characteristics as classified by the Hofmeister series. We will investigate the role of local solvation effects of SCN- and compare to previous studies on I- and IO3- under both bulk and interfacial conditions. A picture emerges where we can correlate local solvation stucture to the Hofmeister series.

Ion adsorption to the water/air Interface↗
DOE OSTI2015

A fundamental study on the [(μ-Cl) 3 Mg 2 (THF) 6 ] + dimer electrolytes for rechargeable Mg batteries

We present a fundamental study on [(μ-Cl) 3 Mg 2 (THF) 6 ] + dimer electrolytes using various physical methods including Subambient Pressure Ionization with Nanoelectrospray Mass spectrometry (SPIN-MS), Raman spectroscopy, 25 Mg{ 1 H} NMR, 27 Al{ 1 H} NMR and electrochemical analysis. For the first time, long time sought THF solvated [MgCl] + species was experimentally characterized by SPIN mass spectrometry in the solution of the Mgdimer containing electrolyte, confirming the mono-Cl- abstraction reaction between MgCl 2 and an Al Lewis acid. Solvated MgCl 2 in the electrolyte was confirmed by Raman spectroscopy. The experimental results establish the previously proposed dimerization equilibrium of solvated [MgCl] + and MgCl2 with [(μ-Cl) 3 Mg 2 (THF)6] + . 25 Mg{ 1 H} NMR, 27 Al{ 1 H} NMR and electrochemical analysis on chloration reaction of [(μ-Cl) 3 Mg 2 (THF) 6 ]AlPh 3 Cl with external Cl- led to further insights on the coordination chemistry of the dimer electrolyte. Finally, a comprehensive mechanism is proposed for the reversible electrochemical Mg deposition and stripping and Mg 2+ and Cl- ion transports of the Mg dimer electrolytes in rechargeable Mg batteries.

25 ENERGY STORAGE↗
DOE OSTI2017

A density functional theory based approach for predicting melting points of ionic liquids

Accurate prediction of melting points of ILs is important both from the fundamental point of view and from the practical perspective for screening ILs with low melting points and broadening their utilization in a wider temperature range. In this work, we present an ab initio approach to calculating melting points of ILs with known crystal structures and illustrate its application for a series of 11 ILs containing imidazolium/pyrrolidinium cations and halide/polyatomic fluoro-containing anions. The melting point is determined as a temperature at which the Gibbs free energy of fusion is zero. The Gibbs free energy of fusion can be expressed through the use of the Born-Fajans-Haber cycle via the lattice free energy of forming a solid IL from gaseous phase ions and the sum of the solvation free energies of ions comprising IL. Dispersion-corrected density functional theory (DFT) involving (semi)local (PBE-D3) and hybrid exchange-correlation (HSE06-D3) functionals is applied to estimate the lattice enthalpy, entropy, and free energy. The ions solvation free energies are calculated with the SMD-generic-IL solvation model at the M06-2X/6-31+G(d) level of theory under standard conditions. The melting points of ILs computed with the HSE06-D3 functional are in good agreement with the experimental data, with a mean absolute error of 30.5 K and a mean relative error of 8.5%. The model is capable of accurately reproducing the trends in melting points upon variation of alkyl substituents in organic cations and replacement one anion by another. The results verify that the lattice energies of ILs containing polyatomic fluoro-containing anions can be approximated reasonably well using the volume-based thermodynamic approach. However, there is no correlation of the computed lattice energies with molecular volume for ILs containing halide anions. Moreover, entropies of solid ILs follow two different linear relationships with molecular volume for halides and polyatomic fluoro-containing anions. As a result, continuous progress in predicting crystal structures of organic salts with halide anions will be a key factor for successful prediction of melting points with no prior knowledge of the crystal structure.

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

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of Δ G ele , Δ G vDW , Δ G polar , Δ G non-polar , T Δ S total , Δ G deform , Δ G con , and Δ G ion . To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (Δ G deform ), the energy cost from translational and rotational entropic contributions ( T Δ S tran+rot ), the total electrostatic interactions (Δ G polar-PB/GB + Δ G ele ) of incorporation, the intramolecular electrostatic interactions (Δ G ele ) and electrostatic polar solvation interactions (Δ G polar-PB/GB ) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (Δ G vDW ), the non-polar solvation interaction (Δ G non-polar ), the vibrational entropic contribution ( T Δ S vib ), and the standard concentration dependent free energies of DOX (Δ G con ) and the ionic solution (Δ G ion ). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX–DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. Here, we have calculated the total binding free energy (BFE) (Δ G t-sim ) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX–DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX–DNA1 and DOX–DNA2 in the 1 : 1 complex is –9.1 and –5.1 kcal mol –1 respectively. The same quantities from the modified MM-PBSA method are –12.74 and –8.35 kcal mol –1 respectively. The value of the total BFE Δ G t-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of –7.7 ± 0.3 kcal mol –1 .

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

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

The identity of metal ions surrounding DNA is key to its biological function and materials applications. In this work, we compare atomistic molecular dynamics simulations of double strand DNA (dsDNA) with four alkaline earth metal ions (Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ ) to elucidate the physical interactions that govern DNA–ion binding. Simulations accurately model the ion–phosphate distance of Mg 2+ and reproduce ion counting experiments for Ca 2+ , Sr 2+ , and Ba 2+ . Our analysis shows that alkaline earth metal ions prefer to bind at the phosphate backbone compared to the major groove and negligible binding occurs in the minor groove. Larger alkaline earth metal ions with variable first solvation shells (Ca 2+ , Sr 2+ , and Ba 2+ ) show both direct and indirect binding, where indirect binding increases with ion size. Mg 2+ does not fit this trend because the strength of its first solvation shell predicts indirect binding only. Ions bound to the phosphate backbone form fewer contacts per ion compared to the major groove. Within the major groove, metal ions preferentially bind to guanine–cystosine base pairs and form simultaneous contacts with the N7 and O6 atoms of guanine. Overall, we find that the interplay among ion size, DNA–ion interaction, and the size and flexibility of the first solvation shell are key to predicting how alkaline earth metal ions interact with DNA.

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

Synthesis, modular composition, and electrochemical properties of lamellar iron sulfides

Transition metal chalcogenides with layered structures are promising materials for energy storage, catalysis, and electronics, among other areas. We have identified a new layered phase of iron sulfide containing interlayer solvated cations. We present an optimized synthesis for the Li + -containing material from an Fe(III) xanthate complex. Structure and composition data indicate the material consists of poorly-ordered iron sulfide layers separated by solvated cations. The lamellar spacing in these materials can be tuned by changing the identity of the cation. Furthermore, the lamellar spacing can also be reversibly tuned by the degree of solvation of the material. The material is electrically conductive and can serve as a pseudocapacitor with comparable performance to commercial materials such as MnO 2 . Finally, these materials also show promise as lithium or sodium ion battery cathodes with good capacity and reversibility.

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

Influence of diluent concentration in localized high concentration electrolytes: elucidation of hidden diluent-Li + interactions and Li + transport mechanism

Localized high concentration electrolytes (LHCE) offer a viable dilution strategy for high concentration electrolytes (HCE) as the dilution process barely impacts the enhanced reductive/oxidative behavior of the HCE formulation but significantly lowers the overall viscosity and, in most cases, increases the ionic conductivity. On the other hand, experimental studies indicate that fluorinated ether electrolytes such as 1,1,2,2-tetrafluoroethylene 2,2,3,3-tetrafluoropropyl ether (TTE) help grow enhanced passivation layers on Ni-rich NMC cathodes. In this work, we study LHCE formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and TTE as the diluent. Here we use molecular dynamics methodologies, and Raman spectra measurements, to address to what extent the diluent content impacts the coordination behavior of the aggregated structures in the electrolyte and to evaluate the Li + transport properties under the influence of an external electric field. In contrast to other fluorinated ethers, we find that TTE interacts with Li + via fluorine atoms, partially limiting the DMC–Li + interactions hence altering the Li + solvation coordination. This competitive interaction with Li+ between the organic solvent and the TTE diluent influences the electrolyte's reductive/oxidative behavior. Nevertheless, the bonding strength of the Li + –F TTE is much weaker than those of the Li + –O DMC and Li + –OFSI – . Therefore, the existence of Li + –F TTE is in a transient state rather than in a steady state. These results provide plausible guiding rules for future dilution strategies of HCE electrolytes. We also demonstrate that Li + ions drift under the electric field's influence via repeated ion dissociation/association processes following a hopping conduction mechanism. Li + ions jump between aggregated networks where Li–O interactions dominate via diluent-enriched phases, a process in which the solvation shells temporarily mutate to a Li–F dominated coordination structure. We expect our results to contribute an improved atomic-level understanding of the solvation structure and dynamics of LHCE electrolytes.

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

Quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) of zwitterionic organic HCl salts

In this work, we benchmark quadrupolar NMR crystallography guided crystal structure prediction (QNMRX-CSP) for determining the crystal structures of two zwitterionic organic HCl salts, L-ornithine HCl ( Orn ) and L-histidine HCl·H 2 O ( Hist ). These salts present an interesting challenge for QNMRX-CSP, as gas-phase geometry optimizations used to generate starting structures for the organic zwitterionic fragments fail to capture their correct solid-state geometries. To overcome this limitation, geometry optimizations using the COSMO water-solvation model are employed to generate initial structural models. Using this approach, QNMRX-CSP yields structural models of the two zwitterionic organic HCl salts that closely match experimentally determined crystal structures. In addition, the application of QNMRX-CSP to Hist represents a further step toward the de novo structural determination of solvated organic HCl salts, as Hist is the first benchmark system of this type to include a water molecule as a component of its crystal structure. This work is significant for its potential application to the structural determination of active pharmaceutical ingredients, which often feature complex organic components and solvated solid forms.

Fleischer, Carl H. [Florida State Univ., Tallahass↗
DOE OSTI2013

Hydration structure of salt solutions from ab initio molecular dynamics

The solvation structures of Na{sup +}, K{sup +}, and Cl{sup -} ions in aqueous solution have been investigated using density functional theory (DFT) based Car-Parrinello (CP) molecular dynamics (MD) simulations. CPMD trajectories were collected for systems containing three NaCl or KCl ion pairs solvated by 122 water molecules using three different but commonly employed density functionals (BLYP, HCTH, and PBE) with electron correlation treated at the level of the generalized gradient approximation (GGA). The effect of including dispersion forces was analyzed through the use of an empirical correction to the DFT-GGA scheme. Special attention was paid to the hydration characteristics, especially the structural properties of the first solvation shell of the ions, which was investigated through ion-water radial distribution functions, coordination numbers, and angular distribution functions. There are significant differences between the present results obtained from CPMD simulations and those provided by classical MD based on either the CHARMM force field or a polarizable model. Overall, the computed structural properties are in fair agreement with the available experimental results. In particular, the observed coordination numbers 5.0-5.5, 6.0-6.4, and 6.0-6.5 for Na{sup +}, K{sup +}, and Cl{sup -}, respectively, are consistent with X-ray and neutron scattering studies but differ somewhat from some of the many other recent computational studies of these important systems. Possible reasons for the differences are discussed.

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

Molecular density functional theory for water with liquid-gas coexistence and correct pressure

The solvation of hydrophobic solutes in water is special because liquid and gas are almost at coexistence. In the common hypernetted chain approximation to integral equations, or equivalently in the homogenous reference fluid of molecular density functional theory, coexistence is not taken into account. Hydration structures and energies of nanometer-scale hydrophobic solutes are thus incorrect. In this article, we propose a bridge functional that corrects this thermodynamic inconsistency by introducing a metastable gas phase for the homogeneous solvent. We show how this can be done by a third order expansion of the functional around the bulk liquid density that imposes the right pressure and the correct second order derivatives. Although this theory is not limited to water, we apply it to study hydrophobic solvation in water at room temperature and pressure and compare the results to all-atom simulations. The solvation free energy of small molecular solutes like n-alkanes and hard sphere solutes whose radii range from angstroms to nanometers is now in quantitative agreement with reference all atom simulations. The macroscopic liquid-gas surface tension predicted by the theory is comparable to experiments. This theory gives an alternative to the empirical hard sphere bridge correction used so far by several authors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

C–IBI: Targeting cumulative coordination within an iterative protocol to derive coarse-grained models of (multi-component) complex fluids

We present a coarse-graining strategy that we test for aqueous mixtures. The method uses pair-wise cumulative coordination as a target function within an iterative Boltzmann inversion (IBI) like protocol. We name this method coordination iterative Boltzmann inversion (C–IBI). While the underlying coarse-grained model is still structure based and, thus, preserves pair-wise solution structure, our method also reproduces solvation thermodynamics of binary and/or ternary mixtures. In addition, we observe much faster convergence within C–IBI compared to IBI. To validate the robustness, we apply C–IBI to study test cases of solvation thermodynamics of aqueous urea and a triglycine solvation in aqueous urea.

74 ATOMIC AND MOLECULAR PHYSICS↗