Swelling, structure, and phase stability of compressible microgels

Structural transitions of ionic microgel solutions driven by circularly polarized electric fields

In this work, a theoretical approach is developed to investigate the structural properties of ionic microgels induced by a circularly polarized (CP) electric field. Following a similar study on chain formation in the presence of linearly polarized fields [T. Colla et al., ACS Nano, 2018, 12, 4321-4337], we propose an effective potential between microgels which incorporates the field-induced interactions via a static, time averaged polarizing charge at the particle surface. In such a coarse-graining framework, the induced dipole interactions are controlled by external parameters such as the field strength and frequency, ionic strength, as well as microgel charge and concentration, thus providing a convenient route to induce different self-assembly scenarios through experimentally adjustable quantities. In contrast to the case of linearly polarized fields, dipole interactions in the case of CP light are purely repulsive in the direction perpendicular to the polarization plane, while featuring an in-plane attractive well. As a result, the CP field induces layering of planar sheets arranged perpendicularly to the field direction, in strong contrast to the chain formation observed in the case of linear polarizations. Depending on the field strength and particle concentration, in-plane crystallization can also take place. Combining molecular dynamics (MD) simulations and the liquid-state hypernetted-chain (HNC) formalism, we herein investigate the emergence of layering formation and in-plane crystal ordering as the dipole strength and microgel concentration are changed over a wide region of parameter space.

A multi-body finite element model for hydrogel packings: linear response to shear

We study a multi-body finite element model of a packing of hydrogel particles using the Flory-Rehner constitutive law to model the deformation of the swollen polymer network. We show that while the dependence of the pressure, Π, on the effective volume fraction, ϕ, is virtually identical to a monolithic Flory material, the shear modulus, μ, behaves in a non-trivial way. μ increases monotonically with Π from zero and remains below about 80% of the monolithic Flory value at the largest Π we study here. The local shear strain in the particles has a large spatial variation. Local strains near the centers of the particles are all roughly equal to the applied shear strain, but the local strains near the contact facets are much smaller and depend on the orientation of the facet. We show that the slip between particles at the facets depends strongly on the orientation of the facet and is, on average, proportional to the component of the applied shear strain resolved onto the facet orientation. This slip screens the stress transmission and results in a reduction of the shear modulus relative to what one would obtain if the particles were welded together at the facet. Surprisingly, given the reduction in the shear modulus arising from the facet slip, and the spatial variations in the local shear strain inside the particles themselves, the deformation of the particle centroids is rather homogeneous with the strains of the Delaunay triangles having fluctuations of only order ±5%. These results should open the way to construction of quantitative estimates of the shear modulus in highly compressed packings via mean-field, effective-medium type approaches.

Nonequilibrium relaxation of soft responsive colloids

Stimuli-responsive macromolecules display large conformational changes during their dynamics, sometimes switching between states. Such a multi-stability is useful for the development of soft functional materials. Here, we introduce a mean-field dynamical density functional theory for a model of responsive colloids to study the nonequilibrium dynamics of a colloidal dispersion in time-dependent external fields, with a focus on the coupling of translational and conformational dynamics during their relaxation. Specifically, we consider soft Gaussian particles with a bimodal size distribution between two confining walls with time-dependent (switching-on and off) external gravitational and osmotic fields. We find a rich relaxation behavior of the systems in excellent agreement with particle-based Brownian dynamics computer simulations. In particular, we find time-asymmetric relaxations of integrated observables (wall pressures, mean size, and liquid center-of-mass) for activation/deactivation of external potentials, respectively, which are tunable by the ratio of translational and conformational diffusion time scales. Our work thus paves the way for studying the nonequilibrium relaxation dynamics of complex soft matter with multiple degrees of freedom and hierarchical relaxations.

External field-driven property localization in liquids of responsive macromolecules

We explore theoretically the effects of external potentials on the spatial distribution of particle properties in a liquid of explicitly responsive macromolecules. In particular, we focus on the bistable particle size as a coarse-grained internal degree of freedom (DoF, or "property"), σ, that moves in a bimodal energy landscape, in order to model the response of a state-switching (big-to-small) macromolecular liquid to external stimuli. We employ a mean-field density functional theory (DFT) that provides the full inhomogeneous equilibrium distributions of a one-component model system of responsive colloids (RCs) interacting with a Gaussian pair potential. For systems confined between two parallel hard walls, we observe and rationalize a significant localization of the big particle state close to the walls, with pressures described by an exact RC wall theorem. Application of more complex external potentials, such as linear (gravitational), osmotic, and Hamaker potentials, promotes even stronger particle size segregation, in which macromolecules of different size are localized in different spatial regions. Importantly, we demonstrate how the degree of responsiveness of the particle size and its coupling to the external potential tune the position-dependent size distribution. The DFT predictions are corroborated by Brownian dynamics simulations. Our study highlights the fact that particle responsiveness can be used to localize liquid properties and therefore helps to control the property- and position-dependent function of macromolecules, e.g., in biomedical applications.

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.

Poly(N-isopropylacrylamide) microgel swelling behavior and suspension structure studied with small-angle neutron scattering

