Coarse-Grained Model of the Dynamics of Electrolyte Solutions

How specific ion effects influence the mechanical behaviors of amide macromolecules? A cross-scale study

The mechanisms of specific ion effects on the properties of amide macromolecules is essential to understanding the evolution of life. Because most biological macromolecules contain both complex hydrophilic and hydrophobic structures, it is challenging to accurately identify the contributions of molecular structure to macroscopic behaviors. Herein, we investigated the influence of specific ion effects on the mechanical behaviors of poly(N-isopropylacrylamide) and neutral polyacrylamide (i.e., PNIPAM and NPAM), through a cross-scale study that includes single-molecule force spectroscopy, molecular dynamics simulation and macro mechanical method. The results indicate that the molecular conformation can be markedly influenced by the hydrophilicity (or hydrophobicity) of both macromolecule chain and ions. An extended chain conformation can be obtained when the side groups and ions are relatively hydrophilic, which can also increase the elasticity of a macromolecule chain and film materials. The relatively hydrophobic components promote the collapse of macromolecule chains and reduce the molecular elasticity. It is believed that the hydrogen bonding intensity between a macromolecule chain and aquated ions controls the chain conformation and the elasticity of molecules and films. This study is not only helpful for understanding the self-assembly mechanism of organisms but also provides a way to associate the molecular properties with the macroscopic performance of materials.

Do Specific Ion Effects on Collective Relaxation Arise from Perturbation of Hydrogen-Bonding Network Structure?

The change in the transport properties (i.e., water diffusivity, shear viscosity, etc.) when adding salts to water has been used to classify ions as either being chaotropic or kosmotropic, a terminology based on the presumption that this phenomenon arises from respective breakdown or enhancement of the hydrogen-bonding network structure. Recent quasi-elastic neutron scattering measurements of the collective structural relaxation time, τC, in aqueous salt solutions were interpreted as confirming this proposed origin of ion effects on the dynamics of water. However, we find similar changes in τC in the same salt solutions based on molecular dynamics (MD) simulations using a coarse-grained water model in which no hydrogen bonding exists, challenging this conventional interpretation of mobility change resulting from the addition of salts to water. A thorough understanding of specific ion effects should be useful in diverse material manufacturing and biomedical applications, where these effects are prevalent, but poorly understood.

Molecular Dynamics Simulation of the Influence of Temperature and Salt on the Dynamic Hydration Layer in a Model Polyzwitterionic Polymer PAEDAPS

We investigate the hydration of poly(3-[2-(acrylamido) ethyldimethylammonio] propanesulfonate) over a range of temperatures in pure water and with the inclusion of 0.1 mol/L NaCl using atomistic molecular dynamics simulation. Drawing on concepts drawn from the field of glass-forming liquids, we use the Debye-Waller parameter () for describing the water mobility gradient around the polybetaine backbone extending to an overall distance ≈18 Å. The water mobility in this layer is defined through the mean-square water molecule displacement at a time on the order of water's β-relaxation time. The brushlike topology of polybetaines leads to two regions in the dynamic hydration layer. The inner region of ≈10.5 Å is explored by pendant group conformational motions, and the outer region of ≈7.5 Å represents an extended layer of reduced water mobility relative to bulk water. The dynamic hydration layer extends far beyond the static hydration layer, adjacent to the polymer.

Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions

Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.

The Contribution of the Ion-Ion and Ion-Solvent Interactions in a Molecular Thermodynamic Treatment of Electrolyte Solutions

Developing molecular equations of state to treat electrolyte solutions is challenging due to the long-range nature of the Coulombic interactions. Seminal approaches commonly used are the mean spherical approximation (MSA) and the Debye-Hückel (DH) theory to account for ion-ion interactions and, often, the Born theory of solvation for ion-solvent interactions. We investigate the accuracy of the MSA and DH approaches using each to calculate the contribution of the ion-ion interactions to the chemical potential of NaCl in water, comparing these with newly computer-generated simulation data; the ion-ion contribution is isolated by selecting an appropriate primitive model with a Lennard-Jones force field to describe the solvent. A study of mixtures with different concentrations and ionic strengths reveals that the calculations from both MSA and DH theories are of similar accuracy, with the MSA approach resulting in marginally better agreement with the simulation data. We also demonstrate that the Born theory provides a good qualitative description of the contribution of the ion-solvent interactions; we employ an explicitly polar water model in these simulations. Quantitative agreement up to moderate salt concentrations and across the relevant range of temperature is achieved by adjusting the Born radius using simulation data of the free energy of solvation. We compute the radial and orientational distribution functions of the systems, thereby providing further insight on the differences observed between the theory and simulation. We thus provide rigorous benchmarks for use of the MSA, DH, and Born theories as perturbation approaches, which will be of value for improving existing models of electrolyte solutions, especially in the context of equations of state.

