Ion permeation inside microgel particles induced by specific interactions: from charge inversion to overcharging

Colloidal Engineering of Microplastic Capture with Biodegradable Soft Dendritic "Microcleaners"

The introduction of colloidal principles that enable efficient microplastic collection from aquatic environments is a goal of great environmental importance. Here, we present a novel method of microplastic (MP) collection using biodegradable hydrogel soft dendritic colloids (hSDCs). These dendritic colloids have abundant nanofibrils and a large surface area, which provide an abundance of interfacial interactions and excellent networking capabilities, allowing for the capture of plastic particles and other contaminants. Here, we show how the polymer composition and morphology of the hSDCs can impact the capture of microplastics modeled by latex microbeads. Additionally, we use colloidal DLVO theory to interpret the capture efficiencies of microbeads of different sizes and surface functional groups. The results demonstrate the microplastic remediation efficiency of hydrogel dendricolloids and highlight the primary factors involved in the microbead interactions and adsorption. On a practical level, the results show that the development of environmentally benign microcleaners based on naturally sourced materials could present a sustainable solution for microplastic cleanup.

Hot Brownian Motion of Thermoresponsive Microgels in Optical Tweezers Shows Discontinuous Volume Phase Transition and Bistability

Microgels are soft microparticles that often exhibit thermoresponsiveness and feature a transformation at a critical temperature, referred to as the volume phase transition temperature. Whether this transformation occurs as a smooth or as a discontinuous one is still a matter of debate. This question can be addressed by studying individual microgels trapped in optical tweezers. For this aim, composite particles are obtained by decorating  Poly-N-isopropylacrylamide (pNIPAM) microgels with iron oxide nanocubes. These composites become self-heating when illuminated by the infrared trapping laser, performing hot Brownian motion within the trap. Above a certain laser power, a single decorated microgel features a volume phase transition that is discontinuous, while the usual continuous sigmoidal-like dependence is recovered after averaging over different microgels. The collective sigmoidal behavior enables the application of a power-to-temperature calibration and provides the effective drag coefficient of the self-heating microgels, thus establishing these composite particles as potential micro-thermometers and micro-heaters. Moreover, the self-heating microgels also exhibit an unexpected and intriguing bistability behavior above the critical temperature, probably due to partial collapses of the microgel. These results set the stage for further studies and the development of applications based on the hot Brownian motion of soft particles.

A self-consistent Ornstein-Zernike jellium for highly charged colloids (microgels) in suspensions with added salt

This work discusses a jellium scheme, built within the framework of the multicomponent Ornstein-Zernike (OZ) equation, which is capable of describing the collective structure of suspensions of highly charged colloids with added salt, even in the presence of finite-size multivalent microions. This approach uses a suitable approximation to decouple the microion-microion correlations from the macroion-microion profiles, which in combination with the methodology from the dressed ion theory (DIT) gives a full account of the electrostatic effective potential among the colloids. The main advantages of the present contribution reside in its ability to manage the short-range potentials and non-linear correlations among the microions, as well as its realistic characterization of the ionic clouds surrounding each macroion. The structure factors predicted by this jellium scheme are contrasted with previously reported experimental results for microgel suspensions with monovalent salts (2019Phys. Rev. E100032602), thus validating its high accuracy in these situations. The present theoretical analysis is then extended to microgel suspensions with multivalent salts, which reveals the prominent influence of the counterion valence on the makeup of the effective potentials. Although the induced differences may be difficult to identify through the mesoscopic structure, our results suggest that the microgel collapsing transition may be used to enhance these distinct effects, thus giving a feasible experimental probe for these phenomena.

Scaling Laws in the Diffusive Release of Neutral Cargo from Hollow Hydrogel Nanoparticles: Paclitaxel-Loaded Poly(4-vinylpyridine)