Microgels are of high interest for applications and as model systems due to their volume response to external stimuli. We use small-angle neutron scattering to measure the form and structure factors of poly(N-isopropylacrylamide) microgels in dilute and concentrated suspensions and find that microgels keep a constant size up to a concentration, above which they deswell. This happens before random-close packing. We emphasize suspension polydispersity must be considered to obtain accurate form and structure factors. Our results are compatible with microgel deswelling triggered by the osmotic pressure set by counterions associated to charged groups in the microgel periphery, which sharply increases when the counterion clouds surrounding the microgels percolate throughout the suspension volume.

Measuring the counterion cloud of soft microgels using SANS with contrast variation

The behavior of microgels and other soft, compressible colloids depends on particle concentration in ways that are absent in their hard-particulate counterparts. For instance, poly-N-isopropylacrylamide (pNIPAM) microgels can spontaneously deswell and reduce suspension polydispersity when concentrated enough. Despite the pNIPAM network in these microgels is neutral, the key to understanding this distinct behavior relies on the existence of peripheric charged groups, responsible for providing colloidal stability when deswollen, and the associated counterion cloud. When in close proximity, clouds of different particles overlap, effectively freeing the associated counterions, which are then able to exert an osmotic pressure that can potentially cause the microgels to decrease their size. Up to now, however, no direct measurement of such an ionic cloud exists, perhaps even also for hard colloids, where it is referred to as an electric double layer. Here, we use small-angle neutron scattering with contrast variation with different ions to isolate the change in the form factor directly related to the counterion cloud, and obtain its radius and width. Our results highlight that the modeling of microgel suspensions must unavoidably and explicitly consider the presence of this cloud, which exists for nearly all microgels synthesized today.

Bragg-Williams Theory for Particles with a Size-Modulating Internal Degree of Freedom

The field of soft matter teems with molecules and aggregates of molecules that have internal size-modulating degrees of freedom. Proteins, peptides, microgels, polymers, micelles, and even some colloids can exist in multiple-often just two dominating-states with different effective sizes, where size can refer to the volume or to the cross-sectional area for particles residing on surfaces. The size-dependence of their accessible states renders the behavior of these particles pressure-sensitive. The Bragg-Williams model is among the most simple mean-field methods to translate the presence of inter-particle interactions into an approximate phase diagram. Here, we extend the Bragg-Williams model to account for the presence of particles that are immersed in a solvent and exist in two distinct states, one occupying a smaller and the other one a larger size. The basis of the extension is a lattice-sublattice approximation that we use to host the two size-differing states. Our model includes particle-solvent interactions that act as an effective surface tension between particles and solvent and are ignorant of the state in which the particles reside. We analyze how the energetic preference of the particles for one or the other state affects the phase diagrams. The possibility of a single phase-two phases-single phase sequence of phase transitions as a function of increasing temperature is demonstrated.

Liquid structure of bistable responsive macromolecules using mean-field density-functional theory

Macromolecular crowding typically applies to biomolecular and polymer-based systems in which the individual particles often feature a two-state folded/unfolded or coil-to-globule transition, such as found for proteins and peptides, DNA and RNA, or supramolecular polymers. Here, we employ a mean-field density functional theory (DFT) of a model of soft and bistable responsive colloids (RCs) in which the size of the macromolecule is explicitly resolved as a degree of freedom living in a bimodal 'Landau' energy landscape (exhibiting big and small states), thus directly responding to the crowding environment. Using this RC-DFT we study the effects of self-crowding on the liquid bulk structure and thermodynamics for different energy barriers and softnesses of the bimodal energy landscape, in conditions close to the coil-to-globule transition. We find substantial crowding effects on the internal distributions, a complex polydispersity behavior, and quasi-universal compression curves for increasing (generalized) packing fractions. Moreover, we uncover distinct signatures of bimodal versus unimodal behavior in the particle compression. Finally, the analysis of the pair structure - derived from the test particle route - reveals that the microstructure of the liquid is quite inhomogeneous due to local depletion effects, tuneable by particle softness.

Active responsive colloids driven by intrinsic dichotomous noise

We study the influence of intrinsic noise on the structure and dynamics of responsive colloids (RCs), which actively change their size and mutual interactions. The colloidal size is explicitly resolved in our RC model as an internal degree of freedom (DOF) in addition to the particle translation. A Hertzian pair potential between the RCs leads to repulsion and shrinking of the particles, resulting in an explicit responsiveness of the system to self-crowding. To render the colloids active, their size is internally driven by a dichotomous noise, randomly switching ("breathing") between growing and shrinking states with a predefined rate, as motivated by recent experiments on synthetic active colloids. The polydispersity of this dichotomous active responsive colloid (D-ARC) model can be tuned by the parameters of the noise. Utilizing stochastic computer simulations, we study crowding effects on the spatial distributions, relaxation times, and self-diffusion of dense suspensions of the D-ARCs. We find a substantial influence of the "built-in" intrinsic noise on the system's behavior, in particular, transitions from unimodal to bimodal size distributions for an increasing colloid density as well as intrinsic noise-modified diffusive translational dynamics. We conclude that controlling the noise of internal DOFs of a macromolecule or cell is a powerful tool for active colloidal materials to enable autonomous changes in the system's collective structure and dynamics towards the adaptation of macroscopic properties to external perturbations.

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.

Active responsive colloids coupled to different thermostats

