Solvation Layer Research Articles

Orientational Ordering in Nano-confined Polar Liquids.

Water and other polar liquids exhibit nanoscale structuring near charged interfaces. When a polar liquid is confined between two charged surfaces, the interfacial solvent layers begin to overlap, resulting in solvation forces. Here, we perform molecular dynamics simulations of polar liquids with different dielectric constants and molecular shapes and sizes confined between charged surfaces, demonstrating strong orientational ordering in the nanoconfined liquids. To rationalize the observed structures, we apply a coarse-grained continuum theory that captures the orientational ordering and solvation forces of those liquids. Our findings reveal the subtle behavior of different nanoconfined polar liquids and establish a simple law for the decay length of the interfacial orientations of the solvents, which depends on their molecular size and polarity. These insights shed light on the nature of solvation forces, which are important in colloid and membrane science, scanning probe microscopy, and nano-electrochemistry.

Examination of Solvent Interactions with Trp-Cage in 1,1,1,3,3,3-Hexafluoro-2-propanol-water at 298 K through MD Simulations and Intermolecular Nuclear Overhauser Effects.

MD simulations of the peptide Trp-cage dissolved in 28% hexafluoroisopropanol (HFIP)-water have been carried out at 298 K with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The work was motivated by the observation that most experimental fluoroalcohol-peptide cross-relaxation terms at 298 K are small, both positive and negative, and not always well predicted from simulations. The cross-relaxation terms for hydrogens of the caged tryptophan residue of Trp-cage are substantially negative, a result consistent with simulations. It was concluded that hexafluoroisopropanol interactions near this part of the peptide are particularly long-lived. While both HFIP and water are present in all regions of the simulation box, the composition of the solvent mixture is not homogeneous throughout the system. HFIP generally accumulates near the peptide surface, while water molecules are preferentially found in regions that are more than 1.5 nm from the surface of the peptide. However, some water remains in higher-than-expected amounts in the solvent layer surrounding 6Trp, 9Asp, Ser13, and Ser14 residues in the helical region of Trp-cage. As observed in simulations of this system at 278 K, HFIP molecules aggregate into clusters that continually form and re-form. Translational diffusion of both HFIP and water appears to be slowed near the surface of the peptide with reduction in diffusion near the 6Trp residue 2- to 3-fold larger than calculated for solvent interactions with other regions of Trp-cage.

Are all poorly soluble drugs dissolved in deep eutectic solvents true solutions?

Hypothesis: The ability of deep eutectic solvents (DES) to enhance solubility of poorly soluble drugs has attracted increasing attention. Researchers have shown that drugs could be dissolved well in DES. In this study, we propose a new existence state of drugs in DES: a quasi-two-phase colloidal system.Experiments: Six poorly soluble drugs were used as the models. The formation of colloidal systems was observed visually by the Tyndall effect and DLS. TEM and SAXS were performed to obtain their structure information. The intermolecular interactions between components were probed via DSC and 1H1H-ROESY. In addition, the properties of colloidal systems were further studied.Findings: Our key finding is that several drugs like lurasidone hydrochloride (LH) could form stable colloids in [Th (thymol)] - [Da (decanoic acid)] DES, resulting from weak interactions between drugs and DES, which is different from the true solution of drugs like ibuprofen where strong interactions were formed. In this LH-DES colloidal system, DES solvation layer was directly observed on the surface of drug particles. In addition, the colloidal system with polydispersity shows superior physical and chemical stability. Different to the prevailing view that substances are fully dissolved in DES, this study discovers another existence state as stable colloidal particles in DES.

Optimizing NaF‐Rich Solid Electrolyte Interphase for Stabilizing Sodium Metal Batteries by Electrolyte Additive

