Fluctuation in electrolyte solutions: the self energy

Surface tension of aqueous electrolyte solutions. A thermomechanical approach

We determine the surface tension of aqueous electrolyte solutions in contact with non-polar dielectric media using a thermomechanical approach, which involves deriving the stress tensor from the thermodynamic potential of an inhomogeneous fluid. To obtain the surface tension, we calculate both the normal and tangential pressures using the components of the stress tensor, recently derived by us [Y. A. Budkov and P. E. Brandyshev, J. Chem. Phys. 159, 174103 (2023)] within the framework of Wang's variational field theory. Using this approach, we derive an analytical expression for the surface tension in the linear approximation. At low ionic concentrations, this expression represents the classical Onsager-Samaras limiting law. By utilizing only one fitting parameter, which is related to the affinity of anions to the dielectric boundary, we successfully approximated experimental data on the surface tension of several aqueous electrolyte solutions. This approximation applies to both the solution-air and solution-dodecane interfaces, covering a wide range of electrolyte concentrations.

Modeling the Electric Double Layer at the Liposome Vesicle via Classical Density Functional Theory: Solution of Poisson's Equations for Curved Membranes

The electric double layer at the liposome vesicle membrane has been investigated by a modified fundamental-measure theory in the framework of the restricted primitive model. An analytical equation has been obtained for the mean electrostatic potential (MEP) by solving Poisson's equation for curved membranes. This study investigates the influence of vesicle size, membrane thickness, surface charges, and electrolyte concentration on the structure, composition, and width of electric double layers (EDLs) on the inner and outer membrane walls. Our findings indicate that a thin and denser layer of ions is formed at the concave wall of the membrane (inner wall) compared to that at the outer membrane. As expected, the width of the diffuse layer decreases with the concentration and surface charge. Also, when the surface charges on both concave and convex walls are the same, the absolute value of MEPs on the inner membrane, concave wall, is greater than that on the convex wall. We have also investigated the diffuse potential, which decreases with concentration, membrane thickness, and cavity size, whereas it increases with surface charges. As we expect, the contact density of counterions at the inner concave wall of the vesicle cavity is always greater than the corresponding value at the convex wall, whereas this trend reverses for co-ions. Also, the contact density of counterions (co-ions) at the inner wall decreases (increases) with cavity size, whereas it increases at the outer wall (decreases). Finally, depletion of co-ions occurs at the membrane walls with enhancement in surface charges.

A Unified Description of Salt Effects on the Liquid-Liquid Phase Separation of Proteins

Protein aggregation via liquid-liquid phase separation (LLPS) is ubiquitous in nature and is intimately connected to many human diseases. Although it is widely known that the addition of salt has crucial impacts on the LLPS of proteins, full understanding of the salt effects remains an outstanding challenge. Here, we develop a molecular theory that systematically incorporates the self-consistent field theory for charged macromolecules into the solution thermodynamics. The electrostatic interaction, hydrophobicity, ion solvation, and translational entropy are included in a unified framework. Our theory fully captures the long-standing puzzles of the nonmonotonic salt concentration dependence and the specific ion effect. We find that proteins show salting-out at low salt concentrations due to ionic screening. The solubility follows the inverse Hofmeister series. In the high salt concentration regime, protein continues salting-out for small ions but turns to salting-in for larger ions, accompanied by the reversal of the Hofmeister series. We reveal that the solubility at high salt concentrations is determined by the competition between the solvation energy and translational entropy of the ion. Furthermore, we derive an analytical criterion for determining the boundary between the salting-in and salting-out regimes, which is in good agreement with experimental results for various proteins and salt ions.

Nature of Overcharging and Charge Inversion in Electrical Double Layers

Understanding overcharging and charge inversion is one of the long-standing challenges in soft matter and biophysics. To study these phenomena, we employ the modified Gaussian renormalized fluctuation theory, which allows for the self-consistent accounting of spatially varying ionic strength as well as the spatial variations in dielectric permittivity and excluded volume effects. The underlying dependence of overcharging on the electrostatic coupling is elucidated by varying the surface charge, counterion valency, and dielectric contrast. Consistent with simulations, three characteristic regimes corresponding to weak, moderate, and strong coupling are identified. Important features like the inversion of zeta potential, crowding, and ionic layering at the surface are successfully captured. For weak coupling, there is no overcharging. In the moderate coupling regime, overcharging increases with the surface charge. Finally, in the strong coupling regime, ionic crowding and saturation in overcharging are observed. Our theory predicts a nonmonotonic dependence of charge inversion on multivalent salt concentration as well as the addition of monovalent salt, in quantitative agreement with experiments.