We introduce a model of active responsive colloids (ARCs) in which an internal degree of freedom (DoF) of a single colloidal particle is "activated" by coupling it to a different thermostat than for the translational DoFs. As for the responsive internal DoF, we consider specifically the size (diameter) of the spherical particles, which is confined by a harmonic parent potential being either entropic or energetic in nature. The ARCs interact via a repulsive Hertzian pair potential, appropriate to model hydrogels or elastic colloids, and are studied for various densities using Brownian dynamics simulations. We tune the internal activity in the nonequilibrium steady state by scanning through a wide range of internal temperatures, both smaller ("colder") and larger ("hotter") than the translational temperature. The results show a rich and intriguing behavior for the emergent property distributions, colloidal pair structure, and the diffusive translational dynamics controlled by the internal activity, substantially depending on whether the internal DoF moves in an entropic or energetic potential. We discuss theoretical thermal limits and phenomenological models which can explain some of the nonequilibrium trends qualitatively. Our study indicates that particle dynamical polydispersity as well as the structure and dynamics of dense macromolecular suspensions can be vastly tuned by internal activity in terms of internal "hot" or "cold" fluctuating states.

How Softness Matters in Soft Nanogels and Nanogel Assemblies

Softness plays a key role in determining the macroscopic properties of colloidal systems, from synthetic nanogels to biological macromolecules, from viruses to star polymers. However, we are missing a way to quantify what the term "softness" means in nanoscience. Having quantitative parameters is fundamental to compare different systems and understand what the consequences of softness on the macroscopic properties are. Here, we propose different quantities that can be measured using scattering methods and microscopy experiments. On the basis of these quantities, we review the recent literature on micro- and nanogels, i.e. cross-linked polymer networks swollen in water, a widely used model system for soft colloids. Applying our criteria, we address the question what makes a nanomaterial soft? We discuss and introduce general criteria to quantify the different definitions of softness for an individual compressible colloid. This is done in terms of the energetic cost associated with the deformation and the capability of the colloid to isotropically deswell. Then, concentrated solutions of soft colloids are considered. New definitions of softness and new parameters, which depend on the particle-to-particle interactions, are introduced in terms of faceting and interpenetration. The influence of the different synthetic routes on the softness of nanogels is discussed. Concentrated solutions of nanogels are considered and we review the recent results in the literature concerning the phase behavior and flow properties of nanogels both in three and two dimensions, in the light of the different parameters we defined. The aim of this review is to look at the results on micro- and nanogels in a more quantitative way that allow us to explain the reported properties in terms of differences in colloidal softness. Furthermore, this review can give researchers dealing with soft colloids quantitative methods to define unambiguously which softness matters in their compound.

Transition in the Glassy Dynamics of Melts of Acid-Hydrolyzed Phytoglycogen Nanoparticles

The deformability, responsiveness, and tunability of soft nanoparticles (NPs) offer unique opportunities to learn about their complex properties and the interactions between particles. In the present study, we provide new insights into the physical properties of phytoglycogen (PG) NPs, which are soft, compact particles with a dendritic architecture that are produced in the kernels of sweet corn. In particular, we study PG NPs modified using acid hydrolysis, which not only reduces their diameter but also alters their stiffness, internal structure, and the interactions between particles in aqueous dispersions. We used steady shear rheology to determine the dependence of the relative zero-shear viscosity η_r of aqueous dispersions of acid-hydrolyzed PG NPs on the effective volume fraction ϕ_eff, which indicated a reduction in stiffness of the particles relative to that of native PG NPs. We quantified this difference by analyzing the nature of the colloidal glasses formed at high ϕ_eff. We measured a smaller value of the fragility index m for acid-hydrolyzed PG NP glasses than that for native PG NP glasses, indicating that acid-hydrolyzed PG NPs form stronger glasses and are therefore softer than native PG NPs. Unlike the native PG NPs, we observed a distinctive change in the character of the glass transition of the acid-hydrolyzed PG NPs as ϕ_eff was increased above ϕ_eff∼1: a crossover in the dependence of η_r on ϕ_eff from Vogel-Fulcher-Tammann behavior to a more gradual, Arrhenius-like behavior. By expressing the steady shear and oscillatory rheology data in terms of generalized Péclet numbers, we obtained collapse of the data onto master curves. We interpret this result in terms of the acid-hydrolyzed PG NPs predominantly interpenetrating neighboring particles at large ϕ_eff, for which fluctuations of the outer chains enhance the mobility of the particles and make α-relaxation times τ_α experimentally accessible.

Changes in the Form Factor and Size Distribution of Nanogels in Crowded Environments

Particle size disparities suppress crystallization. However, soft deformable nanogels can change the size of the larger particles in suspension and crystallize even at a high initial size-polydispersity. Using neutron scattering with contrast variation, the response of individual nanogels in crowded environments was probed, and an increase of the parameter describing size-polydispersity was found, which is often interpreted as deformation. Here, computer simulations are used to generate deformed nanogels and the corresponding form factor. The data are fitted with the spherical model used to analyze scattering data. The fits show the same qualitative increase of the parameter related to the size-polydispersity with increasing particle deformation. Starting from the simulated deformed spheres, we also reproduce experimental scattering data. A further analysis of the particle shows that the size disparities between nanogels do not increase significantly. In contrast, their shapes strongly vary from one nanogel to the other.

Modulating internal transition kinetics of responsive macromolecules by collective crowding