We study the nonequilibrium diffusive release of electroneutral molecular cargo encapsulated inside hollow hydrogel nanoparticles. We propose a theoretical model that includes osmotic, steric, and short-range polymer-cargo attractions to determine the effective cargo-hydrogel interaction, u_eff*, and the effective diffusion coefficient of the cargo inside the polymer network, D_eff*. Using dynamical density functional theory (DDFT), we investigate the scaling of the characteristic release time, τ_1/2, with the key parameters involved in the process, namely, u_eff*, D_eff*, and the swelling ratio. This effort represents a full study of the problem, covering a broad range of cargo sizes and providing predictions for repulsive and attractive polymer shells. Our calculations show that the release time through repulsive polymer networks scales with q²e^βu_eff*/D_eff* for βu_eff* ≫ 1. In this case, the cargo molecules are excluded from the shell of the hydrogel. For attractive shells, the polymer retains the cargo molecules on its internal surface and its interior, and the release time grows exponentially with the attraction strength. The DDFT calculations are compared to an analytical model for the mean first passage time, which provides an excellent quantitative description of the kinetics for both repulsive and attractive shells without fitting parameters. Finally, we apply the method to reproduce experimental results on the release of paclitaxel from hollow poly(4-vinylpyridine) nanoparticles and find that the slow release of the drug can be explained in terms of the strong binding attraction between the drug and the polymer.

Tuning the selective permeability of polydisperse polymer networks

We study the permeability and selectivity ('permselectivity') of model membranes made of polydisperse polymer networks for molecular penetrant transport, using coarse-grained, implicit-solvent computer simulations. In our work, permeability P is determined on the linear-response level using the solution-diffusion model, P = KDin, i.e., by calculating the equilibrium penetrant partition ratio K and penetrant diffusivity Din inside the membrane. We vary two key parameters, namely the network-network interaction, which controls the degree of swelling and collapse of the network, and the network-penetrant interaction, which tunes the selective penetrant uptake and microscopic energy landscape for diffusive transport. We find that the partitioning K covers four orders of magnitude and is a non-monotonic function of the parameters, well interpreted by a second-order virial expansion of the free energy of transferring one penetrant from a reservoir into the membrane. Moreover, we find that the penetrant diffusivity Din in the polydisperse networks, in contrast to highly ordered membrane structures, exhibits relatively simple exponential decays. We propose a semi-empirical scaling law for the penetrant diffusion that describes the simulation data for a wide range of densities and interaction parameters. The resulting permeability P turns out to follow the qualitative behavior (including maximization and minimization) of partitioning. However, partitioning and diffusion are typically anti-correlated, yielding large quantitative cancellations, controlled and fine-tuned by the network density and interactions, as rationalized by our scaling laws. We finally demonstrate that even small changes of network-penetrant interactions, e.g., by half a kBT, modify the permselectivity by almost one order of magnitude.

Volume transition effects on the correlations and effective interactions among highly charged microgels

Recent experimental studies have demonstrated the huge influence that the volume phase transition (VPT) has on the collective structure of highly charged thermo-responsive microgels in aqueous solution with low concentrations of added monovalent salt, thus opening a promising new route for controlling the overall properties of practical colloidal suspensions. We present here an analysis of this structure based on the effective electrostatic potential obtained with the exact methodology of the dressed ion theory (DIT). Starting with a description at the primitive model level, we determine the correlations among the components of our model system (macroions plus monovalent anions and cations) by utilizing the two-density integral equation theory, thus allowing us to consider realistic values for the microgel charges. The resulting microgel structure factors show a good agreement with the reported light scattering measurements, whereas the microscopic pair distributions reveal that in this regime the shrunken states promote an enhanced counterion absorption into the microgels. This packing of counterions inside the microgels induces strongly non-linear correlations among the microions, and in turn provokes a substantial weakening of the microgel-microgel correlations. The ensuing effective interactions are then obtained by contracting the description to the level in which only the macroions are present. We find not only that the magnitude and reach of the corresponding pair potentials are markedly inhibited in the shrunken states, but also that their general form diverges from the conventional screened Coulomb shape. This makes it necessary to rethink the concepts of effective charge and screening length.