Ion correlation-driven like-charge attraction in multivalent salt solutions

The electrostatic double layer force is key to determining the stability and self-assembly of charged colloids and many other soft matter systems. Fully understanding the attractive force between two like-charged surfaces remains a great challenge. Here, we apply the modified Gaussian renormalized fluctuation theory to study ion correlation-driven like-charge attraction in multivalent salt solutions. The effects of spatially varying ion correlations on the structure of overlapping double layers and their free energy are self-consistently accounted for. In the presence of multivalent salts, increasing surface charge or counterion valency leads to a short-range attraction. We demonstrate that although both overcharging and like-charge attraction are outcomes of ion correlation, there is no causal relationship between them. Our theory also captures the non-monotonic dependence of like-charge attraction on multivalent salt concentration. The reduction of attraction at high salt concentrations could be a contributing factor toward the reentrant stability of charged colloidal suspensions. Our theoretical predictions are consistent with the observations reported in experiments and simulations.

Variational field theory of macroscopic forces in coulomb fluids

Based on the variational field theory framework, we extend our previous mean-field formalism [Y. A. Budkov and A. L. Kolesnikov, JStatMech 2022, 053205.2022], taking into account the electrostatic correlations of the ions. We employ a general covariant approach and derive a total stress tensor that considers the electrostatic correlations of ions. This is accomplished through an additional term that depends on the autocorrelation function of the local electric field fluctuations. Utilizing the derived total stress tensor and applying the mechanical equilibrium condition, we establish a general expression for the disjoining pressure of the Coulomb fluids, confined in a pore with a slit-like geometry. Using this equation, we derive an asymptotic expression for the disjoining pressure in a slit-like pore with non-electrified conductive walls. Present theory is the basis for future modeling of the mechanical stresses that occur in electrode pores with conductive charged walls, immersed in liquid phase electrolytes beyond the mean-field theory.

Assessment of the partial saddle point approximation in field-theoretic polymer simulations

Field-theoretic simulations are numerical treatments of polymer field theory models that go beyond the mean-field self-consistent field theory level and have successfully captured a range of mesoscopic phenomena. Inherent in molecularly-based field theories is a "sign problem" associated with complex-valued Hamiltonian functionals. One route to field-theoretic simulations utilizes the complex Langevin (CL) method to importance sample complex-valued field configurations to bypass the sign problem. Although CL is exact in principle, it can be difficult to stabilize in strongly fluctuating systems. An alternate approach for blends or block copolymers with two segment species is to make a "partial saddle point approximation" (PSPA) in which the stiff pressure-like field is constrained to its mean-field value, eliminating the sign problem in the remaining field theory, allowing for traditional (real) sampling methods. The consequences of the PSPA are relatively unknown, and direct comparisons between the two methods are limited. Here, we quantitatively compare thermodynamic observables, order-disorder transitions, and periodic domain sizes predicted by the two approaches for a weakly compressible model of AB diblock copolymers. Using Gaussian fluctuation analysis, we validate our simulation observations, finding that the PSPA incorrectly captures trends in fluctuation corrections to certain thermodynamic observables, microdomain spacing, and location of order-disorder transitions. For incompressible models with contact interactions, we find similar discrepancies between the predictions of CL and PSPA, but these can be minimized by regularization procedures such as Morse calibration. These findings mandate caution in applying the PSPA to broader classes of soft-matter models and systems.

The effect of monomer polarizability on the stability and salt partitioning in model coacervates

Coacervation of charged polymer chains has been a topic of major interest in both polymer and biological sciences, as it is a subset of a phenomenon called liquid-liquid phase separation (LLPS). In this process the polymer-rich phase separates from the polymer-lean supernatant while still maintaining its liquid-like properties. LLPS has been shown to play a crucial role in cellular homeostasis by driving the formation of membraneless organelles. It also has the potential to be harnessed to aid in novel therapeutical applications. Recent studies have demonstrated that there is no one simple mechanism which drives LLPS, which is instead a result of the combined effect of electrostatic, dipolar, hydrophobic, and other weak interactions. Using coarse-grained polymer simulations we investigate the relatively unexplored effects of monomer polarizability and spatially varying dielectric constant on LLPS propensity, and these factors affect the properties of the resulting condensates. In order to produce spatial variations in the dielectric constant, all our simulations include explicit solvent and counterions. We demonstrate that polarizability has only a minor effect on the bulk behaviour of the condensates but plays a major role when ion partitioning and microstructure are considered. We observe that the major contribution comes from the nature of the neutral blocks as endowing them with an induced dipole changes their character from hydrophobic to hydrophilic. We hypothesize that the results of this work can aid in guiding future studies concerned with LLPS by providing a general framework and by highlighting important factors which influence LLPS.