Packing and crowding are used in biology as mechanisms to (self-)regulate internal molecular or cellular processes based on collective signaling. Here, we study how the transition kinetics of an internal "switch" of responsive macromolecules is modified collectively by their spatial packing. We employ Brownian dynamics simulations of a model of Responsive Colloids, in which an explicit internal degree of freedom-here, the particle size-moving in a bimodal energy landscape self-consistently responds to the density fluctuations of the crowded environment. We demonstrate that populations and transition times for the two-state switching kinetics can be tuned over one order of magnitude by "self-crowding." An exponential scaling law derived from a combination of Kramers' and liquid state perturbation theory is in very good agreement with the simulations.

Softness mapping of the concentration dependence of the dynamics in model soft colloidal systems

The dynamics of a series of soft colloids comprised of polystyrene cores with poly(N-isopropylacrylamide) (PNIPAM) coronas was investigated by diffusing wave spectroscopy (DWS). The modulus of the coronas was varied by changing the cross-link density and we were able to interpret the results within a hard-soft mapping framework. The soft, swellable particle properties were modeled using an extended Flory-Rehner theory and a Hertzian pair potential. Following volume fraction jumps, softness effects on the concentration dependence of dynamics were determined, with a 'soft colloids make strong glass-forming liquid'-type of behavior observed close to the nominal glass transition volume fraction, φ_g. Such behavior from the current systems cannot be fully explained by the osmotic deswelling model alone. However, inspired by the soft-hard mapping from Schmiedeberg et al, [Europhys. Lett. 2011, 96(3), 36010] we estimated effective hard-sphere diameters and achieved a successful mapping of the α-relaxation times to a master curve below φ_g. Above φ_g, the curves no longer collapse but show strong deviations from a Vogel-Fulcher type of divergence onto soft jamming plateaux. Our results provide evidence that osmotic deswelling itself cannot fully explain the observed dynamics. Softness also plays an important role in the dynamics of soft, concentrated colloids.

Volume fraction determination of microgel composed of interpenetrating polymer networks of PNIPAM and polyacrylic acid

Interpenetrated polymer network microgels, composed of crosslinked networks of poly(N-isopropylacrylamide) and polyacrylic acid (PAAc), have been investigated through rheological measurements at four different amounts of PAAc. Both PAAc content and crosslinking degree modify particle dimensions, mass and softness, thereby strongly affecting the volume fraction and the system viscosity. Here the volume fraction is derived from the flow curves at low concentrations by fitting the zero-shear viscosity with the Einstein-Batchelor equation which provides a parameterkto shift weight concentration to volume fraction. We find that particles with higher PAAc content and crosslinker are characterized by a greater value ofkand therefore by larger volume fractions when compared to softer particles. The packing fractions obtained from rheological measurements are compared with those from static light scattering for two PAAc contents revealing a good agreement. Moreover, the behaviour of the viscosity as a function of packing fraction, at room temperature, has highlighted an Arrhenius dependence for microgels synthesized with low PAAc content and a Vogel-Fulcher-Tammann dependence for the highest investigated PAAc concentration. A comparison with the hard spheres behaviour indicates a steepest increase of the viscosity with decreasing particles softness. Finally, the volume fraction dependence of the viscosity at a fixed PAAc and at two different temperatures, below and above the volume phase transition, shows a quantitative agreement with the structural relaxation time measured through dynamic light scattering indicating that interpenetrated polymer network microgels softness can be tuned with PAAc and temperature and that, depending on particle softness, two different routes are followed.

Dynamical properties of different models of elastic polymer rings: Confirming the link between deformation and fragility

We report extensive numerical simulations of different models of 2D polymer rings with internal elasticity. We monitor the dynamical behavior of the rings as a function of the packing fraction to address the effects of particle deformation on the collective response of the system. In particular, we compare three different models: (i) a recently investigated model [N. Gnan and E. Zaccarelli, Nat. Phys. 15, 683 (2019)] where an inner Hertzian field providing the internal elasticity acts on the monomers of the ring, (ii) the same model where the effect of such a field on the center of mass is balanced by opposite forces, and (iii) a semi-flexible model where an angular potential between adjacent monomers induces strong particle deformations. By analyzing the dynamics of the three models, we find that in all cases, there exists a direct link between the system fragility and particle asphericity. Among the three, only the first model displays anomalous dynamics in the form of a super-diffusive behavior of the mean-squared displacement and of a compressed exponential relaxation of the density auto-correlation function. We show that this is due to the combination of internal elasticity and the out-of-equilibrium force self-generated by each ring, both of which are necessary ingredients to induce such a peculiar behavior often observed in experiments of colloidal gels. These findings reinforce the role of particle deformation, connected to internal elasticity, in driving the dynamical response of dense soft particles.

Structure and dynamics of responsive colloids with dynamical polydispersity

