Effect of counterion size on polyelectrolyte conformations and thermodynamics.

We present a theoretical model to study the effect of counterion size on the effective charge, size, and thermodynamic behavior of a single, isolated, and flexible polyelectrolyte (PE) chain. We analyze how altering counterion size modifies the energy and entropy contributions to the system, including the ion-pair free energy, excluded volume interactions, entropy of free and condensed ions, and dipolar attraction among monomer-counterion pairs, which result in competing effects challenging intuitive predictions. The PE self-energy is calculated using the Edwards-Muthukumar Hamiltonian, considering a Gaussian monomer distribution for the PE. The condensed ions are assumed to be confined within a cylindrical volume around the PE backbone. The dipolar and excluded volume interactions are described by the second and third virial coefficients. The assumption of freely rotating dipoles results in a first-order coil-globule transition of the PE chain. A more realistic, weaker dipolar attraction, parameterized in our theory, shifts it to a second-order continuous transition. We calculate the size scaling-exponent of the PE and find exponents according to the relative dominance of the electrostatic, excluded volume, or dipolar effects. We further identify the entropy- and energy-driven regimes of the effective charge and conformation of the PE, highlighting the interplay of free ion entropy and ion-pair energy with varying electrostatic strengths. The crossover strength, dependent on the counterion size, indicates that diminishing sizes favor counterion condensation at the expense of free ion entropy. The predictions of the model are consistent with trends in simulations and generalize the findings of the point-like counterion theories.

Driving forces of the complex formation between highly charged disordered proteins.

Highly disordered complexes between oppositely charged intrinsically disordered proteins present a new paradigm of biomolecular interactions. Here, we investigate the driving forces of such interactions for the example of the highly positively charged linker histone H1 and its highly negatively charged chaperone, prothymosin α (ProTα). Temperature-dependent single-molecule Förster resonance energy transfer (FRET) experiments and isothermal titration calorimetry reveal ProTα-H1 binding to be enthalpically unfavorable, and salt-dependent affinity measurements suggest counterion release entropy to be an important thermodynamic driving force. Using single-molecule FRET, we also identify ternary complexes between ProTα and H1 in addition to the heterodimer at equilibrium and show how they contribute to the thermodynamics observed in ensemble experiments. Finally, we explain the observed thermodynamics quantitatively with a mean-field polyelectrolyte theory that treats counterion release explicitly. ProTα-H1 complex formation resembles the interactions between synthetic polyelectrolytes, and the underlying principles are likely to be of broad relevance for interactions between charged biomolecules in general.

Polymer complexation: Partially ionizable asymmetric polyelectrolytes.

Theories of bulk coacervation of oppositely charged polyelectrolytes (PE) obscure single molecule level thermodynamic details, considered significant for coacervate equilibrium, whereas simulations account for only pairwise Coulomb interaction. Also, studies of effects of asymmetry on PE complexation are rare compared to symmetric PEs. We develop a theoretical model, accounting for all entropic and enthalpic contributions at the molecular level, and the mutual segmental screened Coulomb and excluded volume interactions between two asymmetric PEs, by constructing a Hamiltonian following Edwards and Muthukumar. Assuming maximal ion-pairing in the complex, the system free energy comprising configurational entropy of the polyions and free-ion entropy of the small ions is minimized. The effective charge and size of the complex, larger than sub-Gaussian globules as for symmetric chains, increase with asymmetry in polyion length and charge density. The thermodynamic drive for complexation is found to increase with ionizability of symmetric polyions and with a decrease in asymmetry in length for equally ionizable polyions. The crossover Coulomb strength demarcating the ion-pair enthalpy-driven (low strength) and counterion release entropy-driven (high strength) is marginally dependent on the charge density, because so is the degree of counterion condensation, and strongly dependent on the dielectric environment and salt. The key results match the trends in simulations. The framework may provide a direct way to calculate thermodynamic dependencies of complexation on experimental parameters such as electrostatic strength and salt, thus to better analyze and predict observed phenomena for different sets of polymer pairs.

Polyelectrolyte complexation of two oppositely charged symmetric polymers: A minimal theory.

