Osmotic swelling behavior of surface-charged ionic microgels.

In recent years, ionic microgels have garnered much attention due to their unique properties, especially their stimulus-sensitive swelling behavior. The tunable response of these soft, permeable, compressible, charged colloidal particles is increasingly attractive for applications in medicine and biotechnologies, such as controlled drug delivery, tissue engineering, and biosensing. The ability to model and predict variation of the osmotic pressure of a single microgel with respect to changes in particle properties and environmental conditions proves vital to such applications. In this work, we apply both nonlinear Poisson-Boltzmann theory and molecular dynamics simulation to ionic microgels (macroions) in the cell model to compute density profiles of microions (counterions, coions), single-microgel osmotic pressure, and equilibrium swelling ratios of spherical microgels whose fixed charge is confined to the macroion surface. The basis of our approach is an exact theorem that relates the electrostatic component of the osmotic pressure to the microion density profiles. Close agreement between theory and simulation serves as a consistency check to validate our approach. We predict that surface-charged microgels progressively deswell with increasing microgel concentration, starting well below close packing, and with increasing salt concentration, in qualitative agreement with experiments. Comparison with previous results for microgels with fixed charge uniformly distributed over their volume demonstrates that surface-charged microgels deswell more rapidly than volume-charged microgels. We conclude that swelling behavior of ionic microgels in solution is sensitive to the distribution of fixed charge within the polymer-network gel and strongly depends on bulk concentrations of both microgels and salt ions.

Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods.

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion-microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations.

Structural correlations in highly asymmetric binary charged colloidal mixtures.

We explore structural correlations of strongly asymmetric mixtures of binary charged colloids within the primitive model of electrolytes considering large charge and size ratios of 10 and higher. Using computer simulations with explicit microions, we obtain the partial pair correlation functions between the like-charged colloidal macroions. Interestingly the big-small correlation peak amplitude is smaller than that of the big-big and small-small macroion correlation peaks, which is unfamiliar for additive repulsive interactions. Extracting optimal effective microion-averaged pair interactions between the macroions, we find that on top of non-additive Yukawa-like repulsions an additional shifted Gaussian attractive potential between the small macroions is needed to accurately reproduce their correct pair correlations. For small Coulomb couplings, the behavior is reproduced in a coarse-grained theory with microion-averaged effective interactions between the macroions. However, the accuracy of the theory deteriorates with increasing Coulomb coupling. We emphasize the relevance of entropic interactions exerted by the microions on the macroions. Our results are experimentally verifiable in binary mixtures of micron-sized colloids and like-charge nanoparticles.

Beyond simple self-healing: How anisotropic nanogels adapt their shape to their environment.

The response of soft colloids to crowding depends sensitively on the particles' compressibility. Nanogel suspensions provide model systems that are often studied to better understand the properties of soft materials and complex fluids from the formation of colloidal crystals to the flow of viruses, blood, or platelet cells in the body. Large spherical nanogels, when embedded in a matrix of smaller nanogels, have the unique ability to spontaneously deswell to match their size to that of the nanogel composing the matrix. In contrast to hard colloids, this self-healing mechanism allows for crystal formation without giving rise to point defects or dislocations. Here, we show that anisotropic ellipsoidal nanogels adapt both their size and their shape depending on the nature of the particles composing the matrix in which they are embedded. Using small-angle neutron scattering with contrast variation, we show that ellipsoidal nanogels become spherical when embedded in a matrix of spherical nanogels. In contrast, the anisotropy of the ellipsoid is enhanced when they are embedded in a matrix of anisotropic nanogels. Our experimental data are supported by Monte Carlo simulations that reproduce the trend of decreasing aspect ratio of ellipsoidal nanogels with increasing crowding by a matrix of spherical nanogels.

Influence of solvent quality on depletion potentials in colloid-polymer mixtures.