Prestressed Composite Polymer Gels as a Model of the Extracellular-Matrix of Cartilage

Articular cartilage is a composite hydrogel found in animal and human joints, which exhibits unique load-bearing properties that have been challenging to reproduce in synthetic materials and model in molecular dynamics (MD) simulations. We computationally investigate a composite hydrogel that mimics key functional properties of articular cartilage as a potential biomimetic model to investigate its unique load-bearing properties. Specifically, we find that the emergence of prestress in composite gels derives primarily from the stiffness of the polymer matrix and the asymmetry in the enthalpic interactions of the embedded particles and polymer matrix. Our MD simulations of the development of prestress agree qualitatively with osmotic pressure measurements observed in our model composite hydrogel material.

Phosphorylation and sulfation share a common biosynthetic pathway, but extend biochemical and evolutionary diversity of biological macromolecules in distinct ways

Phosphate and sulfate groups are integral to energy metabolism and introduce negative charges into biological macromolecules. One purpose of such modifications is to elicit precise binding/activation of protein partners. The physico-chemical properties of the two groups, while superficially similar, differ in one important respect—the valency of the central (phosphorus or sulfur) atom. This dictates the distinct properties of their respective esters, di-esters and hence their charges, interactions with metal ions and their solubility. These, in turn, determine the contrasting roles for which each group has evolved in biological systems. Biosynthetic links exist between the two modifications; the sulfate donor 3′-phosphoadenosine-5′-phosphosulfate being formed from adenosine triphosphate (ATP) and adenosine phosphosulfate, while the latter is generated from sulfate anions and ATP. Furthermore, phosphorylation, by a xylosyl kinase (Fam20B, glycosaminoglycan xylosylkinase) of the xylose residue of the tetrasaccharide linker region that connects nascent glycosaminoglycan (GAG) chains to their parent proteoglycans, substantially accelerates their biosynthesis. Following observations that GAG chains can enter the cell nucleus, it is hypothesized that sulfated GAGs could influence events in the nucleus, which would complete a feedback loop uniting the complementary anionic modifications of phosphorylation and sulfation through complex, inter-connected signalling networks and warrants further exploration.

Structures of cationic and anionic polyelectrolytes in aqueous solutions: the sign effect

In this study, we use molecular dynamics simulation to explore the structures of anionic and cationic polyelectrolytes in aqueous solutions. We first confirm the significantly stronger solvation effects of single anions compared to cations in water at the fixed ion radii, due to the reversal orientations of asymmetric dipolar H₂O molecules around the ions. Based on this, we demonstrate that the solvation discrepancy of cations/anions and electrostatic correlations of ionic species can synergistically cause the nontrivial structural difference between single anionic and cationic polyelectrolytes. The cationic polyelectrolyte shows an extended structure whereas the anionic polyelectrolyte exhibits a collapsed structure, and their structural differences decline with increasing the counterion size. Furthermore, we corroborate that multiple cationic polyelectrolytes or multiple anionic polyelectrolytes can exhibit largely differential molecular architectures in aqueous solutions. In the solvation dominant regime, the polyelectrolyte solutions exhibit uniform structures; whereas, in the electrostatic correlation dominant regime, the polyelectrolyte solutions exhibit heterogeneous structures, in which the likely charged chains microscopically aggregate through counterion condensations. Increasing the intrinsic chain rigidity causes polyelectrolyte extension and hence moderately weakens the inter-chain clustering. Our work highlights the various, unique structures and molecular architectures of polyelectrolytes in solutions caused by the multi-body correlations between polyelectrolytes, counterions and asymmetric dipolar solvent molecules, which provides insights into the fundamental understanding of ion-containing polymers.

