Search

We investigate the link between the internal microstructure of poly( N -isopropylacrylamide)-poly(ethylene glycol) methyl ether methacrylate (PNIPAM-PEGMA) microgels, their bulk moduli and the rheological response and structural arrangement in dense suspensions. The low degree of crosslinking combined with the increased hydrophilicity induced by the presence of PEGMA results in a diffuse, star-like density profile of the particle and very low values of the bulk modulus in dilute conditions, as determined by small angle neutron scattering (SANS). The ultrasoft nature of the particle is reflected in the changes of the structural arrangement in dense suspensions, which evidence a strong deswelling and a sharp rise of the bulk modulus at moderate packing fractions. At larger packings the single particle morphology and softness saturate, and we observe a structural transition from a dispersion-like to a hydrogel-like behavior. The transition is also reflected in the rheological response in the form of a two-step yielding at large packing fractions, characteristic of systems in which a network structure is present. Our results demonstrate that a knowledge of the internal structure and mechanics of individual microgels is needed to determine and tune the properties of dense suspensions, and optimize their response for applications in biomedicine and as filtration systems.

Deformable colloids and macromolecules adsorb at interfaces as they decrease the interfacial energy between the two media. The deformability, or softness, of these particles plays a pivotal role in the properties of the interface. In this study, we employ a comprehensive in situ approach, combining neutron reflectometry with molecular dynamics simulations, to thoroughly examine the profound influence of softness on the structure of microgel Langmuir monolayers under compression. Lateral compression of both hard and soft microgel particle monolayers induces substantial structural alterations, leading to an amplified protrusion of the microgels into the aqueous phase. However, a critical distinction emerges: hard microgels are pushed away from the interface, in stark contrast to the soft ones, which remain firmly anchored to it. Concurrently, on the air-exposed side of the monolayer, lateral compression induces a flattening of the surface of the hard monolayer. This phenomenon is not observed for the soft particles as the monolayer is already extremely flat even in the absence of compression. These findings significantly advance our understanding of the key role of softness on both the equilibrium phase behavior of the monolayer and its effect when soft colloids are used as stabilizers of responsive interfaces and emulsions.

We investigate by means of small angle neutron scattering experiments and numerical simulations the interactions and inter-particle arrangements of concentrated dispersions of copolymer poly( N -isopropylacrylamide)-poly(ethylene glycol methyl ether methacrylate) (PNIPAM-PEGMA) microgels across the volume phase transition (VPT). The scattering data of moderately concentrated dispersions are accurately modeled at all temperatures by using a star polymer form factor and static structure factors calculated from the effective potential obtained from simulations. Interestingly, for temperatures below the VPT temperature (VPTT), the radius of gyration and blob size of the particles significantly decrease with increasing the effective packing fraction in the non-overlapping regime. This is attributed to the presence of charges in the system associated with the use of an ionic initiator in the synthesis. Simulations using the experimentally corroborated interaction potential are used to explore the state diagram in a wide range of effective packing fractions. Below and slightly above the VPTT, the system undergoes an arrest transition mainly driven by the soft repulsion between the particles. Only well above the VPTT the system is found to phase separate before arresting. Our results highlight the versatility and potential of copolymer PNIPAM-PEGMA microgels to explore different kinds of arrested states balancing attraction and repulsion by changing temperature and packing fraction.

We perform large-scale computer simulations of an off-lattice two-dimensional model of active particles undergoing a motility-induced phase separation (MIPS) to investigate the system's critical behaviour close to the critical point of the MIPS curve. By sampling steady-state configurations for large system sizes and performing finite size scaling analysis we provide exhaustive evidence that the critical behaviour of this active system belongs to the Ising universality class. In addition to the scaling observables that are also typical of passive systems, we study the critical behaviour of the kinetic temperature difference between the two active phases. This quantity, which is always zero in equilibrium, displays instead a critical behavior in the active system which is well described by the same exponent of the order parameter in agreement with mean-field theory.

Recent progress has been made in the numerical modelling of neutral microgel particles with a realistic, disordered structure. In this work we extend this approach to the case of co-polymerised microgels where a thermoresponsive polymer is mixed with acidic groups. We compare the cases where counterions directly interact with microgel charges or are modelled implicitly through a Debye-Hückel description. We do so by performing extensive numerical simulations of single microgels across the volume phase transition (VPT) varying the temperature and the fraction of charged monomers. We find that the presence of charges considerably alters the microgel structure, quantified by the monomer density profiles and by the form factors of the microgels, particularly close to the VPT. We observe significant deviations between the implicit and explicit models, with the latter comparing more favourably to available experiments. In particular, we observe a shift of the VPT temperature to larger values as the amount of charged monomers increases. We also find that below the VPT the microgel-counterion complex is almost neutral, while it develops a net charge above the VPT. Interestingly, under these conditions the collapsed microgel still retains a large amount of counterions inside its structure. Since these interesting features cannot be captured by the implicit model, our results show that it is crucial to explicitly include the counterions in order to realistically model ionic thermoresponsive microgels.

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.