Dynamical polydispersity in single-particle properties, for example a fluctuating particle size, shape, charge density, etc, is intrinsic to responsive colloids (RCs), such as biomacromolecules or microgels, but is typically not resolved in coarse-grained mesoscale simulations. Here, we present Brownian dynamics simulations of suspensions of RCs modeling soft hydrogel colloids, for which the size of the individual particles is an explicitly resolved (Gaussian) degree of freedom and dynamically responds to the local interacting environment. We calculate the liquid structure, emergent size distributions, long-time diffusion, and property (size) relaxation kinetics for a wide range of densities and intrinsic property relaxation times in the canonical ensemble. Comparison to interesting reference cases, such as conventional polydisperse suspensions with a frozen parent distribution, or conventional monodisperse systems interacting with an effective pair potential for one fixed size, shows a significant spread in the structure and dynamics. The differences, most apparent in the high density regimes, are due to many-body correlations and the dynamical coupling between property and translation in RC systems, not explicitly accounted for in the conventional treatments. In particular, the translational diffusion in the RC systems is surprisingly close to the free (single RC) diffusion, mainly due to a cancellation of crowding and size compression effects. We show that an effective monodisperse pair potential can be constructed that describes the many-body correlations reasonably well by convoluting the RC pair potential with the density-dependent emergent size distributions and using a mean effective diffusion constant.

Internal structure of ultralow-crosslinked microgels: From uniform deswelling to phase separation

We perform small angle neutron scattering on ultralow-crosslinked microgels and find that while in certain conditions both the particle size and the characteristic internal length scale change in unison, in other instances this is not the case. We show that nonuniform deswelling depends not only on particle size, but also on the particular way the various contributions to the free energy combine to result in a given size. Only when polymer-solvent demixing strongly competes with ionic or electrostatic effects do we observe nonuniform behavior, reflecting internal microphase separation. The results do not appreciably depend on particle number density; even in concentrated suspensions, we find that at relatively low temperature, where demixing is not very strong, the deswelling behavior is uniform, and that only at sufficiently high temperature, where demixing is very strong, does the microgel structure change akin to internal microphase separation.

Effect of the Cross-Linking Density on the Swelling and Rheological Behavior of Ester-Bridged β-Cyclodextrin Nanosponges

The cross-linking density influences the physicochemical properties of cyclodextrin-based nanosponges (CD-NSs). Although the effect of the cross-linker type and content on the NSs performance has been investigated, a detailed study of the cross-linking density has never been performed. In this contribution, nine ester-bridged NSs based on β-cyclodextrin (β-CD) and different quantities of pyromellitic dianhydride (PMDA), used as a cross-linking agent in stoichiometric proportions of 2, 3, 4, 5, 6, 7, 8, 9, and 10 moles of PMDA for each mole of CD, were synthesized and characterized in terms of swelling and rheological properties. The results, from the swelling experiments, exploiting Flory-Rehner theory, and rheology, strongly showed a cross-linker content-dependent behavior. The study of cross-linking density allowed to shed light on the efficiency of the synthesis reaction methods. Overall, our study demonstrates that by varying the amount of cross-linking agent, the cross-linked structure of the NSs matrix can be controlled effectively. As PMDA βCD-NSs have emerged over the years as a highly versatile class of materials with potential applications in various fields, this study represents the first step towards a full understanding of the correlation between their structure and properties, which is a key requirement to effectively tune their synthesis reaction in view of any specific future application or industrial scale-up.

Structure of the simple harmonic-repulsive system in liquid and glassy states studied by the triple correlation function

An efficient description of the structures of liquids and, in particular, the structural changes that happen with liquids on supercooling remains to be a challenge. The systems composed of soft particles are especially interesting in this context because they often demonstrate non-trivial local orders that do not allow to introduce the concept of the nearest-neighbor shell. For this reason, the use of some methods, developed for the structure analysis of atomic liquids, is questionable for the soft-particle systems. Here we report about our investigations of the structure of the simple harmonic-repulsive liquid in 3D using the triple correlation function (TCF), i.e., the method that does not rely on the nearest neighbor concept. The liquid is considered at reduced pressure (P = 1.8) at which it exhibits remarkable stability against crystallization on cooling. It is demonstrated that the TCF allows addressing the development of the orientational correlations in the structures that do not allow drawing definite conclusions from the studies of the bond-orientational order parameters. Our results demonstrate that the orientational correlations, if measured by the heights of the peaks in the TCF, significantly increase on cooling. This rise in the orientational ordering is not captured properly by the Kirkwood's superposition approximation. Detailed considerations of the peaks' shapes in the TCF suggest the existence of a link between the orientational ordering and the slowdown of the system's dynamics. Our findings support the view that the development of the orientational correlations in liquids may play a significant role in the liquids' dynamics and that the considerations of the pair distribution function may not be sufficient to understand intuitively all the structural changes that happen with liquids on supercooling. In general, our results demonstrate that the considerations of the TCF are useful in the discussions of the liquid's structures beyond the pair density function and interpreting the results obtained with the bond-orientational order parameters.

Osmotic pressure of suspensions comprised of charged microgels

We determine the osmotic pressure of microgel suspensions using membrane osmometry and dialysis, for microgels with different softnesses. Our measurements reveal that the osmotic pressure of solutions of both ionic and neutral microgels is determined by the free ions that leave the microgel periphery to maximize their entropy and not by the translational degrees of freedom of the microgels themselves. Furthermore, up to a given concentration it is energetically favorable for the microgels to maintain a constant volume without appreciable deswelling. The concentration where deswelling starts weakly depends on the crosslinker concentration, which affects the microgel dimension; we explain this by considering the dependence of the osmotic pressure and the microgel bulk modulus on the particle size.