Towards understanding specific ion effects in aqueous media using thermodiffusion

Specific ion effects play an important role in scientific and technological processes. According to Hofmeister, the influence on the hydrogen bond network depends on the ion and leads to a specific order of the ions. Also thermodiffusion the mass transport caused by a temperature gradient is very sensitive to changes of the hydrogen bond network leading to a ranking according to hydrophilicity of the salt. Hence, we investigate various salt solutions in order to compare with the Hofmeister concept. We have studied three different sodium salts in water as a function of temperature (25-45[Formula: see text]C) and concentration (0.5-5 mol kg[Formula: see text]) using Thermal Diffusion Forced Rayleigh Scattering (TDFRS). The three anions studied, carbonate, acetate and thiocyanate, span the entire range of the Hofmeister series from hydrophilic to hydrophobic. We compare the results with the recent measurements of the corresponding potassium salts to see to what extent the cation changes the thermodiffusion of the salt.

Comparative experimental and computational study of synthetic and natural bottlebrush polyelectrolyte solutions

We systematically investigate model synthetic and natural bottlebrush polyelectrolyte solutions through an array of experimental techniques (osmometry and neutron and dynamic light scattering) along with molecular dynamics simulations to characterize and contrast their structures over a wide range of spatial and time scales. In particular, we perform measurements on solutions of aggrecan and the synthetic bottlebrush polymer, poly(sodium acrylate), and simulations of solutions of highly coarse-grained charged bottlebrush molecules having different degrees of side-branch density and inclusion of an explicit solvent and ion hydration effects. While both systems exhibit a general tendency toward supramolecular organization in solution, bottlebrush poly(sodium acrylate) solutions exhibit a distinctive "polyelectrolyte peak" in their structure factor, but no such peak is observed in aggrecan solutions. This qualitative difference in scattering properties, and thus polyelectrolyte solution organization, is attributed to a concerted effect of the bottlebrush polymer topology and the solvation of the polymer backbone and counterions. The coupling of the polyelectrolyte topological structure with the counterion distribution about the charged polymer molecules along with direct polymer segmental hydration makes their solution organization and properties "tunable," a phenomenon that has significant ramifications for biological function and disease as well as for numerous materials applications.

The effects of protein charge patterning on complex coacervation

The complex coacervation of proteins with other macromolecules has applications in protein encapsulation and delivery and for determining the function of cellular coacervates. Theoretical or empirical predictions for protein coacervates would enable the design of these coacervates with tunable and predictable structure-function relationships; unfortunately, no such theories exist. To help establish predictive models, the impact of protein-specific parameters on complex coacervation were probed in this study. The complex coacervation of sequence-specific, polypeptide-tagged, GFP variants and a strong synthetic polyelectrolyte was used to evaluate the effects of protein charge patterning on phase behavior. Phase portraits for the protein coacervates demonstrated that charge patterning dictates the protein's binodal phase boundary. Protein concentrations over 100 mg mL-1 were achieved in the coacervate phase, with concentrations dependent on the tag polypeptide sequence covalently attached to the globular protein domain. In addition to shifting the binodal phase boundary, polypeptide charge patterning provided entropic advantages over isotropically patterned proteins. Together, these results show that modest changes of only a few amino acids in the tag polypeptide sequence alter the coacervation thermodynamics and can be used to tune the phase behavior of polypeptides or proteins of interest.

Observation of General Entropy-Enthalpy Compensation Effect in the Relaxation of Wrinkled Polymer Nanocomposite Films

Surface-textured polymer nanocomposite (PNC) films are utilized in many device applications, and therefore understanding the relaxation behavior of such films is important. By extending an in situ wrinkle relaxation method, we observed that the thermal stability of wrinkled PNC films, both above and below the glass transition temperature (T_g), is proportional to a film's nanoparticle (polymer grafted and bare) concentration, with a slope that changes sign at a compensation temperature (T_comp) that is determined to be in the vicinity of the film's T_g. This provides unambiguous confirmation of entropy-enthalpy compensation (EEC) as a general feature of PNC films, implying that the stability of PNC films changes from being enhanced to becoming diminished by simply passing through this characteristic temperature, a phenomenon having evident practical ramifications. We suggest EEC will also arise in films where residual stresses are associated with the film fabrication process, which is relevant to nanotech device applications.