Electrostatics-Induced Nucleated Conformational Transition of Protein Aggregation

Despite the wide existence of protein aggregation in nature and its intimate connection to many human diseases, the underlying mechanism remains unclear. Here, we develop a molecular theory by systematically incorporating the self-consistent field theory for charged macromolecules into the dilute solution thermodynamics. The kinetic pathway is tracked without any restriction on the morphology of the aggregates. We find that protein aggregation at low salt concentrations is via a two-step nucleated process involving a conformational transition from metastable spherical oligomer to elongated fibril. The scaling analysis elucidates the electrostatic origin of the conformational transition: the fibril enters the screening region much earlier than the spherical aggregate. As salt concentration increases, the classical mode of one-step nucleation corresponding to macroscopic liquid-liquid phase separation is recovered. Our results reveal that the screened electrostatic interaction is essential for the existence of the metastable oligomer and its subsequent conformational transition to fibril. The theoretical predictions of the kinetic pathway and the morphology of the aggregates are in good agreement with the experimental observations of real proteins.

Mathematical modeling of microscale biology: Ion pairing, spatially varying permittivity, and Born energy in glycosaminoglycan brushes

Biological macromolecules including nucleic acids, proteins, and glycosaminoglycans are typically anionic and can span domains of up to hundreds of nanometers and even micron length scales. The structures exist in crowded environments that are dominated by multivalent electrostatic interactions that can be modeled using mean-field continuum approaches that represent underlying molecular nanoscale biophysics. We develop such models for glycosaminoglycan brushes using steady state modified Poisson-Boltzmann models that incorporate important ion-specific (Hofmeister) effects. The results quantify how electroneutrality is attained through ion electrophoresis, spatially-varying permittivity hydration forces, and ion-specific pairing. Brush-salt interfacial profiles of the electrostatic potential as well as bound and unbound ions are characterized for imposed jump conditions across the interface. The models should be applicable to many intrinsically-disordered biophysical environments and are anticipated to provide insight into the design and development of therapeutics and drug-delivery vehicles to improve human health.

Numerical Techniques for Applications of Analytical Theories to Sequence-Dependent Phase Separations of Intrinsically Disordered Proteins

Biomolecular condensates, physically underpinned to a significant extent by liquid-liquid phase separation (LLPS), are now widely recognized by numerous experimental studies to be of fundamental biological, biomedical, and biophysical importance. In the face of experimental discoveries, analytical formulations emerged as a powerful yet tractable tool in recent theoretical investigations of the role of LLPS in the assembly and dissociation of these condensates. The pertinent LLPS often involves, though not exclusively, intrinsically disordered proteins engaging in multivalent interactions that are governed by their amino acid sequences. For researchers interested in applying these theoretical methods, here we provide a practical guide to a set of computational techniques devised for extracting sequence-dependent LLPS properties from analytical formulations. The numerical procedures covered include those for the determination of spinodal and binodal phase boundaries from a general free energy function with examples based on the random phase approximation in polymer theory, construction of tie lines for multiple-component LLPS, and field-theoretic simulation of multiple-chain heteropolymeric systems using complex Langevin dynamics. Since a more accurate physical picture often requires comparing analytical theory against explicit-chain model predictions, a commonly utilized methodology for coarse-grained molecular dynamics simulations of sequence-specific LLPS is also briefly outlined.

Excluded-Volume Interactions in Field-Theoretic Simulations: Multiconvolutions and Model Equivalence

To deal with divergences of functional integrals in field-theoretic simulations (FTS) of complex fluids, the microscopic density is often smeared by being replaced by a convoluted one, typically using a Gaussian masking function. The smearing changes radically the nature of nonbonded interactions of the original microscopic density and results in a regularized model that is free of ultraviolet (UV) divergences. In this work, we first resolve a few fundamental issues related with the use of masking functions for δ-interactions in FTS and then we detail a new methodology that builds on the concept of multiconvoluted inverse potentials and a principle of model equivalence for statistical weights to accommodate more physically relevant interactions in FTS. The capabilities of the new approach are highlighted by examining the Gaussian-regularized Edwards model (GREM) and the Yukawa potential. A successful test calculation of the excess chemical potential of a polymer chain in a good solvent with the GREM illustrates the power of the new theoretical framework.