Phase behavior of ultrasoft spheres show stable bcc lattices

The phase behavior of supersoft spheres is explored using solutions of ultralow cross-linked poly(N-isopropylacrylamide)-based microgels as a model system. For these microgels, the effects of the electric charges on their surfaces can be neglected and therefore only the role of softness on the phase behavior is investigated. The samples show a liquid-to-crystal transition at higher volume fraction with respect to both hard spheres and stiffer microgels. Furthermore, stable body centered cubic (bcc) crystals are observed in addition to the expected face centered cubic (fcc) crystals. Small-angle x-ray and neutron scattering with contrast variation allow the characterization of both the microgel-to-microgel distance and the architecture of single microgels in crowded solutions. The measurements reveal that the stable bcc crystals depend on the interplay between the collapse and the interpenetration of the external shell of the ultralow cross-linked microgels.

Structure and position-dependent properties of inhomogeneous suspensions of responsive colloids

Responsive particles, such as biomacromolecules or hydrogels, display a broad and polymodal distribution of conformations and have thus the ability to change their properties (e.g., size, shape, charge density, etc.) substantially in response to external fields or to their local environment (e.g., mediated by cosolutes or pH). Here we discuss the basic statistical mechanics for a model of responsive colloids (RCs) by introducing an additional "property" degree of freedom as a collective variable in a formal coarse-graining procedure. The latter leads to an additional one-body term in the coarse-grained (CG) free energy, defining a single-particle property distribution for an individual polydisperse RC. We argue that in the equilibrium thermodynamic limit such a CG system of RCs behaves like a conventional polydisperse system of nonresponsive particles. We then illustrate the action of external fields, which impose local (position-dependent) property distributions leading to nontrivial effects on the spatial one-body property and density profiles, even for an ideal (noninteracting) gas of RCs. We finally apply density-functional theory in the local density approximation to discuss the effects of particle interactions for specific examples of (i) a suspension of RCs in an external field linear in both position and property, (ii) a suspension of RCs with highly localized properties (sizes) confined between two walls, and (iii) a two-component suspension where an inhomogeneously distributed (nonresponsive) cosolute component, as found, e.g., in the studies of osmolyte- or salt-induced collapse or swelling transitions of thermosensitive polymers, modifies the local properties and density of the RC liquid.

Flow properties reveal the particle-to-polymer transition of ultra-low crosslinked microgels

Exploiting soft, adaptive microgels as building blocks for soft materials with controlled and predictable viscoelastic properties is of great interest for both industry and fundamental research. Here the flow properties of different poly(N-isopropylacrylamide) (pNIPAM) microgels are compared: regularly crosslinked versus ultra-low crosslinked (ULC) microgels. The latter are the softest microgels that can be produced via precipitation polymerization. The viscosity of ULC microgel suspensions at low concentrations can be described with models typically used for hard spheres and regularly crosslinked microgels. In contrast, at higher concentrations, ULC microgels show a much softer behavior compared to regularly crosslinked microgels. The increase of the storage modulus with concentration discloses that while for regularly crosslinked microgels the flow properties are mainly determined by the more crosslinked core, for ULC microgels the brush-like interaction is dominant at high packing fractions. Both the flow curves and the increase of the storage modulus with concentration indicates that ULC microgels can form glass and even reach an apparent jammed state despite their extreme softness. In contrast, the analysis of oscillatory frequency sweep measurements show that when approaching the glass transition the ultra-low crosslinked microgels behave as the regularly crosslinked microgels. This is consistent with a recent study showing that in this concentration range the equilibrium phase behavior of these ULC microgels is the one expected for regularly crosslinked microgels.

Anomalous behavior and structure of a liquid of particles interacting through the harmonic-repulsive pair potential near the crystallization transition

A characteristic property of many soft matter systems is an ultrasoft effective interaction between their structural units. This softness often leads to complex behavior. In particular, ultrasoft systems under pressure demonstrate polymorphism of complex crystal and quasicrystal structures. Therefore, it is of interest to investigate how different can be the structure of the fluid state in such systems at different pressures. Here we address this issue for a model liquid composed of particles interacting through the harmonic-repulsive pair potential. This system can form different crystal structures as the liquid is cooled. We find that, at certain pressures, the liquid exhibits unusual properties, such as a negative thermal expansion coefficient. Besides, the volume and the potential energy of the system can increase during crystallization. At certain pressures, the system demonstrates high stability against crystallization and it is hardly possible to crystallize it on the timescales of the simulations. To address the liquid's structure at high pressures, we consider the scaled pair distribution function (PDF) and the bond-orientational order (BOO) parameters. The marked change happening with the PDF, as pressure increases, is the splitting of the first peak which is caused by the appearance of non-negligible interactions with the second neighbors and the following rearrangement of the structure. Our findings suggest that non-trivial effects, usually explained by different interactions at different spatial scales, can also be observed in one-component systems with simple one-length-scale ultrasoft repulsive interactions.

Monte Carlo simulation and integral equation study of Hertzian spheres in the low-temperature regime