Unraveling the NaCl Concentration Effect on the First Stages of α-Synuclein Aggregation

Intraneuronal aggregation of the intrinsically disordered protein α-synuclein is at the core of Parkinson's disease and related neurodegenerative disorders. Several reports show that the concentration of salts in the medium heavily affects its aggregation rate and fibril morphology, but a characterization of the individual monomeric conformations underlying these effects is still lacking. In this work, we have applied our α-synuclein-optimized coarse-grained molecular dynamics approach to decipher the structural features of the protein monomer under a range of NaCl concentrations (0.0-1.0 M). The results show that key intramolecular contacts between the terminal domains are lost at intermediate concentrations (leading to extended conformations likely to fibrillate), but recovered at high concentrations (leading to compact conformations likely to evolve toward amorphous aggregates). The pattern of direct interactions of the terminal α-synuclein domains with Na⁺ and Cl^- ions plays a key role in explaining this effect. Our results are consistent with a recent study reporting a fibrillation enhancement at moderate NaCl concentrations but an inhibition at higher concentrations. The present work will contribute to improving our understanding of the structural features of monomeric α-synuclein, determining its NaCl-induced fibrillation propensity and the molecular basis of synucleinopathies, necessary for the future development of disease-halting therapies.

Influence of Kosmotrope and Chaotrope Salts on Water Structural Relaxation

The structural relaxation in water solutions of kosmotrope (structure maker) and chaotrope (structure breaker) salts, namely sodium chloride, potassium chloride, and cesium chloride, were studied through quasielastic neutron scattering measurements. We found that the collective dynamics relaxation time at the structure factor peak obtained using heavy water solutions shows a distinctively different behavior in the kosmotrope as opposed to the chaotrope solutions, increasing with the salt concentration in the former and decreasing in the latter. In both cases the trends are proportional to the concentration dependence of the relative viscosity of the solutions. These results indicate that kosmotropes and chaotropes influence the solution's viscosity by impacting in opposite ways the hydrogen bond network of water, strengthening it in one case and softening it in the other.

Disappearance of the polyelectrolyte peak in salt-free solutions

We investigate the nature of the polyelectrolyte peak in salt-free solutions by molecular dynamics simulations using a minimal model of polyelectrolyte solutions that includes an explicit solvent and counterions and small angle scattering experiments. It is found that the polyelectrolyte peak progressively disappears as the strength of solvation for the charged species is increased and the scattering profiles start to resemble those of neutral polymer solutions. The disappearance of the polyelectrolyte peak coincides with the emergence of attractive interchain interactions over a wide range of length scales. These findings provide insights into the microscopic origin of the polyelectrolyte peak.

Jones R, Lou J, Xia Y. Systematic investigation of synthetic polyelectrolyte bottlebrush solutions by neutron and dynamic light scattering, osmometry, and molecular dynamics simulation

There is a great interest in the synthesis and characterization of polyelectrolytes that mimic naturally occurring bottlebrush polyelectrolytes to capitalize on the unique properties of this class of macromolecules. Charged bottlebrush polymers form the protective mucus layer in the lungs, stomach, and orifices of animals and provide osmotic stabilization and lubrication to joints. In the present work, we systematically investigate bottlebrush poly(sodium acrylates) through a combination of measurements of solution properties (osmometry, small-angle neutron scattering, and dynamic light scattering) and molecular dynamics simulations, where the bottlebrush properties are compared in each case to their linear polymer counterparts. These complementary experimental and computational methods probe vastly different length- and timescales, allowing for a comprehensive characterization of the supermolecular structure and dynamics of synthetic polyelectrolyte bottlebrush molecules in solution.

Electric Double Layers with Surface Charge Regulation Using Density Functional Theory