AbstractNaF‐rich electrode–electrolyte interphases play crucial roles in determining the cycling stability of sodium metal batteries (SMBs) because of their electronic insulation and mechanical stability. In this work, perfluorobenzene (PFB) is proposed as the additive to contribute the formation of NaF‐rich solid electrolyte interphases (SEI). PFB at the periphery of the solvation layer can pull out a part of the EC with the lowest solvation energy by Van der Waals forces, thus allowing more to participate in the Na+ solvation layer and form an anion‐aggregated solvation sheath, thus promoting the decomposition of to produce NaF. In addition, PFB has a higher highest occupied molecular orbital and lower lowest unoccupied molecular orbital energy level, which also preferentially decomposes to produce NaF at both electrodes. Benefiting from the intensified NaF ratio in SEI, the Na||Na symmetric cells with such an electrolyte achieves a superior cycling life over 350 h at 1 mA cm−2, and the Na||Na3V2(PO4)2O2F batteries also realize ultrahigh cycling performance with 88.8% capacity retention after 500 cycles.

Liquid Water: When Hyperpolarizability Fluctuations Boost and Reshape the Second Harmonic Scattering Intensities

Second harmonic scattering (SHS) is a method of choice to investigate the molecular structure of liquids. While a clear interpretation of SHS intensity exists for diluted solutions of dyes, the scattering due to solvents remains difficult to interpret quantitatively. Here, we report a quantum mechanics/molecular mechanics (QM/MM) approach to model the polarization-resolved SHS intensity of liquid water, quantifying different contributions to the signal. We point out that the molecular hyperpolarizability fluctuations and correlations cannot be neglected. The intermolecular orientational and hyperpolarizability correlations up to the third solvation layer strongly increase the scattering intensities and modulate the polarization-resolved oscillation that is predicted here by QM/MM without fitting parameters. Our approach can be generalized to other pure liquids to provide a quantitative interpretation of SHS intensities in terms of short-range molecular ordering.

Multiscale structure and precipitation mechanism of debranched starch precipitated by different alcohols

Alcohol solution is a cheap, simple, and effective precipitating solvent frequently used for separating debranched starch (DBS), yet little is known about the precipitation mechanism of DBS by different alcohols. This study precipitated DBS from pullulanase-hydrolyzed starch using ethanol, n-butanol, and isopentanol. The multiscale structures of DBS were characterized, including chain length, single/double helix, and crystalline. The chain conformation and precipitation mechanism of DBS in different alcohols was investigated using molecular dynamics (MD) simulation. DBS precipitated by n-butanol contained the largest proportion of short chain (DP6–24, 83.2 %), the highest V-type crystallinity (21.1 %), and the largest single-helix content (24.7 %). A single helix conformation of DBS chain was determined in alcohols, where alcohol molecules entered the helix cavity. Intra/inter-molecular hydrogen bonds stabilized the helix, with a large number of hydrogen bonds leading to strong molecular interaction and stable helical structure. The solvent accessible surface area of DBS chain decreased by 7.88–19.32 % in alcohols, and the radial distribution function revealed that the first solvent layer of DBS chain at 0.29 nm was closely related to hydrogen bonding. This study provides a basis for the choice of precipitation solvent for preparing DBS with different chain lengths and physicochemical properties.

Critical role of water in fused silica glass by silica nanopowder sintering: Nanoparticle dispersion

Thermoplastic nanocomposite preparation using a silica nanopowder/polyethylene glycol (PEG) / polyvinyl butyral (PVB) system followed by heat-treatment sintering to obtain fused silica glass has high precision, scalable manufacturing processes, and short processing time. However, conventional solvents are toxic and harmful to the environment, and for this reason, the ethanol/water solvent system was proposed in this study. In this study, the effect of water content on the dispersion effect of nano-silica suspensions was investigated through particle size and zeta potential tests and FTIR spectra analysis; in addition, the study proposes a dispersion theory based on the solvation layer of the ethanol/water competitive adsorption and viscosity effect. Finally, optical photographs and SEM characterization showed that the obtained fused silica glass was transparent and had no seam defects, and the XRD and hydroxyl group tests showed that the glass had a low hydroxyl group content of only 2.55 ppm.

Surface-bubble-modulated liquid chromatography: an experimental strategy for identification of molecular processes of solute retention in reversed-phase separation systems.