As first explained by the classic Asakura-Oosawa (AO) model, effective attractive forces between colloidal particles induced by depletion of nonadsorbing polymers can drive demixing of colloid-polymer mixtures into colloid-rich and colloid-poor phases, with practical relevance for purification of water, stability of foods and pharmaceuticals, and macromolecular crowding in biological cells. By idealizing polymer coils as effective penetrable spheres, the AO model qualitatively captures the influence of polymer depletion on thermodynamic phase behavior of colloidal suspensions. In previous work, we extended the AO model to incorporate aspherical polymer conformations and showed that fluctuating shapes of random-walk coils can significantly modify depletion potentials [W. K. Lim and A. R. Denton, Soft Matter 12, 2247 (2016); J. Chem. Phys. 144, 024904 (2016)]. We further demonstrated that the shapes of polymers in crowded environments sensitively depend on solvent quality [W. J. Davis and A. R. Denton, J. Chem. Phys. 149, 124901 (2018)]. Here, we apply Monte Carlo simulation to analyze the influence of solvent quality on depletion potentials in mixtures of hard-sphere colloids and nonadsorbing polymer coils, modeled as ellipsoids whose principal radii fluctuate according to random-walk statistics. We consider both self-avoiding and non-self-avoiding random walks, corresponding to polymers in good and theta solvents, respectively. Our simulation results demonstrate that depletion of polymers of equal molecular weight induces much stronger attraction between colloids in good solvents than in theta solvents and confirm that depletion interactions are significantly influenced by aspherical polymer conformations.

Osmotic pressure and swelling behavior of ionic microcapsules.

Ionic microcapsules are hollow shells of hydrogel, typically 10-1000 nm in radius, composed of cross-linked polymer networks that become charged and swollen in a good solvent. The ability of microcapsules to swell/deswell in response to changes in external stimuli (e.g., temperature, pH, and ionic strength) suits them to applications, such as drug delivery, biosensing, and catalysis. The equilibrium swelling behavior of ionic microcapsules is determined by a balance of electrostatic and elastic forces. The electrostatic component of the osmotic pressure of a microcapsule-the difference in the pressure between the inside and outside of the particle-plays a vital role in determining the swelling behavior. Within the spherical cell model, we derive exact expressions for the radial pressure profile and for the electrostatic and gel components of the osmotic pressure of a microcapsule, which we compute via Poisson-Boltzmann theory and molecular dynamics simulation. For the gel component, we use the Flory-Rehner theory of polymer networks. By combining the electrostatic and gel components of the osmotic pressure, we compute the equilibrium size of ionic microcapsules as a function of particle concentration, shell thickness, and valence. We predict concentration-driven deswelling at relatively low concentrations at which steric interactions between particles are weak and demonstrate that this response can be attributed to crowding-induced redistribution of counterions. Our approach may help to guide the design and applications of smart stimuli-responsive colloidal particles.

Tailoring the Cavity of Hollow Polyelectrolyte Microgels.

The authors demonstrate how the size and structure of the cavity of hollow charged microgels may be controlled by varying pH and ionic strength. Hollow charged microgels based on N-isopropylacrylamide with ionizable co-monomers (itaconic acid) combine advanced structure with enhanced responsiveness to external stimuli. Structural advantages accrue from the increased surface area provided by the extra internal surface. Extreme sensitivity to pH and ionic strength due to ionizable moieties in the polymer network differentiates these soft colloidal particles from their uncharged counterparts, which sustain a hollow structure only at cross-link densities sufficiently high that stimuli sensitivity is reduced. Using small-angle neutron and light scattering, increased swelling of the network in the charged state accompanied by an expanded internal cavity is observed. Upon addition of salt, the external fuzziness of the microgel surface diminishes while the internal fuzziness grows. These structural changes are interpreted via Poisson-Boltzmann theory in the cell model.

Osmotic pressure of permeable ionic microgels: Poisson-Boltzmann theory and exact statistical mechanical relations in the cell model.

Ionic microgels are soft colloidal particles, composed of crosslinked polymer networks, which ionize and swell when dispersed in a good solvent. Swelling of these permeable, compressible particles involves a balance of electrostatic, elastic, and mixing contributions to the single-particle osmotic pressure. The electrostatic contribution depends on the distributions of mobile counterions and coions and of fixed charge on the polymers. Within the cell model, we employ two complementary methods to derive the electrostatic osmotic pressure of ionic microgels. In Poisson-Boltzmann (PB) theory, we minimize a free energy functional with respect to the electrostatic potential to obtain the bulk pressure. From the pressure tensor, we extract the electrostatic and gel contributions to the total pressure. In a statistical mechanical approach, we vary the free energy with respect to microgel size to obtain exact relations for the microgel electrostatic osmotic pressure. We present results for planar, cylindrical, and spherical geometries. For models of membranes and microgels with fixed charge uniformly distributed over their surface or volume, we derive analogs of the contact value theorem for charged colloids. We validate these relations by solving the PB equation and computing ion densities and osmotic pressures. When implemented within PB theory, the two methods yield identical electrostatic osmotic pressures for surface-charged microgels. For volume-charged microgels, the exact electrostatic osmotic pressure equals the average of the corresponding PB profile over the gel volume. We demonstrate that swelling of ionic microgels depends on the variation of the electrostatic pressure inside the particle and discuss implications for interpreting experiments.