Surprisingly, the local structure of electrolyte solutions in electric double layers is primarily determined by the solvent. This is initially unexpected as the solvent is usually a neutral species and not a subject to dominant Coulombic interactions. Part of the solvent dominance in determining the local structure is simply due to the much larger number of solvent molecules in a typical electrolyte solution.The dominant local packing of solvent then creates a space left for the charged species. Our classical density functional theory work demonstrates that the solvent structural effect strongly couples to the surface chemistry, which governs the charge and potential. In this article we address some outstanding questions relating double layer modeling. Firstly, we address the role of ion-ion correlations that go beyond mean field correlations. Secondly we consider the effects of a density dependent dielectric constant which is crucial in the description of a electrolyte-vapor interface.

Beyond Charge Balance: Counter-Cations in Polyoxometalate Chemistry

Polyoxometalates (POMs) are molecular metal-oxide anions applied in energy conversion and storage, manipulation of biomolecules, catalysis, as well as materials design and assembly. Although often overlooked, the interplay of intrinsically anionic POMs with organic and inorganic cations is crucial to control POM self-assembly, stabilization, solubility, and function. Beyond simple alkali metals and ammonium, chemically diverse cations including dendrimers, polyvalent metals, metal complexes, amphiphiles, and alkaloids allow tailoring properties for known applications, and those yet to be discovered. This review provides an overview of fundamental POM-cation interactions in solution, the resulting solid-state compounds, and behavior and properties that emerge from these POM-cation interactions. We will explore how application-inspired research has exploited cation-controlled design to discover new POM materials, which in turn has led to the quest for fundamental understanding of POM-cation interactions.

Influence of Sodium Salts on the Swelling and Rheology of Hydrophobically Cross-linked Hydrogels Determined by QCM-D

Hydrophobically modified copolymers provide a versatile platform of hydrogel materials for diverse applications, but the influence of salts on the swelling and material properties of this class of hydrogels has not been extensively studied. Here, we investigate model hydrogels with three different sodium salts with anions chosen from the classic Hofmeister series to determine how these counterions influence the swelling and mechanical properties of neutral hydrogels. The gel chosen was based on a statistical copolymer of dimethylacrylamide and 2-(N-ethylperfluorooctane sulfonamido) ethyl acrylate (FOSA). Our measurements utilize a quartz crystal microbalance with dissipation (QCM-D) to quantify both swelling and rheological properties of these gels. We find that a 1 mol/L solution of Na₂SO₄, corresponding to a kosmotropic anion, leads to nearly a 2.6-fold gel deswelling and correspondingly, the complex modulus increases by an order of magnitude under these solution conditions. In contrast, an initial increase in swelling and then a swelling maximum is observed for a 0.02 mol/L concentration in the case of a chaotropic anion, NaClO₄, but the changes in the degree of gel swelling in this system are not directly correlated with changes in the gel shear modulus. The addition of NaBr, an anion salt closer to the middle of the chaotropic to kosmotropic range, leads to hydrogel deswelling where the degree of deswelling and the shear modulus are both nearly independent of salt concentration. Overall, the observed trends are broadly consistent with more kosmotropic ions causing diminished solubility ("salting out") and strongly chaotropic ions causing improved solubility ("salting in"), a trend characteristic of the Hoffmeister series governing the solubility of many proteins and synthetic water-soluble polymers, but trends in the shear stiffness with gel swelling are clearly different from those normally observed in chemically cross-linked gels and are correspondingly difficult to interpret. The salt specificity of swelling and mechanical properties of nonionic hydrogels is important for any potential application in which a wide range of salt concentrations and types are encountered.

Comparative Study of Anomalous Size Dependence of Charged and Neutral Solute Diffusion in Water

We present a comparative study of size dependence of diffusion for charged and neutral solutes in water. Although both show nonmonotonicity of the size dependence of diffusion, their nature and origin are quite different. For neutral solutes, the peak position and the value of diffusion at the maximum are both independent of the solute-water interaction. Interestingly, for charged solutes, with an increase in solute-water interaction strength, the peak position shifts to lower solute sizes and with an increase in charge, it shifts to higher solute sizes. The diffusion value at the peak reduces with an increase in both solute-water interaction and solute charge. We show that all these features observed for charged solutes can be understood in terms of the interplay between ionic and nonionic interactions which is definitely absent for neutral solutes. Some of the earlier studies addressing the nonmonotonicity in diffusion did suggest the interplay between the two interactions to be the cause. However, this is the first time we show that such an interplay gives rise to the nonmonotonicity in the potential energy which is a prerequisite for obtaining the nonmonotonicity in the diffusion. Such nonmonotonicity in the potential energy is absent for neutral solutes.

