Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory.

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here, we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F- over Cl- in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory-experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

Thermodynamics of ion binding and occupancy in potassium channels.

Potassium channels modulate various cellular functions through efficient and selective conduction of K+ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K+ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.

Thermodynamics of Hydration from the Perspective of the Molecular Quasichemical Theory of Solutions.

The quasichemical organization of the potential distribution theorem, molecular quasichemical theory (QCT), enables practical calculations and also provides a conceptual framework for molecular hydration phenomena. QCT can be viewed from multiple perspectives: (a) as a way to regularize an ill-conditioned statistical thermodynamic problem; (b) as an introduction of and emphasis on the neighborship characteristics of a solute of interest; or (c) as a way to include accurate electronic structure descriptions of near-neighbor interactions in defensible statistical thermodynamics by clearly defining neighborship clusters. The theory has been applied to solutes of a wide range of chemical complexity, ranging from ions that interact with water with both long-ranged and chemically intricate short-ranged interactions, to solutes that interact with water solely through traditional van der Waals interations, and including water itself. The solutes range in variety from monatomic ions to chemically heterogeneous macromolecules. A notable feature of QCT is that, in applying the theory to this range of solutes, the theory itself provides guidance on the necessary approximations and simplifications that can facilitate the calculations. In this Perspective, we develop these ideas and document them with examples that reveal the insights that can be extracted using the QCT formulation.

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes.

With a longer-term goal of addressing the comparative behavior of the aqueous halides F-, Cl-, Br-, and I- on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl-, Br-, and I-, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

Hydrophilic Interactions Dominate the Inverse Temperature Dependence of Polypeptide Hydration Free Energies Attributed to Hydrophobicity.

We address the association of the hydrophobic driving forces in protein folding with the inverse temperature dependence of protein hydration, wherein stabilizing hydration effects strengthen with increasing temperature in a physiological range. All-atom calculations of the free energy of hydration of aqueous deca-alanine conformers, holistically including backbone and side-chain interactions together, show that attractive peptide-solvent interactions and the thermal expansion of the solvent dominate the inverse temperature signatures that have been interpreted traditionally as the hydrophobic stabilization of proteins in aqueous solution. Equivalent calculations on a methane solute are also presented as a benchmark for comparison. The present study calls for a reassessment of the forces that stabilize folded protein conformations in aqueous solutions and of the additivity of hydrophobic/hydrophilic contributions.

Molecular Dynamics of Lithium Ion Transport in a Model Solid Electrolyte Interphase.

Li+ transport within a solid electrolyte interphase (SEI) in lithium ion batteries has challenged molecular dynamics (MD) studies due to limited compositional control of that layer. In recent years, experiments and ab initio simulations have identified dilithium ethylene dicarbonate (Li2EDC) as the dominant component of SEI layers. Here, we adopt a parameterized, non-polarizable MD force field for Li2EDC to study transport characteristics of Li+ in this model SEI layer at moderate temperatures over long times. The observed correlations are consistent with recent MD results using a polarizable force field, suggesting that this non-polarizable model is effective for our purposes of investigating Li+ dynamics. Mean-squared displacements distinguish three distinct Li+ transport regimes in EDC - ballistic, trapping, and diffusive. Compared to liquid ethylene carbonate (EC), the nanosecond trapping times in EDC are significantly longer and naturally decrease at higher temperatures. New materials developed for fast-charging Li-ion batteries should have a smaller trapping region. The analyses implemented in this paper can be used for testing transport of Li+ ion in novel battery materials. Non-Gaussian features of van Hove self -correlation functions for Li+ in EDC, along with the mean-squared displacements, are consistent in describing EDC as a glassy material compared with liquid EC. Vibrational modes of Li+ ion, identified by MD, characterize the trapping and are further validated by electronic structure calculations. Some of this work appeared in an extended abstract and has been reproduced with permission from ECS Transactions, 77, 1155-1162 (2017).

Role of Solute Attractive Forces in the Atomic-Scale Theory of Hydrophobic Effects.

The role that van der Waals (vdW) attractive forces play in the hydration and association of atomic hydrophobic solutes such as argon (Ar) in water is reanalyzed using the local molecular field (LMF) theory of those interactions. In this problem, solute vdW attractive forces can reduce or mask hydrophobic interactions as measured by contact peak heights of the ArAr correlation function compared to reference results for purely repulsive core solutes. Nevertheless, both systems exhibit a characteristic hydrophobic inverse temperature behavior in which hydrophobic association becomes stronger with increasing temperature through a moderate temperature range. The new theoretical approximation obtained here is remarkably simple and faithful to the statistical mechanical LMF assessment of the necessary force balance. Our results extend and significantly revise approximations made in a recent application of the LMF approach to this problem and, unexpectedly, support a theory of nearly 40 years ago.

Molecular Simulation Results on Charged Carbon Nanotube Forest-Based Supercapacitors.