Molecular level understanding of the chemistry at the aqueous/hydrophobe interface is crucial to separation processes in aqueous media, such as reversed-phase liquid chromatography (RPLC) and solid-phase extraction (SPE). Despite significant advances in our knowledge of the solute retention mechanism in these reversed-phase systems, direct observation of the behavior of molecules and ions at the interface in reversed-phase systems still remains a major challenge and experimental probing techniques that provide the spatial information of the distribution of molecules and ions are required. This review addresses surface-bubble-modulated liquid chromatography (SBMLC), which has a stationary gas phase in a column packed with hydrophobic porous materials and enables one to observe the molecular distribution in the heterogeneous reversed-phase systems consisting of the bulk liquid phase, the interfacial liquid layer, and the hydrophobic materials. The distribution coefficients of organic compounds referring to their accumulations onto the interface of alkyl- and phenyl-hexyl-bonded silica particles exposed to water or acetonitrile-water and into the bonded layers from the bulk liquid phase are determined by SBMLC. The experimental data obtained by SBMLC show that the water/hydrophobe interface exhibits an accumulation selectivity for organic compounds, which is quite different from that of the interior of the bonded chain layer, and the overall separation selectivity of the reversed-phase systems is determined by the relative sizes of the aqueous/hydrophobe interface and the hydrophobe. The solvent composition and the thickness of the interfacial liquid layer formed on octadecyl-bonded (C18) silica surfaces are also estimated from the bulk liquid phase volume determined by the ion partition method employing small inorganic ions as probes. It is clarified that various hydrophilic organic compounds as well as inorganic ions recognize the interfacial liquid layer formed on the C18-bonded silica surfaces as being different from the bulk liquid phase. The behavior of some solute compounds exhibiting substantially weak retention in RPLC or the so-called negative adsorption, such as urea, sugars, and inorganic ions, can rationally be interpreted with a partition between the bulk liquid phase and the interfacial liquid layer. The spatial distribution of solute molecules and the structural properties of the solvent layer on the C18-bonded layer determined by the liquid chromatographic methods are discussed in comparison to the results obtained by other research groups using molecular simulation methods.

Molecular Dynamic Simulation of Ni-Al Alloy-H2O Reactions Using the ReaxFF Reactive Force Field.

Hydrogen as clean energy can effectively solve the problems of fossil energy shortage and environmental pollution. However, traditional methods of H2 production are generally lacking in application value. The procedure for manufacturing H2 by a reaction between active metals and H2O has received wide attention due to its high efficiency. Profound insights into the mechanism and influencing factors of H2 production from active metals are insufficient. The ReaxFF reaction force field module of the Amsterdam Modeling Suite (AMS) is applied in this paper to simulate the reaction of Ni-Al alloys with H2O. It reveals the reaction route of H2 production at the atomic level. The calculation results show that Al is the most critical active site. Moreover, the H2 production capacity of the alloy varies with the crystal structure and atomic ratio. The H2 production rate decreases due to the influence of the water solvation layer and surface coverage. Oxygen reduces the H2 production capacity because oxygen reduces the active sites for H2O adsorption by forming a stable oxide layer with Al.

Zwitterionic Separator Featured with Superdesolvating Properties for High Performance Lithium-Sulfur Batteries.

Lithium-sulfur chemistry has greatly expanded the boundaries of lithium batteries, but the persistent parasitic reaction of soluble sulfur intermediates with lithium anode remains a primary challenge. Understanding and regulating the solvation structures of lithium ions (Li+) and polysulfides (LiPSs) are critical to addressing the above issues. Herein, inspired by the natural superhydrophilic resistance to contamination, we developed a zwitterionic nanoparticles (ZWP) separator capable of modulating the solvated of Li+ and LiPSs. The dense solvated layer induced by ZWP effectively prevents the movement of LiPSs without compromising Li+ transport. Moreover, the high electrolyte affinity of the ZWP effectively results in minimizing the deposition of LiPSs on the separator. Furthermore, the structure of the solvated Li+ and LiPSs is also unveiled by molecular simulation and nuclear magnetic resonance (NMR). In addition, in situ UV setup proved the ZWP separator can effectively suppress the shuttle of LiPSs. The restricted space formed by the tightly packed ZWP stabilizes the lithium deposition and regulates dendrite growth. Consequently, the performance of lithium-sulfur batteries is significantly improved and good cycle stability is maintained even at high sulfur loadings (5 mg cm-2). This contribution provides a new insight into the rational design of lithium-sulfur battery separators.