We investigate the behavior of Hertzian spheres in the fluid phase and in proximity of the freezing threshold by using Monte Carlo (MC) simulations and integral equation theories, based on the Ornstein-Zernike (OZ) approach. The study is motivated by the importance of the Hertzian model in representing a large class of systems interacting via soft interactions, such as star polymers or microgels. Radial distribution functions, structure factors, and excess entropy clearly show the reentrant behavior typical of the Hertzian fluid, well caught by both MC simulations and OZ theory. Then, we make use of some phenomenological one-phase criteria for testing their reliability in detecting the freezing threshold. All criteria provide evidence of the fluid-solid transition with different degrees of accuracy. This suggests the possibility to adopt these empirical rules to provide a quick and reasonable estimate of the freezing transition in model potentials of wide interest for soft matter systems.

Nanogels and Microgels: From Model Colloids to Applications, Recent Developments, and Future Trends

Nanogels and microgels are soft, deformable, and penetrable objects with an internal gel-like structure that is swollen by the dispersing solvent. Their softness and the potential to respond to external stimuli like temperature, pressure, pH, ionic strength, and different analytes make them interesting as soft model systems in fundamental research as well as for a broad range of applications, in particular in the field of biological applications. Recent tremendous developments in their synthesis open access to systems with complex architectures and compositions allowing for tailoring microgels with specific properties. At the same time state-of-the-art theoretical and simulation approaches offer deeper understanding of the behavior and structure of nano- and microgels under external influences and confinement at interfaces or at high volume fractions. Developments in the experimental analysis of nano- and microgels have become particularly important for structural investigations covering a broad range of length scales relevant to the internal structure, the overall size and shape, and interparticle interactions in concentrated samples. Here we provide an overview of the state-of-the-art, recent developments as well as emerging trends in the field of nano- and microgels. The following aspects build the focus of our discussion: tailoring (multi)functionality through synthesis; the role in biological and biomedical applications; the structure and properties as a model system, e.g., for densely packed arrangements in bulk and at interfaces; as well as the theory and computer simulation.

Spontaneous deswelling of microgels controlled by counterion clouds

Concentrated poly(N-isopropylacrylamide) (pNIPAM) microgel suspensions at a fixed temperature below the deswelling transition of pNIPAM exhibit spontaneous particle deswelling. The microgels deswell before they are in direct contact and in polydisperse suspensions this deswelling is most pronounced for the largest microgel particles; as a consequence, the polydispersity of the suspension is reduced. Recently, we presented a model for this spontaneous deswelling that is based on the presence of counterions originating from charged groups on the surface of the pNIPAM microgels [A. Scotti et al., Proc. Natl. Acad. Sci. USA 113, 5576 (2016)PNASA60027-842410.1073/pnas.1516011113]. Here we present numerical Poisson-Boltzmann calculations of the electrostatic potential and osmotic pressure inside and outside a pNIPAM microgel that could trigger the observed deswelling at high particle concentrations.

Numerical modelling of non-ionic microgels: an overview

Microgels are complex macromolecules. These colloid-sized polymer networks possess internal degrees of freedom and, depending on the polymer(s) they are made of, can acquire a responsiveness to variations of the environment (temperature, pH, salt concentration, etc.). Besides being valuable for many practical applications, microgels are also extremely important to tackle fundamental physics problems. As a result, these last years have seen a rapid development of protocols for the synthesis of microgels, and more and more research has been devoted to the investigation of their bulk properties. However, from a numerical standpoint the picture is more fragmented, as the inherently multi-scale nature of microgels, whose bulk behaviour crucially depends on the microscopic details, cannot be handled at a single level of coarse-graining. Here we present an overview of the methods and models that have been proposed to describe non-ionic microgels at different length-scales, from the atomistic to the single-particle level. We especially focus on monomer-resolved models, as these have the right level of details to capture the most important properties of microgels, responsiveness and softness. We suggest that these microscopic descriptions, if realistic enough, can be employed as starting points to develop the more coarse-grained representations required to investigate the behaviour of bulk suspensions.

A comparison of the network structure and inner dynamics of homogeneously and heterogeneously crosslinked PNIPAM microgels with high crosslinker content

Poly(N-isopropylacrylamide) microgel particles were prepared via a "classical" surfactant-free precipitation polymerization and a continuous monomer feeding approach. It is anticipated that this yields microgel particles with different internal structures, namely a dense core with a fluffy shell for the classical approach and a more even crosslink distribution in the case of the continuous monomer feeding approach. A thorough structural investigation of the resulting microgels with dynamic light scattering, atomic force microscopy and small angle neutron scattering was conducted and related to neutron spin echo spectroscopy data. In this way a link between structural and dynamic features of the internal polymer network was made.

Linear and nonlinear rheology and structural relaxation in dense glassy and jammed soft repulsive pNIPAM microgel suspensions