Electrochemical double-layer capacitances of charged carbon nanotube (CNT) forests with tetraethyl ammonium tetrafluoro borate electrolyte in propylene carbonate are studied on the basis of molecular dynamics simulation. Direct molecular simulation of the filling of pore spaces of the forest is feasible even with realistic, small CNT spacings. The numerical solution of the Poisson equation based on the extracted average charge densities then yields a regular experimental dependence on the width of the pore spaces, in contrast to the anomalous pattern observed in experiments on other carbon materials and also in simulations on planar slot-like pores. The capacitances obtained have realistic magnitudes but are insensitive to electric potential differences between the electrodes in this model. This agrees with previous calculations on CNT forest supercapacitors, but not with experiments which have suggested electrochemical doping for these systems. Those phenomena remain for further theory/modeling work.

Assessment of Simple Models for Molecular Simulation of Ethylene Carbonate and Propylene Carbonate as Solvents for Electrolyte Solutions.

Progress in understanding liquid ethylene carbonate (EC) and propylene carbonate (PC) on the basis of molecular simulation, emphasizing simple models of interatomic forces, is reviewed. Results on the bulk liquids are examined from the perspective of anticipated applications to materials for electrical energy storage devices. Preliminary results on electrochemical double-layer capacitors based on carbon nanotube forests and on model solid-electrolyte interphase (SEI) layers of lithium ion batteries are considered as examples. The basic results discussed suggest that an empirically parameterized, non-polarizable force field can reproduce experimental structural, thermodynamic, and dielectric properties of EC and PC liquids with acceptable accuracy. More sophisticated force fields might include molecular polarizability and Buckingham-model description of inter-atomic overlap repulsions as extensions to Lennard-Jones models of van der Waals interactions. Simple approaches should be similarly successful also for applications to organic molecular ions in EC/PC solutions, but the important case of Li[Formula: see text] deserves special attention because of the particularly strong interactions of that small ion with neighboring solvent molecules. To treat the Li[Formula: see text] ions in liquid EC/PC solutions, we identify interaction models defined by empirically scaled partial charges for ion-solvent interactions. The empirical adjustments use more basic inputs, electronic structure calculations and ab initio molecular dynamics simulations, and also experimental results on Li[Formula: see text] thermodynamics and transport in EC/PC solutions. Application of such models to the mechanism of Li[Formula: see text] transport in glassy SEI models emphasizes the advantage of long time-scale molecular dynamics studies of these non-equilibrium materials.

Anomalous Potential-Dependent Friction on Au(111) Measured by AFM.

We present an exploratory study of the tribological properties between an AFM probe and a Au(111) surface in an aqueous environment while subjected to applied surface potentials. Using a three-electrode setup, the electrical potential and interfacial electric field on a Au(111) working electrode are controlled. Lateral force microscopy is used to measure the friction forces between the AFM probe and the Au surface. As the AFM probe approaches the surface, normal forces are also measured to gain insight into the interfacial forces. When a positive potential is applied to the Au surface, the friction is found to rise sharply at a critical potential and level off at a relatively high value. However, when a negative potential is applied, the friction forces are low, even lower compared to the open circuit potential case. These changes in friction, by a factor of approximately 35, as a function of the applied potential are found to be reversible over multiple cycles. We attribute the origin of the high friction at positive potentials to the formation of a highly confined, ordered icelike water layer at the Au/electrolyte interface that results in effective hydrogen bonding with the AFM probe. At negative potentials, the icelike water layer is disrupted, resulting in the water molecules acting as boundary lubricants and providing low friction. Such friction experiments can provide valuable insight into the structure and properties of water at charged surfaces under various conditions and can potentially impact a variety of technologies relying on molecular-level friction such as MEMs.

Quasi-chemical theory of F-(aq): The "no split occupancies rule" revisited.

We use ab initio molecular dynamics (AIMD) calculations and quasi-chemical theory (QCT) to study the inner-shell structure of F-(aq) and to evaluate that single-ion free energy under standard conditions. Following the "no split occupancies" rule, QCT calculations yield a free energy value of -101 kcal/mol under these conditions, in encouraging agreement with tabulated values (-111 kcal/mol). The AIMD calculations served only to guide the definition of an effective inner-shell constraint. QCT naturally includes quantum mechanical effects that can be concerning in more primitive calculations, including electronic polarizability and induction, electron density transfer, electron correlation, molecular/atomic cooperative interactions generally, molecular flexibility, and zero-point motion. No direct assessment of the contribution of dispersion contributions to the internal energies has been attempted here, however. We anticipate that other aqueous halide ions might be treated successfully with QCT, provided that the structure of the underlying statistical mechanical theory is absorbed, i.e., that the "no split occupancies" rule is recognized.

Scaling Atomic Partial Charges of Carbonate Solvents for Lithium Ion Solvation and Diffusion.