Quantification of a new recreational drug 2-methoxyqualone in human hair using gas chromatography–tandem mass spectrometry

A method for identification and quantification of 2-methoxyqualone, an newly emerging quinazolinone derivative recreational drug, in human scalp hair was established using gas chromatography (GC)–tandem mass spectrometry (MS/MS). In this report, authentic cases are presented, in which the suspects were seized by police security bureau; the police in China requested our laboratory to identify and quantify the involved drug(s) of abuse in the hair samples of the suspects. After washing and cryo-grinding the authentic hair samples, the target compound was extracted with methanol, and the solvent layer was evaporated to dryness. The residue was reconstituted in methanol and analyzed by GC–MS/MS. 2-Methoxyqualone concentrations in the hair were between 35.1 and 116 pg/mg. The calibration curve of the substance in hair samples showed a good linearity in the concentration range of 10–1000 pg/mg (r > 0.998); the extraction recovery rate, 88.8–105.6 %; the interday and intraday precisions and accuracies (biases), not greater than 8.9 %. 2-Methoxyqualone in human hair had good stability under three different storage conditions at room (20 °C), refrigerated (4 °C) and frozen (− 20 °C) temperatures for at least 7 days. In the present report, simple and rapid quantification method for 2-methoxyqualone in human scalp hair have been established using GC–MS/MS and it could successfully be applied to authentic forensic toxicological cases. To our knowledge, this is the first report for quantification of 2-methoxyqualone in human hair samples.

Electrolyte Solvation Structure Regulation Promotes Aluminum–Air Batteries to Approach Theoretical Discharge Capacity

Aluminum–air batteries (AABs) feature high energy density and no limitations for charging facilities but suffer from severe self-corrosion reactions. Suppressing the side reaction between water and aluminum is the fundamental strategy to boost AAB performance. Herein, we first disclose that addition of sodium formate (HCOONa) can largely remold the solvation structure of the classic NaOH electrolyte and greatly reduce the amount of free water, thus constructing AABs close to their theoretical capacity. Ab initio molecular dynamics simulations show that multiple Na+ and H2O form hydrated Na+ ion clusters [Nam(OH2)7m]m+ in the HCOONa–NaOH electrolyte. Meanwhile, HCOO– enters the first solvation layer of Na+ and forms a “stable triangle” with Na+ and H2O. The coordination network formed by Na+, H2O, and HCOO– results in a very low level of free water. Moreover, HCOO– further forms a self-assembled passivation layer on the Al surface. As such, the self-corrosion rate of Al in the HCOONa-based electrolyte is drastically reduced. More encouragingly, the ultrahigh capacity (2767 mAh g–1) delivered by AABs using the hybrid electrolyte reaches 92.9% of its theoretical value, concurrently achieving a high energy density of 3702 Wh kg–1. This work paves the way for constructing durable AABs through an electrolyte regulation strategy.

In situ self-assembled macroporous interconnected nanosheet arrays of Ni-1,3,5-benzenetricarboxylate metal − organic framework on Ti mesh as high-performance oxygen evolution electrodes

Highly efficient metal–organic framework (MOF)-based oxygen evolution reaction (OER) catalysts are desirable for water splitting, but their development remains challenging due to poor accessibility of coordinatively unsaturated metal (cum) sites and low intrinsic activity. A large-area three-dimensional (3-D) macroporous interconnected nanosheet array of Ni-1,3,5-benzenetricarboxylate has been in situ self-assembled on Ti mesh (TM) by using ethanol as the solvent and high-affinity oxide layer on TM to promote in situ nucleation. The obtained nanoarchitecture exhibits much superior catalytic activity compared to most reported catalysts including MOF-based catalysts, other precious-metal-free ones, and Ir/Ru-based ones. Additionally, this electrode undergoes no current decay after 300 cyclic voltammetry (CV) cycles and can maintain at 250 mA cm−2 for over 266 h. The excellent catalytic performance is mainly due to the 3-D macroporous and interconnected nanosheet array structure improving cum site exposure and charge transport and in situ activated cum cations enhancing OH− adsorption. This work not only develops a facile and economical approach to synthesize 3-D macroporous interconnected MOF nanosheet arrays to simultaneously increase the number, exposure, and intrinsic activity of active sites and facilitate charge transport for high-performance eletrocatalysis, but provides scientific insights into the mechanisms for self-assembly of this unique nanoarchitecture and for the high OER performance.