Interplay of Coulomb interaction energy, free ion entropy, and conformational elasticity is a fascinating aspect in polyelectrolytes (PEs). We develop a theory for complexation of two oppositely charged PEs, a process known to be the precursor to the formation of complex coacervates in PE solutions, to explore the underlying thermodynamics of complex formation, at low salts. The theory considers general degrees of solvent polarity and dielectricity within an implicit solvent model, incorporating a varying Coulomb strength. Explicit calculation of the free energy of complexation and its components indicates that the entropy of free counterions and salt ions and the Coulomb enthalpy of bound ion-pairs dictate the equilibrium of PE complexation. This helps decouple the self-consistent dependency of charge and size of the uncomplexed parts of the polyions, derive an analytical expression for charge, and evaluate the free energy components as functions of chain overlap. Complexation is observed to be driven by enthalpy gain at low Coulomb strengths, driven by entropy gain of released counterions but opposed by enthalpy loss due to reduction of ion-pairs at moderate Coulomb strengths, and progressively less favorable due to enthalpy loss at even higher Coulomb strengths. The total free energy of the system is found to decrease linearly with an overlap of chains. Thermodynamic predictions from our model are in good quantitative agreement with simulations in literature.

Understanding the Thermodynamics of the Binding of PAMAM Dendrimers to Graphene: A Combined Analytical and Simulation Study.

We investigate the thermodynamics of the binding of a poly(amidoamine) dendrimer to an uncharged graphene sheet as a function of the pH level using umbrella sampling simulations and a mean-field theory for generations three and four. We find that the dendrimer strongly binds to the graphene sheet ( O (100) kcal/mol) from our potential of mean force (PMF) calculations. In specific, we find that the dendrimer binds the most at neutral pH (∼7) and the least at low pH (∼4). We explain this nonmonotonic nature of the dendrimer's adsorption by studying the interactions contributing to the PMF, i.e., the dendrimer-graphene, dendrimer-water, and dendrimer-ion interactions. We also corroborate our PMF calculations with molecular mechanics generalized Born surface area analysis and free energies obtained from a mean-field theory of Flory-Huggins-Debye-Hückel type [ Muthukumar , M. , J. Chem. Phys. 2010 , 132 , 084901 ], including electrostatic interactions. We find that the van der Waals interactions between the dendrimer and the graphene alone cannot capture the accurate trends in the binding free energies (BEs) as a function of pH. The solvent and the counterions present in the system are also found to have a major influence on these trends. We demonstrate that the dendrimer-graphene and dendrimer-water interactions become favorable, whereas the dendrimer-ion interaction becomes unfavorable, as the dendrimer binds to graphene. These opposing effects lead to the observed nonmonotonicity in the BE trends. Our theoretical model also reproduces these trends in the subinteractions contributing to the PMF. To the best of our knowledge, this is a novel attempt where an equivalence between theory and simulations is made in the context of the dendrimer's adsorption.

Charge-regularized swelling kinetics of polyelectrolyte gels: Elasticity and diffusion.

We apply a recently developed method [S. Sen and A. Kundagrami, J. Chem. Phys. 143, 224904 (2015)], using a phenomenological expression of osmotic stress, as a function of polymer and charge densities, hydrophobicity, and network elasticity for the swelling of spherical polyelectrolyte (PE) gels with fixed and variable charges in a salt-free solvent. This expression of stress is used in the equation of motion of swelling kinetics of spherical PE gels to numerically calculate the spatial profiles for the polymer and free ion densities at different time steps and the time evolution of the size of the gel. We compare the profiles of the same variables obtained from the classical linear theory of elasticity and quantitatively estimate the bulk modulus of the PE gel. Further, we obtain an analytical expression of the elastic modulus from the linearized expression of stress (in the small deformation limit). We find that the estimated bulk modulus of the PE gel decreases with the increase of its effective charge for a fixed degree of deformation during swelling. Finally, we match the gel-front locations with the experimental data, taken from the measurements of charged reversible addition-fragmentation chain transfer gels to show an increase in gel-size with charge and also match the same for PNIPAM (uncharged) and imidazolium-based (charged) minigels, which specifically confirms the decrease of the gel modulus value with the increase of the charge. The agreement between experimental and theoretical results confirms general diffusive behaviour for swelling of PE gels with a decreasing bulk modulus with increasing degree of ionization (charge). The new formalism captures large deformations as well with a significant variation of charge content of the gel. It is found that PE gels with large deformation but same initial size swell faster with a higher charge.

