Contrasting Dielectric Properties of Electrolyte Solutions with Polar and Polarizable Solvents

Predicting surfactant phase behavior with a molecularly informed field theory

HYPOTHESIS

The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations.

SIMULATION WORK

We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements.

FINDINGS

The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.

Understanding specific ion effects and the Hofmeister series

Specific ion effects (SIE), encompassing the Hofmeister Series, have been known for more than 130 years since Hofmeister and Lewith's foundational work. SIEs are ubiquitous and are observed across the medical, biological, chemical and industrial sciences. Nevertheless, no general predictive theory has yet been able to explain ion specificity across these fields; it remains impossible to predict when, how, and to what magnitude, a SIE will be observed. In part, this is due to the complexity of real systems in which ions, counterions, solvents and cosolutes all play varying roles, which give rise to anomalies and reversals in anticipated SIEs. Herein we review the historical explanations for SIE in water and the key ion properties that have been attributed to them. Systems where the Hofmeister series is perturbed or reversed are explored, as is the behaviour of ions at the liquid-vapour interface. We discuss SIEs in mixed electrolytes, nonaqueous solvents, and in highly concentrated electrolyte solutions - exciting frontiers in this field with particular relevance to biological and electrochemical applications. We conclude the perspective by summarising the challenges and opportunities facing this SIE research that highlight potential pathways towards a general predictive theory of SIE.

A Simple Explicit-Solvent Model of Polyampholyte Phase Behaviors and Its Ramifications for Dielectric Effects in Biomolecular Condensates

Biomolecular condensates such as membraneless organelles, underpinned by liquid-liquid phase separation (LLPS), are important for physiological function, with electrostatics, among other interaction types, being a prominent force in their assembly. Charge interactions of intrinsically disordered proteins (IDPs) and other biomolecules are sensitive to the aqueous dielectric environment. Because the relative permittivity of protein is significantly lower than that of water, the interior of an IDP condensate is expected to be a relatively low-dielectric regime, which aside from its possible functional effects on client molecules should facilitate stronger electrostatic interactions among the scaffold IDPs. To gain insight into this LLPS-induced dielectric heterogeneity, addressing in particular whether a low-dielectric condensed phase entails more favorable LLPS than that posited by assuming IDP electrostatic interactions are uniformly modulated by the higher dielectric constant of the pure solvent, we consider a simplified multiple-chain model of polyampholytes immersed in explicit solvents that are either polarizable or possess a permanent dipole. Notably, simulated phase behaviors of these systems exhibit only minor to moderate differences from those obtained using implicit-solvent models with a uniform relative permittivity equals to that of pure solvent. Buttressed by theoretical treatments developed here using random phase approximation and polymer field-theoretic simulations, these observations indicate a partial compensation of effects between favorable solvent-mediated interactions among the polyampholytes in the condensed phase and favorable polyampholyte-solvent interactions in the dilute phase, often netting only a minor enhancement of overall LLPS propensity from the very dielectric heterogeneity that arises from the LLPS itself. Further ramifications of this principle are discussed.

Like dissolves like: A first-principles theory for predicting liquid miscibility and mixture dielectric constant

Liquid mixtures are ubiquitous. Miscibility and dielectric constant are fundamental properties that govern the applications of liquid mixtures. However, despite their importance, miscibility is usually predicted qualitatively based on the vaguely defined polarity of the liquids, and the dielectric constant of the mixture is modeled by introducing mixing rules. Here, we develop a first-principles theory for polar liquid mixtures using a statistical field approach, without resorting to mixing rules. With this theory, we obtain simple expressions for the mixture's dielectric constant and free energy of mixing. The dielectric constant predicted by this theory agrees well with measured data for simple binary mixtures. On the basis of the derived free energy of mixing, we can construct a miscibility map in the parameter space of the dielectric constant and molar volume for each liquid. The predicted miscibility shows remarkable agreement with known data, thus providing a quantitative basis for the empirical "like-dissolves-like" rule.

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. The theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the "accessibility" of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

Statistical field theory of ion–molecular solutions

Schematic representation of the multipolar molecule surrounded by salt ions in a dielectric solvent medium.

Field-Theoretic Study of Salt-Induced Order and Disorder in a Polarizable Diblock Copolymer

We study a salt-doped polarizable symmetric diblock copolymer using a recently developed field theory that self-consistently embeds dielectric response, ion solvation energies, and van der Waals (vdW) attractions via the incorporation of segment polarizabilities and fixed dipoles. This field theory is amenable to direct simulation via the complex Langevin sampling technique and, thus, requires no approximations beyond the phenomenology of the underlying molecular model. We measure the shift in the order-disorder transition (ODT) of a diblock copolymer with salt-loading in field-theoretic simulations and observe rich behavior in which solvation, dilution and charge screening effects compete to determine whether the ordered or disordered phase is stabilized. At low salt concentrations, the salt behaves as a selective solvent, localizing into the high-dielectric domains and stabilizing the ordered phase. At high salt concentrations, however, the salt localization vanishes due to charge screening effects, and the salt behaves as a nonselective solvent that screens vdW attractions and stabilizes the disordered phase.