Deciphering the Properties of Nanoconfined Aqueous Solutions by Vibrational Sum Frequency Generation Spectroscopy.

When confined between walls at nanometer distances, water exhibits surprisingly different properties with reference to bare interfacial water. Based on computer simulations, we demonstrate how vibrational sum frequency generation (VSFG) spectroscopy can be used-even with very mild symmetry breaking-to discriminate multilayer water in wide slit pores from both bilayer and monolayer water confined within molecularly narrow pores. Applying the technique, the VSFG lineshapes of monolayer, bilayer, and multilayer water are found to differ in characteristic ways, which is explained by their distinct density stratifications giving rise to different H-bonding patterns in the respective solvation layers.

Ionic Liquid-Mediated Route to Atomic Layer Deposition of Tin(II) Oxide via a C-C Bond Cleavage Ligand Modification Mechanism.

Atomic layer deposition (ALD) is a technologically important method to grow thin films with high conformality and excellent thickness control from vapor phase precursors. The development of new thermal ALD processes can be limited by precursor reactivity and stability: reaction temperature and precursor design are among the few variables available to achieve higher reactivity in gas-phase reactions, unlike in solution synthesis, where the use of solvent and/or a catalyst can promote a desired reaction. To bridge this synthesis gap between vapor-phase and solution-phase, we demonstrate the use of an ultrathin coating layer of a vapor phase-compatible solvent─an ionic liquid (IL)─on our growth substrate to perform ALD of SnO. Successful SnO deposition is achieved using tin acetylacetonate and water, a process that otherwise would require a stronger counter-reactant such as ozone. The presence of the layer of IL allows a solvent-mediated reaction mechanism to take place on the growth substrate surface. We report a growth per cycle of 0.67 Å/cycle at a deposition temperature of 100 °C in an IL comprising 1-ethyl-3-methylimidazolium hydrogen sulfate. Characterization of the ALD films confirms the SnO film composition, and 1H and 13C NMR are used to probe the solvent-mediated ALD reaction, suggesting a solvent-mediated addition-elimination-type mechanism which breaks a C-C bond in acetylacetonate to form acetone and acetate. Density functional theory calculations show that the IL solvent is beneficial to the proposed solvent-mediated mechanism by lowering the C-C bond cleavage energetics of acetylacetonate compared to the vapor phase. A general class of ligand modification reactions for thermal ALD is thus introduced in this work.

Extending the π-system of inverted solvatochromic phenolate-based dyes shifts their inversion point due to an increment in their sensitivity to solvent polarizability

Two pairs of phenolate-based dyes with a π-connecting system of different lengths were synthesized, and their spectral variations were investigated in pure solvents of variable polarity and methanol-dichloromethane solvent mixtures. All the dyes studied exhibited inverted solvatochromism, and multiparametric analysis of their solvatochromic absorption bands of these dyes revealed that the solvent acidity and polarizability are mainly responsible for this behavior. The extension of the π-connecting system brings out of the shifts of the spectral behavior exhibited by the dyes towards shorter wavelengths with a concomitant displacement of the position where their solvatochromic inversion takes place towards higher solvent polarity due to an increment in the ability of the extended dyes to sense the solvent polarizability, without modifying the specific solute–solvent interactions in the first or the second solvation layer, a result confirmed by theoretical calculations.

Numerical study of microscopic particle arrangement of suspension flow in a narrow channel for the estimation of macroscopic rheological properties