Effective colloid-colloid interactions can be tailored through the addition of a complex cosolute. Here we investigate the case of a cosolute made by self-assembling patchy particles. Depending on the valence, these particles can form either polymer chains or branched structures. We numerically calculate the effective potential V eff between two colloids immersed in a suspension of reversible patchy particles, exploring a wide region of the cosolute phase diagram and the role of valence. In addition to well-known excluded volume and depletion effects, we find that, under appropriate conditions, V eff is completely attractive but shows an oscillatory character. In the case of polymerizing cosolute, this results from the fact that chains are efficiently confined by the colloids through the onset of local order. This argument is then generalized to the case of particles with higher valence, under the condition that they are still able to maintain a fully bonded organization upon confinement. The resulting effective potentials are relevant for understanding the behavior of complex mixtures in crowded environments, but may also be exploited for tuning colloidal self-assembly at preferred target distances in order to build desired superstructures.

We report numerical results for the effective potential arising between two colloids immersed in a self-assembling cosolute which forms reversible clusters. The potential is evaluated at cosolute state points with different densities and temperatures but with the same connectivity properties. We find that the range of the resulting effective potential is controlled only by the cosolute thermal correlation length rather than by its connectivity length. We discuss the significant differences from previous results focusing on the cosolute forming irreversible clusters and we show that the irreversible bond case represents a singular limit which cannot be accessed in equilibrium by continuously increasing the bond lifetime.

We investigate binary mixtures of large colloids interacting through soft potentials with small, ideal depletants. We show that softness has a dramatic effect on the resulting colloid-colloid effective potential when the depletant-to-colloid size ratio q is small, with significant consequences on the colloidal phase behaviour. We provide an exact relationship that allows us to obtain the effective pair potential for any type of colloid-depletant interaction in the case of ideal depletants, without having to rely on complicated and expensive full-mixture simulations. We also show that soft repulsion among depletants further enhances the tendency of colloids to aggregate. Our theoretical and numerical results demonstrate that--in the limit of small q--soft mixtures cannot be mapped onto hard systems and hence soft depletion is not a mere extension of the widely used Asakura-Oosawa potential.

The aging dynamics of a colloidal glass has been studied by multiangle dynamic light scattering, neutron spin echo, X-ray photon correlation spectroscopy and molecular dynamics simulations. The two relaxation processes, microscopic (fast) and structural (slow), have been investigated in an unprecedentedly wide range of time and length scales covering both ergodic and nonergodic regimes. The microscopic relaxation time remains diffusive at all length scales across the glass transition scaling with wavevector Q as Q(-2). The length-scale dependence of structural relaxation time changes from diffusive, characterized by a Q(-2)-dependence in the early stages of aging, to a Q(-1)-dependence in the full aging regime which marks a discontinuous hopping dynamics. Both regimes are associated with a stretched behaviour of the correlation functions. We expect these findings to provide a general description of both relaxations across the glass transition.

We investigate the dynamics of polydisperse hard spheres at high packing fractions ϕ. We use extensive numerical simulations based on an experimentally-realistic particle size distribution (PSD) and compare to commonly-used PSDs such as Gaussian or top hat distribution. We find that the mode of kinetic arrest depends on the PSD's shape and not only on its variance. For the experimentally-realistic PSD we find ageing dynamics even though the density correlators decay fully to zero for ϕ ≥ 0.59. We observe substantial decoupling of the dynamics of the smallest and largest particles. While the smallest particles remain diffusive in all our simulations, a power-law describes the largest-particle diffusion, suggesting an ideal arrest at ϕc ∼ 0.588. The latter is however averted just before ϕc, due to the presence of the mobile smallest particles. In addition, we identify that a partial aging mechanism is at work, whose effects are most pronounced for the largest particles. By comparing our results with recent experimental observations of ergodic behavior up to ϕ ∼ 0.6 in a hard-sphere system, we argue that this is an effect of polydispersity, which smears out the glass transition.

We report event-driven molecular dynamics simulations of the irreversible gelation of hard ellipsoids of revolution containing several associating groups, characterizing how the cluster size distribution evolves as a function of the extent of reaction, both below and above the gel point. We find that over a very large interval of values of the extent of reaction, parameter-free mean-field predictions are extremely accurate, providing evidence that in this model the Ginzburg zone near the gel point, where non-mean field effects are important, is very limited. We also find that the Flory's hypothesis for the post-gelation regime properly describes the connectivity of the clusters even if the long-time limit of the extent of reaction does not reach the fully reacted state. This study shows that irreversibly aggregating asymmetric hard-core patchy particles may provide a close realization of the mean-field model, for which available theoretical predictions may help control the structure and the connectivity of the gel state. Besides chemical gels, the model is relevant to network-forming soft materials like systems with bioselective interactions, functionalized molecules and patchy colloids.