Kinetics of swelling of polyelectrolyte gels: Fixed degree of ionization.

The swelling kinetics of uncharged and charged polymer (polyelectrolyte) gels in salt-free conditions is studied in one dimension by solving the constitutive equation of motion (Newton's law for the elementary gel volume) of the displacement variable by two theoretical methods: one in which the classical definition of stress is used with the bulk modulus taken as a parameter, and the other in which a phenomenological expression of the osmotic stress as a function of polymer density and degree of ionization is taken as an input to the dynamics. The time-evolution profiles for spatially varying polymer density and stress, along with the location of the gel-solvent interface, are obtained from the two methods. We show that both the polymer density (volume fraction) and stress inside the gel follow expected behaviours of being maximum for the uniformly shrunken gel, and relaxing slowly to the lowest values as the gel approaches equilibrium. We further show that, by comparing the temporal profiles of the gel-solvent interface and other variables between the two methods, one may attempt to assign an effective bulk modulus to the polyelectrolyte gel as a function of the degree of ionization and other parameters of the gel such as hydrophobicity, cross-link density, and the temperature. The major result we get is that the effective bulk modulus of a polyelectrolyte gel increases monotonically with its degree of ionization. In the process of identifying the parameters for a monotonic swelling, we calculated using a well-known expression of the free energy the equilibrium results of two-phase co-existence and the critical point of a polyelectrolyte gel with a fixed degree of ionization.

Effective charge and coil-globule transition of a polyelectrolyte chain.

Considering the adsorption of counterions on an isolated polyelectrolyte (PE) chain and using a variational theory, phase boundaries and the critical point for the first-order coil-globule transition are calculated. The transition is induced cooperatively by counterion adsorption and chain conformations and the calculation is done self-consistently. The size of the PE chain is a single-valued function of charge. The discontinuous transition of the coil size is accompanied by a discontinuous transition of the charge. Phase boundaries for the coil-globule transitions induced by both Coulomb strength (inverse temperature or dielectric constant) and ionic strength (salt) show that the PE chain collapses at a substantially lower Coulomb strength in the presence of salt. In the expanded state of the coil, an analytical formula is derived for the effective charge of the chain for conditions where the coupling between chain conformations and counterion adsorption is weak. In general, the dielectric heterogeneity of the solvent close to the polymer backbone is found to play a crucial role in the charge regularization and the chain collapse.

Charge regularization in phase separating polyelectrolyte solutions.

Theoretical investigations of phase separation in polyelectrolyte solutions have so far assumed that the effective charge of the polyelectrolyte chains is fixed. The ability of the polyelectrolyte chains to self-regulate their effective charge due to the self-consistent coupling between ionization equilibrium and polymer conformations, depending on the dielectric constant, temperature, and polymer concentration, affects the critical phenomena and phase transitions drastically. By considering salt-free polyelectrolyte solutions, we show that the daughter phases have different polymer charges from that of the mother phase. The critical point is also altered significantly by the charge self-regularization of the polymer chains. This work extends the progress made so far in the theory of phase separation of strong polyelectrolyte solutions to a higher level of understanding by considering chains which can self-regulate their charge.

Counterion adsorption on flexible polyelectrolytes: comparison of theories.

Counterion adsorption on a flexible polyelectrolyte chain in a spherical cavity is considered by taking a "permuted" charge distribution on the chain so that the "adsorbed" counterions are allowed to move along the backbone. We compute the degree of ionization by using self-consistent field theory (SCFT) and compare with the previously developed variational theory. Analysis of various contributions to the free energy in both theories reveals that the equilibrium degree of ionization is attained mainly as an interplay of the adsorption energy of counterions on the backbone, the translational entropy of the small ions, and their correlated density fluctuations. Degree of ionization computed from SCFT is significantly lower than that from the variational formalism. The difference is entirely due to the density fluctuations of the small ions in the system, which are accounted for in the variational procedure. When these fluctuations are deliberately suppressed in the truncated variational procedure, there emerges a remarkable quantitative agreement in the various contributing factors to the equilibrium degree of ionization, in spite of the fundamental differences in the approximations and computational procedures used in these two schemes. Furthermore, it is found that the total free energies from the truncated variational procedure and the SCFT are in quantitative agreement at low monomer densities and differ from each other at higher monomer densities. The disagreement at higher monomer densities is due to the inability of the variational calculation to accurately compute the solvent entropy at higher concentrations. A comparison of electrostatic energies (which are relatively small) reveals that the Debye-Hückel estimate used in the variational theory is an overestimation of electrostatic energy as compared to the Poisson-Boltzmann estimate. Nevertheless, since the significant effects from density fluctuations of small ions are not captured by the SCFT, and due to the close agreement between SCFT and the other contributing factors in the more transparent variational procedure, the latter is a better computational tool for obtaining the degree of ionization.

Theory of competitive counterion adsorption on flexible polyelectrolytes: divalent salts.

The counterion distribution around an isolated flexible polyelectrolyte in the presence of a divalent salt is evaluated using the adsorption model [M. Muthukumar, J. Chem. Phys. 120, 9343 (2004)] that considers the Bjerrum length, salt concentration, and local dielectric heterogeneity as physical variables in the system. Self-consistent calculations of effective charge and size of the polymer show that divalent counterions replace condensed monovalent counterions in competitive adsorption. The theory further predicts that at modest physical conditions for a flexible polyelectrolytes such as sodium polystyrene sulfonate in aqueous solutions polymer charge is compensated and reversed with increasing divalent salt. Consequently, the polyelectrolyte shrinks and reswells. Lower temperatures and higher degrees of dielectric heterogeneity between chain backbone and solvent enhance condensation of all species of ions. Complete diagrams of states for the effective charge calculated as functions of the Coulomb strength and salt concentration suggest that (a) overcharging requires a minimum Coulomb strength and (b) progressively higher presence of salt recharges the polymer due to either electrostatic screening (for low Coulomb strengths) or coion condensation (for high Coulomb strengths). Consideration of ion-bridging by divalent counterions leads to a first-order collapse of polyelectrolytes in modest presence of divalent salts and at higher Coulomb strengths. The authors' theoretical predictions are in agreement with the generic results from experiments and simulations.

Continuum theory of polymer crystallization.

We present a kinetic model of crystal growth of polymers of finite molecular weight. Experiments help to classify polymer crystallization broadly into two kinetic regimes. One is observed in melts or in high molar mass polymer solutions and is dominated by nucleation control with G approximately exp(1/TDeltaT), where G is the growth rate and DeltaT is the supercooling. The other is observed in low molar mass solutions (as well as for small molecules) and is diffusion controlled with G approximately DeltaT, for small DeltaT. Our model unifies these two regimes in a single formalism. The model accounts for the accumulation of polymer chains near the growth front and invokes an entropic barrier theory to recover both limits of nucleation and diffusion control. The basic theory applies to both melts and solutions, and we numerically calculate the growth details of a single crystal in a dilute solution. The effects of molecular weight and concentration are also determined considering conventional polymer dynamics. Our theory shows that entropic considerations, in addition to the traditional energetic arguments, can capture general trends of a vast range of phenomenology. Unifying ideas on crystallization from small molecules and from flexible polymer chains emerge from our theory.

Structure of twist-grain-boundary-C phases.

We study properties of the Renn-Lubensky twist-grain-boundary-C (TGB(C)) phase, with layer normal rotating in a plane perpendicular to the pitch axis, and the Bordeaux TGB(C) phase, with the layer normal rotating on a cone parallel to the pitch axis near the upper critical twist k(c2) marking the transition to the cholesteric phase. We introduce a generalized model free energy for the smectic-C phase that allows either TGB(C) phase to be stable, and we calculate k(c2) and the order-parameter profile, which shows only modest spatial variation, for both phases.