Identifying Water-Anion Correlated Motion in Aqueous Solutions through Van Hove Functions

Electrolyte solutions are ubiquitous in materials in daily use and in biological systems. However, the understanding of their molecular and ionic dynamics, particularly those of their correlated motions, are elusive despite extensive experimental, theoretical, and numerical studies. Here we report the real-space observations of the molecular/ionic-correlated dynamics of aqueous salt (NaCl, NaBr, and NaI) solutions using the Van Hove functions obtained by high-resolution inelastic X-ray scattering measurement and molecular dynamics simulation. Our results directly depict the distance-dependent dynamics of aqueous salt solutions on the picosecond time scale and identify the changes in the anion-water correlations. This study demonstrates the capability of the real-space Van Hove function analysis to describe the local correlated dynamics in aqueous salt solutions.

Hydration of Hofmeister ions

Water dissolves salt into ions and then hydrates the ions to form an aqueous solution. Hydration of ions deforms the hydrogen bonding network and triggers the solution with what the pure water never shows such as conductivity, molecular diffusivity, thermal stability, surface stress, solubility, and viscosity, having enormous impact to many branches in biochemistry, chemistry, physics, and energy and environmental industry sectors. However, regulations for the solute-solute-solvent interactions are still open for exploration. From the perspective of the screened ionic polarization and O:H-O bond relaxation, this treatise features the recent progress and a perspective in understanding the hydration dynamics of Hofmeister ions in the typical YI, NaX, ZX₂, and NaT salt solutions (Y = Li, Na, K, Rb, Cs; X = F, Cl, Br, I; Z = Mg, Ca, Ba, Sr; T = ClO₄, NO₃, HSO₄, SCN). Phonon spectrometric analysis turned out the f(C) number fraction of bonds transition from the mode of deionized water to the hydrating. The linear f(C) ∝ C form features the invariant hydration volume of small cations that are fully-screened by their hydration H₂O dipoles. The nonlinear f(C) ∝ 1 - exp.(-C/C₀) form describes that the number insufficiency of the ordered hydrating H₂O dipoles partially screens the anions. Molecular anions show stronger yet shorter electric field of dipoles. The screened ionic polarization, inter-solute interaction, and O:H-O bond transition unify the solution conductivity, surface stress, viscosity, and critical energies for phase transition.

Polyelectrolyte association and solvation

There has been significant interest in the tendency of highly charged particles having the same charge to form dynamic clusters in solution, but an accepted theoretical framework that can account for this ubiquitous phenomenon has been slow to develop. The theoretical difficulties are especially great for flexible polyelectrolytes due to the additional complex coupling between the polyelectrolyte chain configurations and the spatial distribution of the ionic species in solution. For highly charged polyelectrolytes, this leads to the formation of a diffuse "polarizable" cloud of counter-ions around these polymers, an effect having significant implications for the function of proteins and other natural occurring polyelectrolytes, as emphasized long ago by Kirkwood and co-workers. To investigate this phenomenon, we perform molecular dynamics simulations of a minimal model of polyelectrolyte solutions that includes an explicit solvent and counter-ions, where the relative affinity of the counter-ions and the polymer for the solvent is tunable through the variation of the relative strength of the dispersion interactions of the polymer and ions. In particular, we find that these dispersion interactions can greatly influence the nature of the association between the polyelectrolyte chains under salt-free conditions. We calculate static and dynamic correlation functions to quantify the equilibrium structure and dynamics of these complex liquids. Based on our coarse-grained model of polyelectrolyte solutions, we identify conditions in which three distinct types of polyelectrolyte association arise. We rationalize these types of polyelectrolyte association based on the impact of the selective solvent affinity on the charge distribution and polymer solvation in these solutions. Our findings demonstrate the essential role of the solvent in the description of the polyelectrolyte solutions, as well as providing a guideline for the development of a more predictive theory of the properties of the thermodynamic and transport properties of these complex fluids.