Self-Consistent Description of Vapor-Liquid Interface in Ionic Fluids

Inhomogeneity of ion correlation widely exists in many physicochemical, soft matter, and biological systems. Here, we apply the modified Gaussian renormalized fluctuation theory to study the classic example of the vapor-liquid interface of ionic fluids. The ion correlation is decomposed into a short-range contribution associated with the local electrostatic environment and a long-range contribution accounting for the spatially varying ionic strength and dielectric permittivity. For symmetric salt, both the coexistence curve and the interfacial tension predicted by our theory are in quantitative agreement with simulation data reported in the literature. Furthermore, we provide the first theoretical prediction of interfacial structure for asymmetric salt, highlighting the importance of capturing local charge separation.

Analytical Formulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range Interactions

A theory for sequence-dependent liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) in the study of biomolecular condensates is formulated by extending the random phase approximation (RPA) and field-theoretic simulation (FTS) of heteropolymers with spatially long-range Coulomb interactions to include the fundamental effects of short-range, hydrophobic-like interactions between amino acid residues. To this end, short-range effects are modeled by Yukawa interactions between multiple nonelectrostatic charges derived from an eigenvalue decomposition of pairwise residue-residue contact energies. Chain excluded volume is afforded by incompressibility constraints. A mean-field approximation leads to an effective Flory-Huggins χ parameter, which, in conjunction with RPA, accounts for the contact-interaction effects of amino acid composition and the sequence-pattern effects of long-range electrostatics in IDP LLPS, whereas FTS based on the formulation provides full sequence dependence for both short- and long-range interactions. This general approach is illustrated here by applications to variants of a natural IDP in the context of several different amino-acid interaction schemes as well as a set of different model hydrophobic-polar sequences sharing the same composition. Effectiveness of the methodology is verified by coarse-grained explicit-chain molecular dynamics simulations.

Electrostatic Correlation Induced Ion Condensation and Charge Inversion in Multivalent Electrolytes

The study of the electrical double layer lies at the heart of colloidal and interfacial science. The standard mean-field Poisson-Boltzmann (PB) theory is incapable of modeling many phenomena originating from ion correlation. An important example is charge inversion or overcharging of electrical double layers in multivalent electrolyte solutions. Existing theories aiming to include correlations cannot capture the non-monotonic dependence of charge inversion on salt concentration because they have not systematically accounted for the inhomogeneous nature of correlations from surface to the bulk and the excluded volume effect of ions and solvent molecules. In this work, we modify the Gaussian renormalized fluctuation theory by including the excluded volume effect to study ion condensation and charge inversion. A boundary layer approach is developed to accurately model the giant difference in ion correlations between the condensed layer near the surface and the diffuse layer outside. The theory is used to study charge inversion in multivalent electrolytes and their mixtures. We predict a surface charge induced formation of a three-dimensional condensed layer, which is necessary but not sufficient for charge inversion. The value of the effective surface potential is found to depend non-monotonically on the bulk salt concentration. Our results also show a non-monotonic reduction in charge inversion in monovalent and multivalent electrolyte mixtures. Our work is the first to qualitatively reproduce experimental and simulation observations and explains the underlying physics.

Association of two polyelectrolytes in salt solutions

The association of polyelectrolytes (PEs) in solution affects a wealth of structural and dynamic behaviors, and is also fundamentally important for an understanding of protein association and aggregation. Here, we theoretically study the association of two PE chains by addressing the stability and morphology of the non-spherical associates. Our theory predicts that an elongated pearl-necklace (PN) associate can be stable at high salt concentrations due to the screened electrostatic repulsion. This contradicts the implication of scaling theory. In addition, there is no one-to-one correspondence between the morphology of the associate and its constituting unimers, which is demonstrated by the existence of different association modes.

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

Electric-field-induced oscillations in ionic fluids: a unified formulation of modified Poisson-Nernst-Planck models and its relevance to correlation function analysis