Modeling deswelling, thermodynamics, structure, and dynamics in ionic microgel suspensions.

Ionic microgel particles in a good solvent swell to an equilibrium size determined by a balance of electrostatic and elastic forces. When crowded, ionic microgels deswell owing to a redistribution of microions inside and outside the particles. The concentration-dependent deswelling affects the interactions between the microgels and, consequently, the suspension properties. We present a comprehensive theoretical study of crowding effects on thermodynamic, structural, and dynamic properties of weakly cross-linked ionic microgels in a good solvent. The microgels are modeled as microion- and solvent-permeable colloidal spheres with fixed charge uniformly distributed over the polymer gel backbone, whose elastic and solvent-interaction free energies are described using the Flory-Rehner theory. Two mean-field methods for calculating the crowding-dependent microgel radius are investigated and combined with calculations of the net microgel charge characterizing the electrostatic part of an effective microgel pair potential, with charge renormalization accounted for. Using this effective pair potential, thermodynamic and static suspension properties are calculated, including the osmotic pressure and microgel pair distribution function. The latter is used in our calculations of dynamic suspension properties, where we account for hydrodynamic interactions. Results for diffusion and rheological properties are presented over ranges of microgel concentration and charge. We show that deswelling mildly enhances self- diffusion and collective diffusion and the osmotic pressure, lowers the suspension viscosity, and significantly shifts the suspension crystallization point to higher concentrations. This paper presents a bottom-up approach to efficiently computing suspension properties of crowded ionic microgels using single-particle characteristics.

Concentration-dependent swelling and structure of ionic microgels: simulation and theory of a coarse-grained model.

We study swelling and structural properties of ionic microgel suspensions within a comprehensive coarse-grained model that combines the polymeric and colloidal natures of microgels as permeable, compressible, charged spheres governed by effective interparticle interactions. The model synthesizes the Flory-Rehner theory of cross-linked polymer gels, the Hertz continuum theory of effective elastic interactions, and a theory of density-dependent effective electrostatic interactions. Implementing the model using Monte Carlo simulation and thermodynamic perturbation theory, we compute equilibrium particle size distributions, swelling ratios, volume fractions, net valences, radial distribution functions, and static structure factors as functions of concentration. Trial Monte Carlo moves comprising particle displacements and size variations are accepted or rejected based on the total change in elastic and electrostatic energies. The theory combines first-order thermodynamic perturbation and variational free energy approximations. For illustrative system parameters, theory and simulation agree closely at concentrations ranging from dilute to beyond particle overlap. With increasing concentration, as microgels deswell, we predict a decrease in the net valence and an unusual saturation of pair correlations. Comparison with experimental data for deionized, aqueous suspensions of PNIPAM particles demonstrates the capacity of the coarse-grained model to predict and interpret measured swelling behavior.

Influence of solvent quality on conformations of crowded polymers.

The structure and function of polymers in confined environments, e.g., biopolymers in the cytoplasm of a cell, are strongly affected by macromolecular crowding. To explore the influence of solvent quality on conformations of crowded polymers, we model polymers as penetrable ellipsoids, whose shape fluctuations are governed by the statistics of self-avoiding walks, appropriate for a polymer in a good solvent. Within this coarse-grained model, we perform Monte Carlo simulations of mixtures of polymers and hard-nanosphere crowders, including trial changes in polymer size and shape. Penetration of polymers by crowders is incorporated via a free energy cost predicted by polymer field theory. To analyze the impact of crowding on polymer conformations in different solvents, we compute the average polymer shape distributions, radius of gyration, volume, and asphericity over ranges of the polymer-to-crowder size ratio and crowder volume fraction. The simulation results are accurately predicted by a free-volume theory of polymer crowding. Comparison of results for polymers in good and theta solvents indicates that excluded-volume interactions between polymer segments significantly affect crowding, especially in the limit of crowders much smaller than polymers. Our approach may help to motivate future experimental studies of polymers in crowded environments, with possible relevance for drug delivery and gene therapy.