Lithium-ion solvation and diffusion properties in ethylene carbonate (EC) and propylene carbonate (PC) were studied by molecular simulation, experiments, and electronic structure calculations. Studies carried out in water provide a reference for interpretation. Classical molecular dynamics simulation results are compared to ab initio molecular dynamics to assess nonpolarizable force field parameters for solvation structure of the carbonate solvents. Quasi-chemical theory (QCT) was adapted to take advantage of fourfold occupancy of the near-neighbor solvation structure observed in simulations and used to calculate solvation free energies. The computed free energy for transfer of Li+ to PC from water, based on electronic structure calculations with cluster-QCT, agrees with the experimental value. The simulation-based direct-QCT results with scaled partial charges agree with the electronic structure-based QCT values. The computed Li+/PF6- transference numbers of 0.35/0.65 (EC) and 0.31/0.69 (PC) agree well with NMR experimental values of 0.31/0.69 (EC) and 0.34/0.66 (PC) and similar values obtained here with impedance spectroscopy. These combined results demonstrate that solvent partial charges can be scaled in systems dominated by strong electrostatic interactions to achieve trends in ion solvation and transport properties that are comparable to ab initio and experimental results. Thus, the results support the use of scaled partial charges in simple, nonpolarizable force fields in future studies of these electrolyte solutions.

Statistical Analyses of Hydrophobic Interactions: A Mini-Review.

This review focuses on the striking recent progress in solving for hydrophobic interactions between small inert molecules. We discuss several new understandings. First, the inverse temperature phenomenology of hydrophobic interactions, i.e., strengthening of hydrophobic bonds with increasing temperature, is decisively exhibited by hydrophobic interactions between atomic-scale hard sphere solutes in water. Second, inclusion of attractive interactions associated with atomic-size hydrophobic reference cases leads to substantial, nontrivial corrections to reference results for purely repulsive solutes. Hydrophobic bonds are weakened by adding solute dispersion forces to treatment of reference cases. The classic statistical mechanical theory for those corrections is not accurate in this application, but molecular quasi-chemical theory shows promise. Finally, because of the masking roles of excluded volume and attractive interactions, comparisons that do not discriminate the different possibilities face an interpretive danger.

Dielectric Relaxation of Ethylene Carbonate and Propylene Carbonate from Molecular Dynamics Simulations.

Ethylene carbonate (EC) and propylene carbonate (PC) are widely used solvents in lithium (Li)-ion batteries and supercapacitors. Ion dissolution and diffusion in those media are correlated with solvent dielectric responses. Here, we use all-atom molecular dynamics simulations of the pure solvents to calculate dielectric constants and relaxation times, and molecular mobilities. The computed results are compared with limited available experiments to assist more exhaustive studies of these important characteristics. The observed agreement is encouraging and provides guidance for further validation of force-field simulation models for EC and PC solvents.

Hydration of Kr(aq) in Dilute and Concentrated Solutions.

Molecular dynamics simulations of water with both multi-Kr and single Kr atomic solutes are carried out to implement quasi-chemical theory evaluation of the hydration free energy of Kr(aq). This approach obtains free energy differences reflecting Kr-Kr interactions at higher concentrations. Those differences are negative changes in hydration free energies with increasing concentrations at constant pressure. The changes are due to a slight reduction of packing contributions in the higher concentration case. The observed Kr-Kr distributions, analyzed with the extrapolation procedure of Krüger et al., yield a modestly attractive osmotic second virial coefficient, B2 ≈ -60 cm(3)/mol. The thermodynamic analysis interconnecting these two approaches shows that they are closely consistent with each other, providing support for both approaches.

Molecular-scale hydrophobic interactions between hard-sphere reference solutes are attractive and endothermic.

The osmotic second virial coefficients, B2, for atomic-sized hard spheres in water are attractive (B2 < 0) and become more attractive with increasing temperature (ΔB2/ΔT < 0) in the temperature range 300 K ≤ T ≤ 360 K. Thus, these hydrophobic interactions are attractive and endothermic at moderate temperatures. Hydrophobic interactions between atomic-sized hard spheres in water are more attractive than predicted by the available statistical mechanical theory. These results constitute an initial step toward detailed molecular theory of additional intermolecular interaction features, specifically, attractive interactions associated with hydrophobic solutes.

Interfaces of propylene carbonate.

Propylene carbonate (PC) wets graphite with a contact angle of 31° at ambient conditions. Molecular dynamics simulations agree with this contact angle after 40% reduction of the strength of graphite-C atom Lennard-Jones interactions with the solvent, relative to the models used initially. A simulated nano-scale PC droplet on graphite displays a pronounced layering tendency and an Aztex pyramid structure for the droplet. Extrapolation of the computed tensions of PC liquid-vapor interface estimates the critical temperature of PC accurately to about 3%. PC molecules lie flat on the PC liquid-vapor surface and tend to project the propyl carbon toward the vapor phase. For close PC neighbors in liquid PC, an important packing motif stacks carbonate planes with the outer oxygen of one molecule snuggled into the positively charged propyl end of another molecule so that neighboring molecule dipole moments are approximately antiparallel. The calculated thermal expansion coefficient and the dielectric constants for liquid PC agree well with experiment. The distribution of PC molecule binding energies is closely Gaussian. Evaluation of the density of the coexisting vapor then permits estimation of the packing contribution to the PC chemical potential and that contribution is about two thirds of the magnitude of the contributions due to attractive interactions, with opposite sign.