Complex Coacervation in Polyelectrolytes from a Coarse-Grained Model

Complex coacervation refers to the formation of distinct liquid phases that arise when polyelectrolytes are mixed under appropriate polymer and salt concentrations. Molecular-level studies of coacervation have been limited. In this work, a coarse-grained model of the polymers and the corresponding counterions is proposed and used to simulate coacervation as a function of polymer length and overall salt concentration. Several sampling methods are used to determine the phase behavior of the underlying polymers. In particular, the results of simulations in different ensembles are shown to be consistent and to reproduce a number of phenomena observed in experiments, including the disruption of complexation by increasing ionic strength or by decreasing molecular weight. The coacervate concentrations determined from phase behavior calculations are then used to examine the rheology of the corresponding materials. By relying on long dynamic simulations, we are able to generate the dynamic response of the material in the form of dynamic moduli as a function of frequency, which are also found to compare favorably with experimental measurements.

Structure and rheology of polyelectrolyte complex coacervates

Scattering investigations of the structure and chain conformations, and the rheological properties of polyelectrolyte complexes (PECs) comprising model polyelectrolytes are presented. The use of charged polypeptides - (poly)-lysine and (poly)-glutamic acid with identical backbones allowed for facile tuning of the system parameters, including chain length, side-chain functionality, and chirality. Systematic studies using small-angle X-ray scattering (SAXS) of liquid PEC coacervates revealed a physical description of these materials as strongly screened semidilute polyelectrolyte solutions comprising oppositely charged chains. At the same time, solid PECs were found to be composed of hydrogen-bonding driven stiff ladder-like structures. While the coacervates behaved akin to semidilute polyelectrolyte solutions upon addition of salt, the solids were largely unaffected by it. Rheology measurements of PEC coacervates revealed a terminal relaxation regime, with an unusual plateauing of the storage modulus at low oscillation frequencies. The plateau may be ascribed to a combination of instrumental limitations and the long-range electrostatic interactions contributing to weak energy storage modes. Excellent superposition of the dynamic moduli was achieved by a time-salt superposition. The shift factors, however, varied more strongly than previously reported with added salt concentration.

Ion Correlation Effects in Salt-Doped Block Copolymers

We apply classical density functional theory to study how salt changes the microphase morphology of diblock copolymers. Polymers are freely jointed and one monomer type favorably interacts with ions, to account for the selective solvation that arises from different dielectric constants of the microphases. By including correlations from liquid state theory of an unbound reference fluid, the theory can treat chain behavior, microphase separation, ion correlations, and preferential solvation, at the same coarse-grained level. We show good agreement with molecular dynamics simulations.

The Influence of Polymer and Ion Solvation on the Conformational Properties of Flexible Polyelectrolytes

The study of the coupling between the conformational properties of a polyelectrolyte chain and the distribution of counter-ions surrounding the chain is important in developing predictive theories for more complex polymer materials, such as polyelectrolyte gels. We investigated the influence of solvent affinity to counter-ions and the polyelectrolyte backbone on the conformational properties of highly charged flexible polymer chains using molecular dynamics simulations that include both ions and an explicit solvent. We find that the solvation of the polyelectrolyte backbone can be achieved by either increasing the solvent affinity for the polyelectrolyte segments or by increasing the solvent affinity for the counter-ions. However, these two mechanisms influence the conformational properties of the polyelectrolyte chain in rather different ways, suggesting the inadequacy of polyelectrolyte solution models that treat the solvent as a continuum medium.

Communication: Counter-ion solvation and anomalous low-angle scattering in salt-free polyelectrolyte solutions

We investigate the influence of counter-ion solvation on the homogeneity of salt-free polyelectrolyte solutions based on a coarse-grained model that includes an explicit solvent. We show that the solvation of the counter-ions can cause a transformation between a nearly homogeneous to a non-uniform polymer solution, in which there is both a chain clustering and the formation of large charge-free domains, i.e., "voids." The emergence of these heterogeneous structures induced by counter-ion solvation is accompanied by the localization and formation of counter-ion rich domains that are symptomatic of emergent effective long-range attractive interchain interactions.