We theoretically investigate an electric-field-driven system of charged spheres as a primitive model of concentrated electrolytes under an applied electric field. First, we provide a unified formulation for the stochastic charge and density dynamics of the electric-field-driven primitive model using the stochastic density functional theory (DFT). The stochastic DFT integrates the four frameworks (the equilibrium and dynamic DFTs, the liquid state theory and the field-theoretic approach), which allows us to justify in a unified manner various modifications previously made for the Poisson-Nernst-Planck model. Next, we consider stationary density-density and charge-charge correlation functions of the primitive model with a static electric field. We predict an electric-field-induced synchronization between emergences of density and charge oscillations. We are mainly concerned with the emergence of stripe states formed by segregation bands transverse to the external field, thereby demonstrating the following: (i) the electric-field-induced crossover occurs prior to the conventional Kirkwood crossover without an applied electric field, and (ii) the ion concentration dependence of the decay lengths at the onset of oscillations bears a similarity to the underscreening behavior found by recent simulation and theoretical studies on equilibrium electrolytes. Also, the 2D inverse Fourier transform of the correlation function illustrates the existence of stripe states beyond the electric-field-induced Kirkwood crossover.

Field theory description of ion association in re-entrant phase separation of polyampholytes

Phase separation of several different overall neutral polyampholyte species (with zero net charge) is studied in solution with two oppositely charged ion species that can form ion-pairs through an association reaction. A field theory description of the system, that treats polyampholyte charge sequence dependent electrostatic interactions as well as excluded volume effects, is hereby given. Interestingly, analysis of the model using random phase approximation and field theoretic simulation consistently show evidence of a re-entrant polyampholyte phase separation at high ion concentrations when there is an overall decrease of volume upon ion-association. As an illustration of the ramifications of our theoretical framework, several polyampholyte concentration vs ion concentration phase diagrams under constant temperature conditions are presented to elucidate the dependence of phase separation behavior on polyampholyte sequence charge pattern as well as ion-pair dissociation constant, volumetric effects on ion association, solvent quality, and temperature.

Effects of Local Order Parameter Dependent Transport Coefficient in Diblock Copolymers Under Applied Electric Fields

We present an approach for constructing thermodynamically consistent time-dependent models relevant to thin films of diblock copolymers in applied electric fields. The approach is based on the principles of linear irreversible thermodynamics and in this work, it is applied to study the effects of electric  fields on thin  films of incompressible diblock copolymers. Enforcement of local incompressibility constraint at all times leads to a local order parameter dependent transport coefficient in the model for the diblock copolymers. The dependence of transport coefficient on the local order parameter is used to relate it with diffusion constant of Rouse chains and leads to sensitivity of the model to initial conditions. Also, transient behavior is found to be a affected when compared with an ad hoc model assuming a constant transport coefficient. Numerical results such as electric  field induced alignment of lamellae domains due to the  field are found to be in qualitative agreements with experiments.This approach opens up a systematic way of developing kinetic models for simulating effects of electrolytes added to thin  films containing diblock copolymers in the presence of applied electric  fields.

Acid-Base Equilibrium and Dielectric Environment Regulate Charge in Supramolecular Nanofibers

Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C₁₆ − V₂A₂E₂. In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment.

Theory of freezing point depression in charged porous media

Freezing in charged porous media can induce significant pressure and cause damage to tissues and functional materials. We formulate a thermodynamically consistent theory to model freezing phenomena inside charged heterogeneous porous space. Two regimes are distinguished: free ions in open pore space lead to negligible effects of freezing point depression and pressure. On the other hand, if nanofluidic salt trapping happens, subsequent ice formation is suppressed due to the high concentration of ions in the electrolyte. In this case our theory predicts that freezing starts at a significantly lower temperature compared to pure water. In one dimension, as the temperature goes even lower, ice continuously grows until the salt concentration reaches saturation, all ions precipitate to form salt crystals, and freezing completes. Enormous pressure can be generated if initial salt concentration is high before salt entrapment. We show modifications to the classical nucleation theory due to the trapped salt ions. Interestingly, although the freezing process is enormously changed by trapped salts, our analysis shows that the Gibbs-Thompson equation on confined melting point shift is not affected by the presence of the electrolyte.

Competition between Born solvation, dielectric exclusion, and Coulomb attraction in spherical nanopores