Structure and stability of charged colloid-nanoparticle mixtures.

Physical properties of colloidal materials can be modified by addition of nanoparticles. Within a model of like-charged mixtures of particles governed by effective electrostatic interactions, we explore the influence of charged nanoparticles on the structure and thermodynamic phase stability of charge-stabilized colloidal suspensions. Focusing on salt-free mixtures of particles of high size and charge asymmetry, interacting via repulsive Yukawa effective pair potentials, we perform molecular dynamics simulations and compute radial distribution functions and static structure factors. Analysis of these structural properties indicates that increasing the charge and concentration of nanoparticles progressively weakens correlations between charged colloids. We show that addition of charged nanoparticles to a suspension of like-charged colloids can induce a colloidal crystal to melt and can facilitate aggregation of a fluid suspension due to attractive van der Waals interactions. We attribute the destabilizing influence of charged nanoparticles to enhanced screening of electrostatic interactions, which weakens repulsion between charged colloids. This interpretation is consistent with recent predictions of an effective interaction theory of charged colloid-nanoparticle mixtures.

Influence of polymer shape on depletion potentials and crowding in colloid-polymer mixtures.

Depletion-induced interactions between colloids in colloid-polymer mixtures depend in range and strength on size, shape, and concentration of depletants. Crowding by colloids in turn affects shapes of polymer coils, such as biopolymers in biological cells. By simulating hard-sphere colloids and random-walk polymers, modeled as fluctuating ellipsoids, we compute depletion-induced potentials and polymer shape distributions. Comparing results with exact density-functional theory calculations, molecular simulations, and experiments, we show that polymer shape fluctuations play an important role in depletion and crowding phenomena.

Swelling, structure, and phase stability of compressible microgels.

Microgels are soft colloidal particles that, when dispersed in a solvent, swell and deswell in response to changes in environmental conditions, such as temperature, concentration, and pH. Using Monte Carlo simulation, we model bulk suspensions of microgels that interact via Hertzian elastic interparticle forces and can expand or contract via trial moves that allow particles to change size in accordance with the Flory-Rehner free energy of cross-linked polymer gels. We monitor the influence of particle compressibility, size fluctuations, and concentration on bulk structural and thermal properties by computing particle swelling ratios, radial distribution functions, static structure factors, osmotic pressures, and freezing densities. For microgels in the nanoscale size range, particle compressibility and associated size fluctuations suppress crystallization, shifting the freezing transition to a higher density than for the hard-sphere fluid. As densities increase beyond close packing, microgels progressively deswell, while their intrinsic size distribution grows increasingly polydisperse.

Counterion-induced swelling of ionic microgels.

Ionic microgel particles, when dispersed in a solvent, swell to equilibrium sizes that are governed by a balance between electrostatic and elastic forces. Tuning of particle size by varying external stimuli, such as pH, salt concentration, and temperature, has relevance for drug delivery, microfluidics, and filtration. To model swelling of ionic microgels, we derive a statistical mechanical theorem, which proves exact within the cell model, for the electrostatic contribution to the osmotic pressure inside a permeable colloidal macroion. Applying the theorem, we demonstrate how the distribution of counterions within an ionic microgel determines the internal osmotic pressure. By combining the electrostatic pressure, which we compute via both Poisson-Boltzmann theory and molecular dynamics simulation, with the elastic pressure, modeled via the Flory-Rehner theory of swollen polymer networks, we show how deswelling of ionic microgels with increasing concentration of particles can result from a redistribution of counterions that reduces electrostatic pressure. A linearized approximation for the electrostatic pressure, which proves remarkably accurate, provides physical insight and greatly eases numerical calculations for practical applications. Comparing with experiments, we explain why soft particles in deionized suspensions deswell upon increasing concentration and why this effect may be suppressed at higher ionic strength. The failure of the uniform ideal-gas approximation to adequately account for counterion-induced deswelling below close packing of microgels is attributed to neglect of spatial variation of the counterion density profile and the electrostatic pressure of incompletely neutralized macroions.

Depletion-induced forces and crowding in polymer-nanoparticle mixtures: Role of polymer shape fluctuations and penetrability.