Coarse-grained Monte Carlo simulations of nanogel-polyelectrolyte complexes: electrostatic effects

Coarse-grained Monte-Carlo simulations of nanogel-polyelectrolyte complexes have been carried out. The results presented here capture two phenomena reported in experiments with real complexes: (i) the reduction in size after absorbing just a few chains and (ii) the charge inversion detected through electrophoretic mobility data. Our simulations reveal that charge inversion occurs if the polyelectrolyte charge is large enough. In addition, the distribution of chains inside the nanogel strongly depends on whether charge inversion takes place. It should also be stressed that the chain topology has little influence on most of the properties studied here.

Accounting for effective interactions among charged microgels

We introduce a theoretical approach to describe structural correlations among charged permeable spheres at finite particle concentrations. This theory explicitly accounts for correlations among microions and between microions and macroions and allows for the proposal of an effective interaction among macroions that successfully captures structural correlations observed in poly-N-isopropyl acrylamide microgel systems. In our description the bare charge is fixed and independent of the microgel size, the microgel concentration, and the ionic strength, which contrasts with results obtained using linear response approximations, where the bare charge needs to be adapted to properly account for microgel correlations obtained at different conditions.

Tuning the Permeability of Dense Membranes by Shaping Nanoscale Potentials

Permeability is one of the most fundamental transport properties in soft matter physics, material engineering, and nanofluidics. Here, we report by means of Langevin simulations of ideal penetrants in a nanoscale membrane made of a fixed lattice of attractive interaction sites, how the permeability can be massively tuned, even minimized or maximized, by tailoring the potential energy landscape for the diffusing penetrants, depending on the membrane attraction, topology, and density. Supported by limiting scaling theories we demonstrate that the observed nonmonotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anticorrelation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.

Nonequilibrium Uptake Kinetics of Molecular Cargo into Hollow Hydrogels Tuned by Electrosteric Interactions

Hollow hydrogels represent excellent nano- and microcarriers due to their ability to encapsulate and release large amounts of cargo molecules (cosolutes) such as reactants, drugs, and proteins. In this work, we use a combination of a phenomenological effective cosolute-hydrogel interaction potential and dynamic density functional theory to investigate the full nonequilibrium encapsulation kinetics of charged and dipolar cosolutes by an isolated charged hollow hydrogel immersed in a 1:1 electrolyte aqueous solution. Our analysis covers a broad spectrum of cosolute valences ( z_c) and electric dipole moments (μ_c), as well as hydrogel swelling states and hydrogel charge densities. Our calculations show that, close to the collapsed state, the polar cosolutes are predominantly precluded and the encapsulation process is strongly hindered by the excluded-volume interaction exerted by the polymer network. Different equilibrium and kinetic sorption regimes (interface versus interior) are found depending on the value and sign of z_c and the value of μ_c. For cosolutes of the same sign of charge as the gel, the superposition of steric and electrostatic repulsion leads to an "interaction-controlled" encapsulation process, in which the characteristic time to fill the empty core of the hydrogel grows exponentially with z_c. On the other hand, for cosolutes oppositely charged to the gel, we find a "diffusion-controlled" kinetic regime, where cosolutes tend to rapidly absorb into the hydrogel membrane and the encapsulation rate depends only on the cosolute diffusion time across the membrane. Finally, we find that increasing μ_c promotes the appearance of metastable and stable surface adsorption states. For large enough μ_c, the kinetics enters an "adsorption-hindered diffusion", where the enhanced surface adsorption imposes a barrier and slows down the uptake. Our study represents the first attempt to systematically describe how the swelling state of the hydrogel and other leading physical interaction parameters determine the encapsulation kinetics and the final equilibrium distribution of polar molecular cargo.

Tuning the collapse transition of weakly charged polymers by ion-specific screening and adsorption