The recent measurement of a very low dielectric constant, ε, of water confined in nanometric slit pores leads us to reconsider the physical basis of ion partitioning into nanopores. For confined ions in chemical equilibrium with a bulk of dielectric constant ε_{b}>ε, three physical mechanisms, at the origin of ion exclusion in nanopores, are expected to be modified due to this dielectric mismatch: dielectric exclusion at the water-pore interface (with membrane dielectric constant, ε_{m}<ε), the solvation energy related to the difference in Debye-Hückel screening parameters in the pore, κ, and in the bulk κ_{b}, and the classical Born solvation self-energy proportional to ε^{-1}-ε_{b}^{-1}. Our goal is to clarify the interplay between these three mechanisms and investigate the role played by the Born contribution in ionic liquid-vapor (LV) phase separation in confined geometries. We first compute analytically the potential of mean force (PMF) of an ion of radius R_{i} located at the center of a nanometric spherical pore of radius R. Computing the variational grand potential for a solution of confined ions, we then deduce the partition coefficients of ions in the pore versus R and the bulk electrolyte concentration ρ_{b}. We show how the ionic LV transition, directly induced by the abrupt change of the dielectric contribution of the PMF with κ, is favored by the Born self-energy and explore the decrease of the concentration in the pore with ε both in the vapor and liquid states. Phase diagrams are established for various parameter values and we show that a signature of this phase transition can be detected by monitoring the total osmotic pressure as a function of R. For charged nanopores, these exclusion effects compete with the electrostatic attraction that imposes the entry of counterions into the pore to enforce electroneutrality. This study will therefore help in deciphering the respective roles of the Born self-energy and dielectric mismatch in experiments and simulations of ionic transport through nanopores.

Extracting free energies of counterion binding to polyelectrolytes by molecular dynamics simulations

We use all-atom molecular dynamics simulations to extract ΔG^eff, the free energy of binding of potassium ions K⁺ to the partially charged polyelectrolyte poly(acrylic acid), or PAA, in dilute regimes. Upon increasing the charge fraction of PAA, the chains adopt more extended conformations, and simultaneously, potassium ions bind more strongly (i.e., with more negative ΔG^eff) to the highly charged chains to relieve electrostatic repulsions between charged monomers along the chains. We compare the simulation results with the predictions of a model that describes potassium binding to PAA chains as a reversible reaction whose binding free energy (ΔG^eff) is adjusted from its intrinsic value (ΔG) by electrostatic correlations, captured by a random phase approximation. The bare or intrinsic binding free energy ΔG, which is an input in the model, depends on the binding species and is obtained from the radial distribution function of K⁺ around the charged monomer of a singly charged, short PAA chain in dilute solutions. We find that the model yields semi-quantitative predictions for ΔG^eff and the degree of potassium binding to PAA chains, α, as a function of PAA charge fraction without using fitting parameters.

Pimples reduce and dimples enhance flat dielectric surface image repulsion

In solid-liquid, or liquid-liquid, interfaces with dielectric contrast, charged particles interact with the induced polarization charge of the interface. These interactions contribute to an effective self-energy of the bulk ions and mediate ion-ion interactions. For flat interfaces, the self-energy and the mediated interactions are neatly constructed by the image charge method. For other geometries, explicit results are scarce and the problem must be treated via approximations or direct computation. The case of interfaces with roughness is of great practical importance. This article provides analytical results, valid to first-order in perturbation theory, for the self-energy of particles near rough substrates. Explicit formulas are provided for the case of a sinusoidal deformation of a flat surface. Generic deformations can be treated by superposition. In addition to results for the self-energy, the surface polarization charge is presented as a quadrature. The interaction between an ion and the deformed surface is modified by the change in relative distance as well as by the local curvature of the surface. Solid walls, with a lower dielectric constant than the liquid, repel all ions. We show that the repulsion is reduced by local convexity and enhanced by concavity; dimples are more repulsive than pimples.

Ion transport in electrolytes of dielectric nanodevices

Ion transport in electrolytes with nanoscale confinements is of great importance in many fields such as nanofluidics and electrochemical energy devices. The mobility and conductance for ions are often described by the classical Debye-Hückel-Onsager (DHO) theory but this theory fails for ions near dielectric interfaces. We propose a generalized DHO theory by using the Wentzel-Kramers-Brillouin techniques for the solution of the Onsager-Fuoss equation with variable coefficients. The theory allows to quantitatively measure physical quantities of ion transport in nanodevices and is demonstrated to well explain the abnormal increase or decrease of the ionic mobility tuned via the dielectric mismatch. By numerical calculations, our theory unravels the crucial role of the size of confinements and the ionic concentration on the ion transport, and demonstrates that the dielectric polarization can provide a giant enhancement on the conductance of electrolytes in nanodevices. This mechanism provides a practical guide for related nanoscale technologies with controllable transport properties.

Field-theoretic simulations beyond δ-interactions: Overcoming the inverse potential problem in auxiliary field models