Depletion forces and macromolecular crowding govern the structure and function of biopolymers in biological cells and the properties of polymer nanocomposite materials. To isolate and analyze the influence of polymer shape fluctuations and penetrability on depletion-induced interactions and crowding by nanoparticles, we model polymers as effective penetrable ellipsoids, whose shapes fluctuate according to the probability distributions of the eigenvalues of the gyration tensor of an ideal random walk. Within this model, we apply Monte Carlo simulation methods to compute the depletion-induced potential of mean force between hard nanospheres and crowding-induced shape distributions of polymers in the protein limit, in which polymer coils can be easily penetrated by smaller nanospheres. By comparing depletion potentials from simulations of ellipsoidal and spherical polymer models with predictions of polymer field theory and free-volume theory, we show that polymer depletion-induced interactions and crowding depend sensitively on polymer shapes and penetrability, with important implications for bulk thermodynamic phase behavior.

Ion density deviations in semipermeable ionic microcapsules.

By implementing the nonlinear Poisson-Boltzmann theory in a cell model, we theoretically investigate the influence of polyelectrolye gel permeability on ion densities and pH deviations inside the cavities of ionic microcapsules. Our calculations show that variations in permeability of a charged capsule shell cause a redistribution of ion densities within the capsule, which ultimately affects the pH deviation and Donnan potential induced by the electric field of the shell. We find that semipermeable capsules can induce larger pH deviations inside their cavities that can permeable capsules. Furthermore, with increasing capsule charge, the influence of permeability on pH deviations progressively increases. Our theory, while providing a self-consistent method for modeling the influence of permeability on fundamental properties of ionic microgels, makes predictions of practical significance for the design of microcapsules loaded with fluorescent dyes, which can serve as biosensors for diagnostic purposes.

Structure and osmotic pressure of ionic microgel dispersions.

We investigate structural and thermodynamic properties of aqueous dispersions of ionic microgels--soft colloidal gel particles that exhibit unusual phase behavior. Starting from a coarse-grained model of microgel macroions as charged spheres that are permeable to microions, we perform simulations and theoretical calculations using two complementary implementations of Poisson-Boltzmann (PB) theory. Within a one-component model, based on a linear-screening approximation for effective electrostatic pair interactions, we perform molecular dynamics simulations to compute macroion-macroion radial distribution functions, static structure factors, and macroion contributions to the osmotic pressure. For the same model, using a variational approximation for the free energy, we compute both macroion and microion contributions to the osmotic pressure. Within a spherical cell model, which neglects macroion correlations, we solve the nonlinear PB equation to compute microion distributions and osmotic pressures. By comparing the one-component and cell model implementations of PB theory, we demonstrate that the linear-screening approximation is valid for moderately charged microgels. By further comparing cell model predictions with simulation data for osmotic pressure, we chart the cell model's limits in predicting osmotic pressures of salty dispersions.

Phase separation and the 'coffee-ring' effect in polymer-nanocrystal mixtures.

The coupling between the 'coffee-ring' effect and liquid-liquid phase separation is examined for ternary mixtures of solvent, polymer and semiconductor nanocrystal. Specifically, we study mixtures of toluene, polystyrene (PS) and colloidal silicon nanocrystals (SiNCs) using real-space imaging and spectroscopic techniques to resolve the kinetic morphology of the drying front for varied molecular weight of the PS. Our results demonstrate that the size of the polymer has a significant impact on both phase-separation and drying, and we relate these observations to simulations, measured and predicted binodal curves, and the observed shape of the flow field at the contact line. The results inform a deposition process that reduces the influence of drying instabilities for low-molecular-weight polymers while paving the way for more detailed and generic computational descriptions of drying instabilities in complex fluids.

Ion density deviations in polyelectrolyte microcapsules: influence on biosensors.

Polyelectrolyte microcapsules loaded with fluorescent dyes have been proposed as biosensors to monitor local pH and ionic strength for diagnostic purposes. In the case of charged microcapsules, however, the local electric field can cause deviations of ion densities inside the cavities, potentially resulting in misdiagnosis of some diseases. Using nonlinear Poisson-Boltzmann theory, we systematically investigate these deviations induced by charged microcapsules. Our results show that the microcapsule charge density, as well as the capsule and salt concentrations, contribute to deviations of local ion concentrations and pH. Our findings are relevant for applications of polyelectrolyte microcapsules with encapsulated ion-sensitive dyes as biosensors.