The experimentally observed swelling and collapse response of weakly charged polymers to the addition of specific salts displays quite convoluted behavior that is not easy to categorize. Here we use a minimalistic implicit-solvent/explicit-salt simulation model with a focus on ion-specific interactions between ions and a single weakly charged polyelectrolyte to qualitatively explain the observed effects. In particular, we demonstrate ion-specific screening and bridging effects cause collapse at low salt concentrations whereas the same strong ion-specific direct interactions drive re-entrant swelling at high concentrations. Consistently with experiments, a distinct salt concentration at which the salting-out power of anions inverts from the reverse to direct Hofmeister series is observed. At this so called isospheric point, the ion-specific effects vanish. Furthermore, with additional simplifying assumptions, an ion-specific mean-field model is developed for the collapse transition which quantitatively agrees with the simulations. Our work demonstrates the sensitivity of the structural behavior of charged polymers to the addition of specific salt beyond simple screening and shall be useful for further guidance of experiments.

Conformation change of an isotactic poly (N-isopropylacrylamide) membrane: Molecular dynamics

In this work, isotactic Poly (N-Isopropylacrylamide)-PNIPAM-in neat water and in electrolyte solutions is studied by means of molecular dynamics simulations. This is done for an infinitely diluted oligomer and for an assembly of several PNIPAM chains arranged into a planar membrane configuration with a core-shell morphology. We employed two different force fields, AMBER (assisted model building with energy refinement) and OPLS-AA (all atom - optimized potentials for liquid simulations) in combination with extended simple point charge water. Despite the more water insoluble character of isotactic oligomers, our results support the existence of a coil to globule transition for the isolated 30-mer. This may imply the existence of an oligomer rich phase of coil-like structures in equilibrium with a water rich phase for temperatures close but below the coil to globule transition temperature, T_Θ. However, the obtained coil structure is much more compact than that corresponding to the syndiotactic chain. Our estimations of T_Θ are (308±5) K and (303±5) K for AMBER and OPLS-AA, respectively. The membrane configuration allows one to include chain-chain interactions, to follow density profiles of water, polymer, and solutes, and accessing the membrane-water interface tension. Results show gradual shrinking and swelling of the membrane by switching temperature above and below T_Θ, as well as the increase and decrease of the membrane-water interface tension. Finally, concentration profiles for 1M NaCl and 1M NaI electrolytes are shown, depicting a strong salting-out effect for NaCl and a much lighter effect for NaI, in good qualitative agreement with experiments.

Maximizing the absorption of small cosolutes inside neutral hydrogels: steric exclusion versus hydrophobic adhesion

In this work the equilibrium absorption of nanometric cosolutes (which could represent drugs, reactants, small globular proteins and other kind of biomacromolecules) inside neutral hydrogels is studied. We specially focus on exploring, for different swelling states, the competition between the steric exclusion induced by the cross-linked polymer network constituting the hydrogel, and the solvent-induced short-range hydrophobic attraction between the polymer chains and the cosolute particle. For this purpose, the cosolute partition coefficient is calculated by means of coarse-grained grand canonical Monte Carlo simulations, and the results are compared to theoretical predictions based on the calculation of the excluded and binding volume around the polymer chains. For small hydrophobic attractions or large cosolute sizes, the steric repulsion dominates, and the partition coefficient decreases monotonically with the polymer volume fraction, ϕ_m. However, for large enough hydrophobic attraction strength, the interplay between hydrophobic adhesion and the steric exclusion leads to a maximum in the partition coefficient at certain intermediate polymer density. Good qualitative and quantitative agreement is achieved between simulation results and theoretical predictions in the limit of small ϕ_m, pointing out the importance of geometrical aspects of the cross-linked polymer network, even for hydrogels in the swollen state. In addition, the theory is able to predict analytically the onset of the maximum formation in terms of the details of the cosolute-monomer pair interaction, in good agreement with simulations too. Finally, the effect of the many-body attractions between the cosolute and multiple polymer chains is quantified. The results clearly show that these many-body attractions play a very relevant role determining the cosolute binding, enhancing its absorption in more than one order of magnitude.