Modern field-theoretic simulations of complex fluids and polymers are constructed around a particle-to-field transformation that brings an inverse potential u^-1 in the model equations. This has restricted the application of the framework to systems characterized by relatively simple pairwise interatomic interactions; for example, excluded volume effects are treated through the use of δ-function interactions. In this study, we first review available nonbonded pair interactions in field-theoretic models and propose a classification. Then, we outline the inverse potential problem and present an alternative approach on the basis of a saddle-point approximation, enabling the use of a richer set of pair interaction functions. We test our approach by using as an example the Morse potential, which finds extensive applications in particle-based simulations, and we calibrate u^-1 with results from a molecular dynamics simulation. The u^-1 thus obtained is consistent with the field-theoretic model equations, and when used in stand-alone self-consistent field simulations, it produces the correct fluid structure starting from a random initial state of the density field.

Looping-in complexation and ion partitioning in nonstoichiometric polyelectrolyte mixtures

A wide variety of intracellular membraneless compartments are formed via liquid-liquid phase separation of charged proteins and nucleic acids. Understanding the stability of these compartments, while accounting for the compositional heterogeneity intrinsic to cellular environments, poses a daunting challenge. We combined experimental and theoretical efforts to study the effects of nonstoichiometric mixing on coacervation behavior and accurately measured the concentrations of polyelectrolytes and small ions in the coacervate and supernatant phases. For synthetic polyacrylamides and polypeptides/DNA, with unequal mixing stoichiometry, we report a general "looping-in" phenomenon found around physiological salt concentrations, where the polymer concentrations in the coacervate initially increase with salt addition before subsequently decreasing. This looping-in behavior is captured by a molecular model that considers reversible ion binding and electrostatic interactions. Further analysis in the low-salt regime shows that the looping-in phenomenon originates from the translational entropy of counterions that are needed to neutralize nonstoichiometric coacervates.

Weakening of Solvation-Induced Ordering by Composition Fluctuation in Salt-Doped Block Polymers

The spontaneous ordering of block polymers doped with ions is affected by both selective solvation and long-range Coulombic interaction. The mean-field treatment was recently shown to overestimate the solvation-induced ordering, requiring a large solvation radius to fit experimental phase diagrams, which may be relieved by including composition fluctuations. Treating the composition fluctuations in such systems is challenging because of the need of resolving heterogeneous dielectric profile that couples with the ordering itself. Starting from a minimal model, we develop a Landau-Brazovskiĭ expansion for the free energy of salt-doped block polymer near the ordering transition. It is found that the wavelength for typical composition fluctuations first decreases with salt doping, due to Coulombic interaction, then increases due to ionic solvation. Two mechanisms that weaken the solvation-enhanced ordering are identified: the Brazovskiĭ-type composition fluctuation that stabilizes disordered phase, and the coupling between mismatch in dispersion interaction and the dielectric permittivity through monomeric polarizability.

Molecularly Informed Field Theories from Bottom-up Coarse-Graining

Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.

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.

A Perspective on the Design of Ion-Containing Polymers for Polymer Electrolyte Applications

Ion-containing polymers have numerous potential applications as energy storage and conversion devices, water purification membranes, and gas separation membranes, to name a few. Given the low dielectric constant of the media, ions and charges on polymers in a molten state interact strongly producing large effects on chain statistics, thermodynamics, and diffusion properties. Here, we discuss recent research accomplishments on the effects of ionic correlation and dielectric heterogeneity on the phase behavior of ion-containing polymers. Progress made in studying ion transport properties in these material systems is also highlighted. Charged block copolymers (BCPs), among all kinds of ion-containing polymers, have a particular advantage owing to their robust mechanical support and ion conducting paths provided by the segregation of the neutral and charged blocks. Coulombic interactions among the charges play a critical role in determining the phase segregation in charged BCPs and the domain size of charge-rich regions. We show that strongly charged BCPs display ordered phases as a result of electrostatic interactions alone. In addition, bulky charge-containing side groups attached to the charged block lead to the formation of morphologies that provide continuous channels and better dissociation for ion conduction purposes. Finally, a few avenues for designing ion-containing polymers for energy applications are discussed.

Subcompartmentalization of polyampholyte species in organelle-like condensates is promoted by charge-pattern mismatch and strong excluded-volume interaction

Polyampholyte field theory and explicit-chain molecular dynamics models of sequence-specific phase separation of a system with two intrinsically disordered protein (IDP) species indicate consistently that a substantial polymer excluded volume and a significant mismatch of the IDP sequence charge patterns can act in concert, but not in isolation, to demix the two IDP species upon condensation. This finding reveals an energetic-geometric interplay in a stochastic, "fuzzy" molecular recognition mechanism that may facilitate subcompartmentalization of membraneless organelles.

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.

Salt-Induced Liquid-Liquid Phase Separation: Combined Experimental and Theoretical Investigation of Water-Acetonitrile-Salt Mixtures