In a narrow channel, the apparent relative viscosity of a suspension with finite-size particles is strongly dependent on its microscopic particle arrangement. Relative viscosity increases when suspended particles flow near the channel wall; thus, a suspension in a narrow channel does not always exhibit the same rheological properties even if the concentration is the same. In this study, we focus on the inertia and concentration of particles in a narrow channel and consider their effects on the microscopic particle arrangement and macroscopic suspension rheology. Two-dimensional pressure-driven suspension flow simulations were performed using a two-way coupling scheme, and normalized particle density distribution (PDD) were implemented to consider their particle arrangements. The results demonstrated that the velocity profiles for the particle suspension were changed by the Reynolds number and particle concentration because of the interactions between particles according to the power-law index. These changes affected the particle equilibrium positions in the channel, and the subsequent changes in solvent layer thickness caused changes in the macroscopic apparent viscosity. The behavior of microscopic particles played important roles in determining macroscopic rheology. Thus, we have confirmed that a normalized PDD can be used to estimate and assess the macroscopic rheology of a suspension.

Confinement Effects on Distribution and Transport of Neutral Solutes in a Small Hydrophobic Nanopore.

Molecular dynamics simulations are used to study confinement effects in small cylindrical silica pores with extended hydrophobic surface functionalization as realized, for example, in reversed-phase liquid chromatography (RPLC) columns. In particular, we use a 6 nm cylindrical and a 10 nm slit pore bearing the same C18 stationary phase to compare the conditions inside the smaller-than-average pores within an RPLC column to column-averaged properties. Two small, neutral, apolar to moderately polar solutes are used to assess the consequences of spatial confinement for typical RPLC analytes with water (W)-acetonitrile (ACN) mobile phases at W/ACN ratios between 70/30 and 10/90 (v/v). The simulated data show that true bulk liquid behavior, as observed over an extended center region in the 10 nm slit pore, is not recovered within the 6 nm cylindrical pore. Instead, the ACN-enriched solvent layer around the C18 chain ends (the ACN ditch), a general feature of hydrophobic interfaces equilibrated with aqueous-organic liquids, extends over the entire pore lumen of the small cylindrical pore. This renders the entire pore a highly hydrophobic environment, where, contrary to column-averaged behavior, neither the local nor the pore-averaged sorption and diffusion of analytes scales directly with the W/ACN ratio of the mobile phase. Additionally, the solute polarity-related discrimination between analytes is enhanced. The consequences of local ACN ditch overlap in RPLC columns are reminiscent of ion transport in porous media with charged surfaces, where electrical double-layer overlap occurring locally in smaller pores leads to discrimination between co- and counterionic species.

Load‐Dependent Friction Hysteresis for Graphitic Surfaces in n‐Hexadecane

AbstractSliding‐induced friction behavior of a single‐asperity silica probe against few‐layer (FL) graphene and bulk graphite is measured in the presence of n‐hexadecane using an atomic force microscope (AFM). The load‐dependent nanoscale friction measurements display friction hysteresis, i.e., higher friction forces during unloading of the contact than loading, at a given normal load. However, unlike hysteresis in friction of graphene measured in ambient, several unique trends are observed when the contact is immersed in n‐hexadecane. First, the friction hysteresis is measured up to a transition load and beyond which is found negligible; second, a similar behavior is observed on bulk graphite; and third, a friction strengthening of the contact persisted up to several nanometers. Quasi‐static force‐separation curves identify up to four layers of n‐hexadecane solvation layers on graphitic surfaces. Molecular dynamic simulations illustrate that the solvated n‐hexadecane molecules within the contact carry the probe load and determine the generated contact area, affecting friction hysteresis. Further, during AFM probe sliding, instead of a pucker, a molecular pile‐up of n‐hexadecane, in front of the tip is observed. These findings provide new perspectives on understanding of the dissipation mechanisms of graphene that predominantly are surrounded by structured liquid molecules.

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model.

Molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. In this work, we develop and assess two multiscale solvation models for computing solvation thermodynamic properties. The new SLIC/CDC model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with new statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). Given the structures of 500 solutes, the SLIC/CDC model predicts Gibbs energies of solvation in water with an average accuracy better than 1 kcal/mol, when compared to experimental measurements, and better than 0.8 kcal/mol, when compared to explicit-solvent molecular dynamics simulations. The individual SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The previously published SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent-accessible surface area (SASA) nonpolar energy mode. With our new, improved parametrization method, the SLIC/SASA model now predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems, compared to experimental and explicit-solvent molecular dynamics, and better than 1.6 kcal/mol average accuracy in ionic liquids, compared to explicit-solvent molecular dynamics. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.