We present an integrated experimental and quantitative theoretical study of the mechanics of self-crosslinked, slightly charged, repulsive pNIPAM microgel suspensions over a very wide range of concentrations (c) that span the fluid, glassy and putative "soft jammed" regimes. In the glassy regime we measure a linear elastic dynamic shear modulus over 3 decades which follows an apparent power law concentration dependence G' ∼ c5.64, a variation that appears distinct from prior studies of crosslinked ionic microgel suspensions. At very high concentrations there is a sharp crossover to a nearly linear growth of the modulus. To theoretically understand these observations, we formulate an approach to address all three regimes within a single conceptual Brownian dynamics framework. A minimalist single particle description is constructed that allows microgel size to vary with concentration due to steric de-swelling effects. Using a Hertzian repulsion interparticle potential and a suite of statistical mechanical theories, quantitative predictions under quiescent conditions of microgel collective structure, dynamic localization length, elastic modulus, and the structural relaxation time are made. Based on a constant inter-particle repulsion strength parameter which is determined by requiring the theory to reproduce the linear elastic shear modulus over the entire concentration regime, we demonstrate good agreement between theory and experiment. Testable predictions are then made. We also measured nonlinear rheological properties with a focus on the yield stress and strain. A theoretical analysis with no adjustable parameters predicts how the quiescent structural relaxation time changes under deformation, and how the yield stress and strain change as a function of concentration. Reasonable agreement with our observations is obtained. To the best of our knowledge, this is the first attempt to quantitatively understand structure, quiescent relaxation and shear elasticity, and nonlinear yielding of dense microgel suspensions using microscopic force based theoretical methods that include activated hopping processes. We expect our approach will be useful for other soft polymeric particle suspensions in the core-shell family.

Unusual polysaccharide rheology of aqueous dispersions of soft phytoglycogen nanoparticles

Phytoglycogen is a natural polysaccharide produced in the form of dense, 35 nm diameter nanoparticles by some varieties of plants such as sweet corn. The highly-branched, dendrimeric structure of phytoglycogen leads to interesting and useful properties such as softness and deformability of the particles, and a strong interaction with water. These properties make the particles ideal for use as unique additives in personal care, nutrition and biomedical formulations. In the present study, we describe rheology measurements of aqueous dispersions of phytoglycogen nanoparticles. The viscosity of the dispersions remained Newtonian up to large concentrations (∼20% w/w). For higher concentrations, the zero-shear viscosity increased dramatically, reaching values that exceeded that of the water solvent by six orders of magnitude at a concentration of 30% w/w and were well described by the Vogel-Fulcher-Tammann relation of glassy dynamics. The very large values of the zero-shear viscosity are coupled with significant deformation of the soft nanoparticles. We quantified the softness of the particles by performing osmotic pressure measurements on concentrated dispersions, obtaining a value of 15 kPa for the compressional modulus. For the most concentrated samples, we observed flow at stresses less than the apparent yield stress value determined by fitting the high strain rate data to the Herschel-Bulkley model. This behavior, similar to that of star polymer glasses, suggests the possibility of a hairy colloid particle geometry. Remarkably, phytoglycogen nanoparticles dispersed in water provide a very simple experimental realization of glass-forming dispersions of soft colloidal particles that can be used to validate theoretical models in detail.

Quenching of fully symmetric mixtures of oppositely charged microgels: the role of soft stiffness

Using molecular dynamics simulations, we investigate the self-assembly of a coarse-grained binary system of oppositely charged microgels, symmetric in size and concentration. The microgel pair interactions are described by an effective pair potential which implicitly accounts for the averaged ionic contributions, in addition to a short-range elastic repulsion that accounts for the overlapping of the polymer chains, the latter being described by the Hertzian interaction. Particular emphasis is placed on the role played by the strength of the soft repulsive interaction on the resulting particle aggregation. It is found that the possibility of particle inter-penetration in oppositely charged soft particles results in a much wider variety of cluster morphologies in comparison with their hard-spheres counterparts. Specifically, the softness of the steric interactions enhances the competition between repulsive and attractive electrostatic interactions, leading to the formation of aggregates that are comprised of strongly bounded charged particles displaying a low degree of charge ordering.

Ionic size effects on the Poisson-Boltzmann theory

In this paper, we develop a simple theory to study the effects of ionic size on ionic distributions around a charged spherical particle. We include a correction to the regular Poisson-Boltzmann equation in order to take into account the size of ions in a mean-field regime. The results are compared with Monte Carlo simulations and a density functional theory based on the fundamental measure approach and a second-order bulk expansion which accounts for electrostatic correlations. The agreement is very good even for multivalent ions. Our results show that the theory can be applied with very good accuracy in the description of ions with highly effective ionic radii and low concentration, interacting with a colloid or a nanoparticle in an electrolyte solution.

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.

Self-Assembly of Ionic Microgels Driven by an Alternating Electric Field: Theory, Simulations, and Experiments

The structural properties of a system of ionic microgels under the influence of an alternating electric field are investigated both theoretically and experimentally. This combined investigation aims to shed light on the structural transitions that can be induced by changing either the driving frequency or the strength of the applied field, which range from string-like formation along the field to crystal-like structures across the orthogonal plane. In order to highlight the physical mechanisms responsible for the observed particle self-assembly, we develop a coarse-grained description, in which effective interactions among the charged microgels are induced by both equilibrium ionic distributions and their time-averaged hydrodynamic responses to the applied field. These contributions are modeled by the buildup of an effective dipole moment at the microgels backbones, which is partially screened by their ionic double layer. We show that this description is able to capture the structural properties of this system, allowing for very good agreement with the experimental results. The model coarse-graining parameters are indirectly obtained via the measured pair distribution functions and then further assigned with a clear physical interpretation, allowing us to highlight the main physical mechanisms accounting for the observed self-assembly behavior.

A phenomenological order approach to the volume phase transition in microgel particles

A phenomenological insight into the volume transition of microgel particles via nematic-like ordering behavior through solvation processes of a polymer matrix.