Salt-induced liquid-liquid phase separation in liquid mixtures is a common phenomenon in nature and in various applications, such as in separation and extraction of chemicals. Here, we present results of a systematic investigation of the phase behaviors in water-acetonitrile-salt mixtures using a combination of experiment and theory. We obtain complete ternary phase diagrams for nine representative salts in water-acetonitrile mixtures by cloud point and component analysis. We construct a thermodynamic free energy model by accounting for the nonideal mixing of the liquids, ion hydration, electrostatic interactions, and Born energy. Our theory yields phase diagrams in good agreement with the experimental data. By comparing the contributions due to the electrostatic interaction, Born energy, and hydration, we find that hydration is the main driving force for the liquid-liquid separation and is a major contributor to the specific ion effects. Our theory highlights the important role of entropy in the hydration driving force. We discuss the implications of our findings in the context of salting-out assisted liquid-liquid extraction and make suggestions for selecting salt ions to optimize the separation performance.

Comparative roles of charge, π, and hydrophobic interactions in sequence-dependent phase separation of intrinsically disordered proteins

Endeavoring toward a transferable, predictive coarse-grained explicit-chain model for biomolecular condensates underlain by liquid-liquid phase separation (LLPS) of proteins, we conducted multiple-chain simulations of the N-terminal intrinsically disordered region (IDR) of DEAD-box helicase Ddx4, as a test case, to assess roles of electrostatic, hydrophobic, cation-π, and aromatic interactions in amino acid sequence-dependent LLPS. We evaluated three different residue-residue interaction schemes with a shared electrostatic potential. Neither a common hydrophobicity scheme nor one augmented with arginine/lysine-aromatic cation-π interactions consistently accounted for available experimental LLPS data on the wild-type, a charge-scrambled, a phenylalanine-to-alanine (FtoA), and an arginine-to-lysine (RtoK) mutant of Ddx4 IDR. In contrast, interactions based on contact statistics among folded globular protein structures reproduce the overall experimental trend, including that the RtoK mutant has a much diminished LLPS propensity. Consistency between simulation and experiment was also found for RtoK mutants of P-granule protein LAF-1, underscoring that, to a degree, important LLPS-driving π-related interactions are embodied in classical statistical potentials. Further elucidation is necessary, however, especially of phenylalanine's role in condensate assembly because experiments on FtoA and tyrosine-to-phenylalanine mutants suggest that LLPS-driving phenylalanine interactions are significantly weaker than posited by common statistical potentials. Protein-protein electrostatic interactions are modulated by relative permittivity, which in general depends on aqueous protein concentration. Analytical theory suggests that this dependence entails enhanced interprotein interactions in the condensed phase but more favorable protein-solvent interactions in the dilute phase. The opposing trends lead to only a modest overall impact on LLPS.

Modulating Microphase Separation of Lamellae-Forming Diblock Copolymers via Ionic Junctions

We present a molecular dynamics simulation study investigating the phase behavior of lamellae-forming diblock copolymers with a single ionic junction on the backbone. Our results show qualitative agreement with experimental findings regarding enhanced microphase separation with the introduction of an ionic junction at the conjunction point, while further revealing nonmonotonic changes in domain spacing and order-disorder transition as a function of the electrostatic interaction strength. This highlights the dominant roles of entropic and binding effects of counterions under weak and strong ionic correlations, respectively. The location of the ionic junction is found to effectively modulate the charge distribution and chain conformation in the ordered domains; its presence in the middle of a block promotes folding of the block, leading to a smaller domain size. These findings demonstrate the interplay of ionic coupling with steric hindrance and chain end effects, which enhances our understanding of the delicate control over the microphase domain features.

Enhancement of the ion concentration in a salt solution near a wall due to electrical image potentials and enhancement of surface tension due to the presence of salt

Effects of electrical image potentials on the salt ion concentration near a solid wall are studied using a one-loop approximation treatment of the grand canonical partition function, which is the Debye-Hückel approximation. Electrical image potentials resulting from both metallic and dielectric walls of dielectric constant larger than that of water near the wall are considered. Our treatment of this problem supports the conclusions of an earlier publication by one of the authors which shows that near a solid wall there should be a high concentration of ions, resulting from image potentials. We have also applied our treatment to the increase of the surface tension of a liquid that occurs when salt is dissolved in the liquid. Our treatment gives the -c_{s}log(c_{s}) dependence of the surface tension found by Onsager and Samarasa in the small c_{s} limit, where c_{s} is the salt concentration.