On the physics of both surface overcharging and charge reversal at heterophase interfaces

A series of Monte Carlo simulations are employed to reveal the physics of both surface overcharging and charge reversal at a negatively charged dielectric interface exposed to a bulk solution containing a +2:−1 electrolyte in the absence and presence of a monovalent salt.

Catalyzed Bimolecular Reactions in Responsive Nanoreactors

We describe a general theory for surface-catalyzed bimolecular reactions in responsive nanoreactors, catalytically active nanoparticles coated by a stimuli-responsive "gating" shell, whose permeability controls the activity of the process. We address two archetypal scenarios encountered in this system: the first, where two species diffusing from a bulk solution react at the catalyst's surface, and the second, where only one of the reactants diffuses from the bulk while the other is produced at the nanoparticle surface, e.g., by light conversion. We find that in both scenarios the total catalytic rate has the same mathematical structure, once diffusion rates are properly redefined. Moreover, the diffusional fluxes of the different reactants are strongly coupled, providing a behavior richer than that arising in unimolecular reactions. We also show that, in stark contrast to bulk reactions, the identification of a limiting reactant is not simply determined by the relative bulk concentrations but is controlled by the nanoreactor shell permeability. Finally, we describe an application of our theory by analyzing experimental data on the reaction between hexacyanoferrate(III) and borohydride ions in responsive hydrogel-based core-shell nanoreactors.

Competition between excluded-volume and electrostatic interactions for nanogel swelling: effects of the counterion valence and nanogel charge

The equilibrium distribution of monovalent and trivalent ions within a thermo-responsive charged nanogel is investigated using Monte Carlo simulations and Ornstein–Zernike equations.

Swelling of ionic microgel particles in the presence of excluded-volume interactions: a density functional approach

In this work a new density functional theory framework is developed to predict the salt-concentration dependent swelling state of charged microgels and the local concentration of monovalent ions inside and outside the microgel.

Temperature-sensitive nanogels in the presence of salt: explicit coarse-grained simulations

In this work, coarse-grained simulations of two charged thermo-shrinking nanogels (with degrees of ionization of 0.125 and 0.250) in the presence of 1:1 and 3:1 electrolytes have been explicitly performed through the bead-spring model of polyelectrolyte. In a first set of simulations, salt concentrations for 1:1 and 3:1 electrolytes ranged from 1 to 100 mM and from 0.167 to 16.7 mM, respectively, whereas temperature remained fixed at a value for which hydrophobic forces were negligible in our case (288 K). The sizes of swollen nanogels are smaller when trivalent cations are present, but they do not change significantly in the range of concentrations of 3:1 electrolyte studied here. It should be also stressed that trivalent cations neutralize the nanogel charge more efficiently. According to these results the electrostatic repulsion plays an important role. In a second set of simulations, the temperature varied from 288 to 333 K to study the effect of salt on the thermal response when hydrophobic forces are not negligible. For the nanogels with the lowest degree of ionization, the behavior of the radius with increasing the temperature can be described by a sigmoid function, which shifts towards lower temperatures in the presence of salt. This shift is more clearly observed for trivalent cations, even at low concentrations. For the nanogels with the highest degree of ionization, the effect of additional electrolyte is also noticeable. In this case, hydrophobic forces are not the only responsible for their shrinkage in the presence of trivalent cations. The surface electrostatic potential and the concentration of salt cations inside the nanogel have been computed from simulations and a modified Poisson-Boltzmann (PB) cell model. The thermosensitivity in size have certain influence on the sensitivity of these properties to temperature changes. The rich behavior of the surface electrostatic potential and the uptake of salt cations are successfully predicted by the modified PB cell model proposed (at least qualitatively). Particularly, the model is able to predict how the retention of salt cations depends on their charge and the ionic valence when nanogels shrink.

Excluded volume effects on ionic partitioning in gels and microgels: a simulation study

An analytical expression accounting for excluded volume effects on ionic partitioning in gels is tested through Monte Carlo simulations.