In Silico Synthesis of Microgel Particles

Revisiting the density profile of the fuzzy sphere model for microgel colloids

Common neutral polymer microgels exhibit an inhomogeneous density profile with a gradual decay that is commonly described using the fuzzy sphere model. The model is based on the idea of convolving the collapsed solid sphere profile with a Gaussian to describe inhomogeneous swelling of the microgel in a good solvent. Here we show that the corresponding density profile in real space used in several recent works - such as in super-resolution microscopy - is different from the fuzzy sphere model, and we explain how to correctly transition between reciprocal space modelling to real space. Our work aims to clarify the application of the model so that errors can be avoided in the future. Our discussion is also crucial when comparing alternative real-space models for the density profile with the established fuzzy sphere model.

Core-shell nanogels: the effects of morphology, electro- and magnetostatic interactions

We study the influence of core-shell morphology on the structural characteristics of nanogels. Using computer simulations, we examine three different types of systems, distinguished by their intermonomer interactions: those with excluded volume only; those with charged monomers and excluded volume; and those with excluded volume combined with a certain number of magnetised nanoparticles incorporated within the nanogel. We observe that if the polymers in the shell are short and dense, they tend to penetrate the core. This effect of backfolding is enhanced in charged nanogels, regardless of whether all monomers are charged, or only the core or shell ones. The presence of an experimentally available amount of magnetic nanoparticles in a gel, on the one hand, does not lead to any significant morphological changes. On the other hand, the morphology of the nanogel with magnetic particles has an impact on its magnetic susceptibility. Particular growth of the magnetic response is observed if a long shell of a nanogel is functionalised.

Topological Programmability of Isomerizable Polymers

Topology isomerizable networks (TINs) can be programmed into numerous polymers exhibiting unique and spatially defined (thermo-) mechanical properties. However, capturing the dynamics in topological transformations and revealing the intrinsic mechanisms of mechanical property modulation at the microscopic level is a significant challenge. Here, we use a combination of coarse-grained molecular dynamics simulations and reaction kinetic theory to reveal the impact of dynamic bond exchange reactions on the topology of branched chains. We find that, the grafted units follow a geometric distribution with a converged uniformity, which depends solely on the average grafted units of branched chains. Furthermore, we demonstrate that the topological structure can lead to spontaneous modulation of mechanical properties. The theoretical framework provides a research paradigm for studying the topology and mechanical properties of TINs.

Real and In Silico Microgels Show Comparable Bulk Moduli Below and Above the Volume Phase Transition

The compressibility of soft colloids influences their phase behavior and flow properties, especially in concentrated suspensions. Particle compressibility, which is proportional to the reciprocal of the bulk modulus K, is a key parameter for soft polymer-based particles that can be compressed in crowded environments. Here, microgels with different degrees of cross-linking, i.e., softness, are investigated below and above their volume phase transition temperature (VPTT). By combining molecular dynamics simulations with small-angle neutron scattering with contrast variation, a change in the particle bulk moduli of two orders of magnitude is observed. The degree of cross-linking has a significant impact on the bulk modulus of the swollen microgel, while above the VPTT the values of K are almost independent of the cross-linking density. The excellent agreement between experimental results and simulations also highlight that the model microgels from computer simulations possess both the internal architecture and the elastic properties of real polymeric networks. This paves the way to a systematic use of simulations to investigate the behavior of dense microgel suspensions below and above their VPTT.

Softness matters: effects of compression on the behavior of adsorbed microgels at interfaces

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.

Auxetic polymer networks: The role of crosslinking, density, and disorder

Low-crosslinked polymer networks have recently been found to behave auxetically when subjected to small tensions, that is, their Poisson's ratio ν becomes negative. In addition, for specific state points, numerical simulations revealed that diamond-like networks reach the limit of mechanical stability, exhibiting values of ν = -1, a condition that we define as hyper-auxeticity. This behavior is interesting per se for its consequences in materials science but is also appealing for fundamental physics because the mechanical instability is accompanied by evidence of criticality. In this work, we deepen our understanding of this phenomenon by performing a large set of equilibrium and stress-strain simulations in combination with phenomenological elasticity theory. The two approaches are found to be in good agreement, confirming the above results. We also extend our investigations to disordered polymer networks and find that the hyper-auxetic behavior also holds in this case, still manifesting a similar critical-like behavior as in the diamond one. Finally, we highlight the role of the number density, which is found to be a relevant control parameter determining the elastic properties of the system. The validity of the results under disordered conditions paves the way for an experimental investigation of this phenomenon in real systems, such as hydrogels.

Toward a Unified Description of the Electrostatic Assembly of Microgels and Nanoparticles

The interplay of soft responsive particles, such as microgels, with nanoparticles (NPs) yields highly versatile complexes that show great potential for applications, ranging from plasmonic sensing to catalysis and drug delivery. However, the microgel-NP assembly process has not been investigated so far at the microscopic level, thus hindering the possibility of designing such hybrid systems a priori. In this work, we combine state-of-the-art numerical simulations with experiments to elucidate the fundamental mechanisms taking place when microgel-NP assembly is controlled by electrostatic interactions and the associated effects on the structure of the resulting complexes. We find a general behavior where, by increasing the number of interacting NPs, the microgel deswells up to a minimum size after which a plateau behavior occurs. This occurs either when NPs are mainly adsorbed to the microgel corona via the folding of the more external chains or when NPs penetrate inside the microgel, thereby inducing a collective reorganization of the polymer network. By varying microgel properties, such as fraction of cross-linkers or charge, as well as NP size and charge, we further show that the microgel deswelling curves can be rescaled onto a single master curve, for both experiments and simulations, demonstrating that the process is entirely controlled by the charge of the whole microgel-NP complex. Our results thus have a direct relevance in fundamental materials science and offer novel tools to tailor the nanofabrication of hybrid devices of technological interest.

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.

Interactions between interfaces dictate stimuli-responsive emulsion behaviour

Stimuli-responsive emulsions offer a dual advantage, combining long-term storage with controlled release triggered by external cues such as pH or temperature changes. This study establishes that thermo-responsive emulsion behaviour is primarily determined by interactions between, rather than within, interfaces. Consequently, the stability of these emulsions is intricately tied to the nature of the stabilizing microgel particles - whether they are more polymeric or colloidal, and the morphology they assume at the liquid interface. The colloidal properties of the microgels provide the foundation for the long-term stability of Pickering emulsions. However, limited deformability can lead to non-responsive emulsions. Conversely, the polymeric properties of the microgels enable them to spread and flatten at the liquid interface, enabling stimuli-responsive behaviour. Furthermore, microgels shared between two emulsion droplets in flocculated emulsions facilitate stimuli-responsiveness, regardless of their internal architecture. This underscores the pivotal role of microgel morphology and the forces they exert on liquid interfaces in the control and design of stimuli-responsive emulsions and interfaces.

Concentration and temperature dependent interactions and state diagram of dispersions of copolymer microgels

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.

Buckling and Interfacial Deformation of Fluorescent Poly(N-isopropylacrylamide) Microgel Capsules

Hollow microgels are fascinating model systems at the crossover between polymer vesicles, emulsions, and colloids as they deform, interpenetrate, and eventually shrink at higher volume fraction or when subjected to an external stress. Here, we introduce a system consisting of microgels with a micrometer-sized cavity enabling a straightforward characterization in situ using fluorescence microscopy techniques. Similarly to elastic capsules, these systems are found to reversibly buckle above a critical osmotic pressure, conversely to smaller hollow microgels, which were previously reported to deswell at high volume fraction. Simulations performed on monomer-resolved in silico hollow microgels confirm the buckling transition and show that the presented microgels can be described with a thin shell model theory. When brought to an interface, these microgels, that we define as microgel capsules, strongly deform and we thus propose to utilize them to locally probe interfacial properties within a theoretical framework adapted from the Johnson-Kendall-Roberts (JKR) theory. Besides their capability to sense their environment and to address fundamental questions on the elasticity and permeability of microgel systems, microgel capsules can be further envisioned as model systems mimicking anisotropic responsive biological systems such as red blood and epithelial cells thanks to the possibility offered by microgels to be synthesized with custom-designed properties.

Structure of swollen hollow polyelectrolyte nanogels with inhomogeneous cross-link distribution

HYPOTHESIS

Recently, it has become possible to synthesize hollow polyelectrolyte nano- and microgels. The shell permeability can be controlled by external stimuli, while the cavity can serve as a storage place for guest molecules. However, there is a lack of a detailed understanding at the molecular level regarding the role of the network topology, inhomogeneities of the distribution of cross-links, and the impact of the electrostatics on the structural response of hollow microgel to external stimuli. To bridge these gaps, molecular dynamics (MD) of computer simulations are used.

EXPERIMENTS

Here, we propose a fresh methodology to create realistic hollow microgel particles in silico. The technique involves a gradual change in the average local length of subchains depending on the distance to the center of mass of the microgel particles resulting in the microgels with a non-uniform distribution of cross-linking species. In particular, a series of microgels with (i) a highly cross-linked inner part of the shell and gradually decreased cross-linker concentration towards the periphery, (ii) microgels with loosely cross-linked inner and outer parts, as well as (iii) microgels with a more-or-less homogeneous structure, have been created and validated. Counterions and salt ions are taken into account explicitly, and electrostatic interactions are described by the Coulomb potential.

FINDINGS

Our studies reveal a strong dependence of the microgel swelling response on the network topology. Simple redistribution of cross-links plays a significant role in the structure of the microgels, including cavity size, microgel size, fuzziness, and extension of the inner and outer parts of the microgels. Our results indicate the possibilities of qualitative justification of the structure of the hollow microgels in the experiments by measuring the relative change in the size of the sacrificial core to the size of the cavity and by estimation of a power law function, [Formula: see text] , of the hydrodynamic radius of the hollow microgels as a function of added salt concentration.

Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks

The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.

Membrane-Mediated Interactions Between Nonspherical Elastic Particles

The transport of particles across lipid-bilayer membranes is important for biological cells to exchange information and material with their environment. Large particles often get wrapped by membranes, a process which has been intensively investigated in the case of hard particles. However, many particles in vivo and in vitro are deformable, e.g., vesicles, filamentous viruses, macromolecular condensates, polymer-grafted nanoparticles, and microgels. Vesicles may serve as a generic model system for deformable particles. Here, we study nonspherical vesicles with various sizes, shapes, and elastic properties at initially planar lipid-bilayer membranes. Using the Helfrich Hamiltonian, triangulated membranes, and energy minimization, we predict the interplay of vesicle shapes and wrapping states. Increasing particle softness enhances the stability of shallow-wrapped and deep-wrapped states over nonwrapped and complete-wrapped states. The free membrane mediates an interaction between partial-wrapped vesicles. For the pair interaction between deep-wrapped vesicles, we predict repulsion. For shallow-wrapped vesicles, we predict attraction for tip-to-tip orientation and repulsion for side-by-side orientation. Our predictions may guide the design and fabrication of deformable particles for efficient use in medical applications, such as targeted drug delivery.

Probing Temperature Responsivity of Microgels and Its Interplay with a Solid Surface by Super-Resolution Microscopy and Numerical Simulations

Super-resolution microscopy has become a powerful tool to investigate the internal structure of complex colloidal and polymeric systems, such as microgels, at the nanometer scale. An interesting feature of this method is the possibility of monitoring microgel response to temperature changes in situ. However, when performing advanced microscopy experiments, interactions between the particle and the environment can be important. Often microgels are deposited on a substrate, since they have to remain still for several minutes during the experiment. This study uses direct stochastic optical reconstruction microscopy (dSTORM) and advanced coarse-grained molecular dynamics simulations to investigate how individual microgels anchored on hydrophilic and hydrophobic surfaces undergo their volume phase transition with temperature. We find that, in the presence of a hydrophilic substrate, the structure of the microgel is unperturbed and the resulting density profiles quantitatively agree with simulations performed under bulk conditions. Instead, when a hydrophobic surface is used, the microgel spreads at the interface and an interesting competition between the two hydrophobic strengths,monomer-monomer vs monomer-surface,comes into play at high temperatures. The robust agreement between experiments and simulations makes the present study a fundamental step to establish this high-resolution monitoring technique as a platform for investigating more complex systems, these being either macromolecules with peculiar internal structure or nanocomplexes where molecules of interest can be encapsulated in the microgel network and controllably released with temperature.

Ionisation and swelling behaviour of weak polyampholyte core-shell networks - a Monte Carlo study

The network charge of polyampholyte microgels can be tuned by varying the pH of the surrounding solution, and a charge reversal from a positively charged microgel at low pH to a negatively charged microgel at high pH can be achieved. In a titration experiment, it is difficult to tell apart the ionisation of the acidic and basic monomers in the network and to determine the distribution of charges in the network, whereas using Metropolis Monte Carlo simulations, both the degree of ionisation and the distribution of ionised monomers can be determined separately for both species. Building on our earlier work on alternating polyampholyte microgels, we now investigated the pH-dependent ionisation and the swelling behaviour of polyampholyte core-shell microgels under good solvent conditions. For this purpose, we performed Metropolis Monte Carlo simulations for a bead-spring model using the constant-pH method. As in our previous study on alternating microgels, the width of the U-shaped curve of the microgels volume as a function of pH depends on the relative dissociation constants of acid and base, and the microgel volume can be approximated by a linear function of the total network charge. Due to the spatial separation of acid and base in core-shell systems, the ionisation is less enhanced compared to a microgel with an alternating distribution of the two species. Nevertheless, we still see an influence of the presence of one species on the ionisation behaviour of the other species under good solvent conditions. Furthermore, the isoelectric point is shifted towards higher pH, which is caused by a higher charge density in the core compared to that in the shell. Added salt changes the Donnan equilibrium, which determines the counterion distribution within and outside of the microgel. At the same time, it contributes to the electrostatic screening of the network charges, leading to a narrowing of the U-shaped volume transition curve.

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.

Molecular dynamics study of the swelling and osmotic properties of compact nanogel particles

Owing to their great importance in materials science and other fields, we investigate the solution and osmotic properties of uncharged compact nanogel particles over a wide range of solvent quality and particle concentration by molecular dynamics (MD) simulations. We characterize the osmotic pressure by estimating the second and third virial coefficients, and by extension, we identify the θ-point where the second virial coefficient vanishes. Calculations of the structure factor indicate that these particles are similar to macrogels in that the particle-like scattering profile disappears at moderate concentrations. We also find that improving the solvent quality enhances the spatial segmental uniformity, while significant heterogeneous structure arises near the θ-point. Well below the θ-point where the second osmotic virial coefficient vanishes, these heterogeneous structures become less prevalent as the particles tend to collapse. We also investigate the degree of swelling and structure of compact nanogel particles with a variable excluded volume interaction and gel particle concentration. The osmotic modulus and the scaling exponents in good and θ-point conditions of these gels are characteristic of interacting randomly branched polymers, i.e., "lattice animals".

Designing Enhanced Entropy Binding in Single-Chain Nanoparticles

Single-chain nanoparticles (SCNPs) are a new class of bio- and soft-matter polymeric objects in which a fraction of the monomers are able to form equivalently intra- or interpolymer bonds. Here we numerically show that a fully entropic gas-liquid phase separation can take place in SCNP systems. Control over the discontinuous (first-order) change-from a phase of independent diluted (fully-bonded) polymers to a phase in which polymers entropically bind to each other to form a (fully-bonded) polymer network-can be achieved by a judicious design of the patterns of reactive monomers along the polymer chain. Such a sensitivity arises from a delicate balance between the distinct entropic contributions controlling the binding.

In-situ study of the impact of temperature and architecture on the interfacial structure of microgels

The structural characterization of microgels at interfaces is fundamental to understand both their 2D phase behavior and their role as stabilizers that enable emulsions to be broken on demand. However, this characterization is usually limited by available experimental techniques, which do not allow a direct investigation at interfaces. To overcome this difficulty, here we employ neutron reflectometry, which allows us to probe the structure and responsiveness of the microgels in-situ at the air-water interface. We investigate two types of microgels with different cross-link density, thus having different softness and deformability, both below and above their volume phase transition temperature, by combining experiments with computer simulations of in silico synthesized microgels. We find that temperature only affects the portion of microgels in water, while the strongest effect of the microgels softness is observed in their ability to protrude into the air. In particular, standard microgels have an apparent contact angle of few degrees, while ultra-low cross-linked microgels form a flat polymeric layer with zero contact angle. Altogether, this study provides an in-depth microscopic description of how different microgel architectures affect their arrangements at interfaces, and will be the foundation for a better understanding of their phase behavior and assembly.

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.

Effect of network topology and crosslinker reactivity on microgel structure and ordering at liquid-liquid interface

Polymer microgels synthesized in silico were studied at a liquid-liquid interface via mesoscopic computer simulations and compared to microgels with ideal (diamond-like) structure. The effect of crosslinkers reactivity ratio on the single particle morphology at the interface and monolayer behavior was examined. It was demonstrated that single particles deform into an explicit core-corona morphology when adsorbed at the interface. An increase in the crosslinker reactivity ratio decreased both the deformation ratio and the ratio between the core and corona sizes. Meanwhile, the compression of microgel monolayers revealed the existence of five distinct interparticle contact regimes, which have been observed experimentally in the literature. The crosslinker reactivity ratio appeared to define the compression range in these regimes and the sharpness of the transition between them. In particular, the higher the crosslinker reactivity ratio, the smaller the corona, and in turn, the narrower the range of the intermediate regime comprising both core-core and corona-corona contacts. The obtained results demonstrate that the more realistic model of microgels synthesized via precipitation polymerization allows for a more accurate prediction of the properties of the microgels at a liquid-liquid interface in comparison to the conventional diamond-like lattice model.

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.

Link between Morphology, Structure, and Interactions of Composite Microgels

We combine small-angle scattering experiments and simulations to investigate the internal structure and interactions of composite poly(N-isopropylacrylamide)-poly(ethylene glycol) (PNIPAM-PEG) microgels. At low temperatures the experimentally determined form factors and the simulated density profiles indicate a loose internal particle structure with an extended corona that can be modeled as a starlike object. With increasing temperature across the volumetric phase transition, the form factor develops an inflection that, using simulations, is interpreted as arising from a conformation in which PEG chains are incorporated in the interior of the PNIPAM network. This gives rise to a peculiar density profile characterized by two dense, separated regions, at odds with configurations in which the PEG chains reside on the surface of the PNIPAM core. The conformation of the PEG chains also have profound effects on the interparticle interactions: Although chains on the surface reduce the solvophobic attraction typically experienced by PNIPAM particles at high temperatures, PEG chains inside the PNIPAM network shift the onset of attractive interaction at even lower temperatures. Our results show that by tuning the morphology of the composite microgels, we can qualitatively change both their structure and their mutual interactions, opening the way to explore new collective behaviors of these objects.

Onset of criticality in hyper-auxetic polymer networks

Against common sense, auxetic materials expand or contract perpendicularly when stretched or compressed, respectively, by uniaxial strain, being characterized by a negative Poisson's ratio ν. The amount of deformation in response to the applied force can be at most equal to the imposed one, so that ν = - 1 is the lowest bound for the mechanical stability of solids, a condition here defined as "hyper-auxeticity". In this work, we numerically show that ultra-low-crosslinked polymer networks under tension display hyper-auxetic behavior at a finite crosslinker concentration. At this point, the nearby mechanical instability triggers the onset of a critical-like transition between two states of different densities. This phenomenon displays similar features as well as important differences with respect to gas-liquid phase separation. Since our model is able to faithfully describe real-world hydrogels, the present results can be readily tested in laboratory experiments, paving the way to explore this unconventional phase behavior.

On the role of softness in ionic microgel interactions

Thermoresponsive microgels are a popular model system to study phase transitions in soft matter, because temperature directly controls their volume fraction. Ionic microgels are additionally pH-responsive and possess a rich phase diagram. Although effective interaction potentials between microgel particles have been proposed, these have never been fully tested, leading to a gap in our understanding of the link between single-particle and collective properties. To help resolve this gap, four sets of ionic microgels with varying crosslinker density were synthesised and characterised using light scattering techniques and confocal microscopy. The resultant structural and dynamical information was used to investigate how particle softness affects the phase behaviour of ionic microgels and to validate the proposed interaction potential. We find that the architecture of the microgel plays a marked role in its phase behaviour. Rather than the ionic charges, it is the dangling ends which drive phase transitions and interactions at low concentration. Comparison to theory underlines the need for a refined theoretical model which takes into consideration these close-contact interactions.

Structure and conformational properties of ideal nanogel particles in athermal solutions

We investigate the conformational properties of "ideal" nanogel particles having a lattice network topology by molecular dynamics simulations to quantify the influence of polymer topology on the solution properties of this type of branched molecular architecture. In particular, we calculate the mass scaling of the radius of gyration (Rg), the hydrodynamic radius, as well as the intrinsic viscosity with the variation of the degree of branching, the length of the chains between the branched points, and the average mesh size within these nanogel particles under good solvent conditions. We find competing trends between the molecular characteristics, where an increase in mesh size or degree of branching results in the emergence of particle-like characteristics, while an increase in the chain length enhances linear polymer-like characteristics. This crossover between these limiting behaviors is also apparent in our calculation of the form factor, P(q), for these structures. Specifically, a primary scattering peak emerges, characterizing the overall nanogel particle size. Moreover, a distinct power-law regime emerges in P(q) at length scales larger than the chain size but smaller than Rg of the nanogel particle, and the Rg mass scaling exponent progressively approaches zero as the mesh size increases, the same scaling as for an infinite network of Gaussian chains. The "fuzzy sphere" model does not capture this feature, and we propose an extension to this popular model. These structural features become more pronounced for values of molecular parameters that enhance the localization of the branching segments within the nanogel particle.

Sharing the Load: Stress Redistribution Governs Fracture of Polymer Double Networks

The stress response of polymer double networks depends not only on the properties of the constituent networks but also on the interactions arising between them. Here, we demonstrate, via coarse-grained simulations, that both their global stress response and their microscopic fracture mechanics are governed by load sharing through these internetwork interactions. By comparing our results with affine predictions, where stress redistribution is by definition homogeneous, we show that stress redistribution is highly inhomogeneous. In particular, the affine prediction overestimates the fraction of broken chains by almost an order of magnitude. Furthermore, homogeneous stress distribution predicts a single fracture process, while in our simulations, fracture of sacrificial chains takes place in two steps governed by load sharing within a network and between networks, respectively. Our results thus provide a detailed microscopic picture of how inhomogeneous stress redistribution after rupture of chains governs the fracture of polymer double networks.

Two-step deswelling in the Volume Phase Transition of thermoresponsive microgels

Microgels, colloidal-scale polymer networks, are the prototype soft colloids. When the constituent polymers are thermoresponsive, they undergo a volume phase transition (VPT) from a swollen to a collapsed state at a characteristic temperature, close to ambient one, of great appeal for several applications. To describe this phenomenon, microgels are usually treated as neutral, but here we show that electrostatics needs to be taken into account. In particular, deswelling occurs via a two-step, rather than a homogeneous, particle collapse, mainly driven by peripheral charges located on the microgel corona, for which we also establish a unifying framework encompassing all studied microgels. Our work thus provides a change of perspective to describe these fascinating systems.

Thermoresponsive microgels are one of the most investigated types of soft colloids, thanks to their ability to undergo a Volume Phase Transition (VPT) close to ambient temperature. However, this fundamental phenomenon still lacks a detailed microscopic understanding, particularly regarding the presence and the role of charges in the deswelling process. This is particularly important for the widely used poly(N-isopropylacrylamide)–based microgels, where the constituent monomers are neutral but charged groups arise due to the initiator molecules used in the synthesis. Here, we address this point combining experiments with state-of-the-art simulations to show that the microgel collapse does not happen in a homogeneous fashion, but through a two-step mechanism, entirely attributable to electrostatic effects. The signature of this phenomenon is the emergence of a minimum in the ratio between gyration and hydrodynamic radii at the VPT. Thanks to simulations of microgels with different cross-linker concentrations, charge contents, and charge distributions, we provide evidence that peripheral charges arising from the synthesis are responsible for this behavior and we further build a universal master curve able to predict the two-step deswelling. Our results have direct relevance on fundamental soft condensed matter science and on applications where microgels are involved, ranging from materials to biomedical technologies.

Effect of Internal Architecture on the Assembly of Soft Particles at Fluid Interfaces

Monolayers of soft colloidal particles confined at fluid interfaces are at the core of a broad range of technological processes, from the stabilization of responsive foams and emulsions to advanced lithographic techniques. However, establishing a fundamental relation between their internal architecture, which is controlled during synthesis, and their structural and mechanical properties upon interfacial confinement remains an elusive task. To address this open issue, which defines the monolayer's properties, we synthesize core-shell microgels, whose soft core can be chemically degraded in a controlled fashion. This strategy allows us to obtain a series of particles ranging from analogues of standard batch-synthesized microgels to completely hollow ones after total core removal. Combined experimental and numerical results show that our hollow particles have a thin and deformable shell, leading to a temperature-responsive collapse of the internal cavity and a complete flattening after adsorption at a fluid interface. Mechanical characterization shows that a critical degree of core removal is required to obtain soft disk-like particles at an oil-water interface, which present a distinct response to compression. At low packing fractions, the mechanical response of the monolayer is dominated by the outer polymer chains forming a corona surrounding the particles within the interfacial plane, regardless of the presence of a core. By contrast, at high compression, the absence of a core enables the particles to deform in the direction orthogonal to the interface and to be continuously compressed without altering the monolayer structure. These findings show how fine, single-particle architectural control during synthesis can be engineered to determine the interfacial behavior of microgels, enabling one to link particle conformation with the resulting material properties.

Interface-induced hysteretic volume phase transition of microgels: simulation and experiment

Thermo-responsive microgel particles can exhibit a drastic volume shrinkage upon increasing the solvent temperature. Recently we found that the spreading of poly(N-isopropylacrylamide) (PNiPAm) microgels at a liquid interface under the influence of surface tension hinders the temperature-induced volume phase transition. In addition, we observed a hysteresis behavior upon temperature cycling, i.e. a different evolution in microgel size and shape depending on whether the microgel was initially adsorbed to the interface in expanded or collapsed state. Here, we model the volume phase transition of such microgels at an air/water interface by monomer-resolved Brownian dynamics simulations and compare the observed behavior with experiments. We reproduce the experimentally observed hysteresis in the microgel dimensions upon temperature variation. Our simulations did not observe any hysteresis for microgels dispersed in the bulk liquid, suggesting that it results from the distinct interfacial morphology of the microgel adsorbed at the liquid interface. An initially collapsed microgel brought to the interface and subjected to subsequent swelling and collapsing (resp. cooling and heating) will end up in a larger size than it had in the original collapsed state. Further temperature cycling, however, only shows a much reduced hysteresis, in agreement with our experimental observations. We attribute the hysteretic behavior to a kinetically trapped initial collapsed configuration, which relaxes upon expanding in the swollen state. We find a similar behavior for linear PNiPAm chains adsorbed to an interface. Our combined experimental - simulation investigation provides new insights into the volume phase transition of PNiPAm materials adsorbed to liquid interfaces.

Characterization of the volume fraction of soft deformable microgels by means of small-angle neutron scattering with contrast variation

The volume occupied by colloids in a suspension - namely the volume fraction - is the thermodynamic variable that determines the phase behavior of these systems. While for hard incompressible spheres this quantity is well defined, for soft compressible colloids such as microgels - polymeric crosslinked networks swollen in a good solvent - the determination of the real volume occupied by these particles in solution is particularly challenging. This fact depends on two aspects: first the surface and, therefore, the volume of the microgels is hard to define properly given their external fuzziness; second, microgels can osmotically deswell, deform or interpenetrate their neighbors, i.e. change their shape and size depending on the solution concentration. Here, the form factors of few hydrogenated microgels embedded in a matrix of deuterated but otherwise identical microgels are measured using small-angle neutron scattering with contrast variation. From the analysis of the scattering data, the variation of the volume of the microgels as a function of concentration is obtained and used to calculate the real microgel volume fraction in solution. Soft neutral microgels are shown to facet already at low concentrations while in contrast, harder microgels maintain their shape and change their volume.

Soft Particles at Liquid Interfaces: From Molecular Particle Architecture to Collective Phase Behavior

Soft particles such as microgels can undergo significant and anisotropic deformations when adsorbed to a liquid interface. This, in turn, leads to a complex phase behavior upon compression. To date, experimental efforts have predominantly provided phenomenological links between microgel structure and resulting interfacial behavior, while simulations have not been entirely successful in reproducing experiments or predicting the minimal requirements for the desired phase behavior. Here, we develop a multiscale framework to link the molecular particle architecture to the resulting interfacial morphology and, ultimately, to the collective interfacial phase behavior. To this end, we investigate interfacial morphologies of different poly(N-isopropylacrylamide) particle systems using phase-contrast atomic force microscopy and correlate the distinct interfacial morphology with their bulk molecular architecture. We subsequently introduce a new coarse-grained simulation method that uses augmented potentials to translate this interfacial morphology into the resulting phase behavior upon compression. The main novelty of this method is the possibility to efficiently encode multibody interactions, the effects of which are key to distinguishing between heterostructural (anisotropic collapse) and isostructural (isotropic collapse) phase transitions. Our approach allows us to qualitatively resolve existing discrepancies between experiments and simulations. Notably, we demonstrate the first in silico account of the two-dimensional isostructural transition, which is frequently found in experiments but elusive in simulations. In addition, we provide the first experimental demonstration of a heterostructural transition to a chain phase in a single-component system, which has been theoretically predicted decades ago. Overall, our multiscale framework provides a phenomenological bridge between physicochemical soft-particle characteristics at the molecular scale and nanoscale and the collective self-assembly phenomenology at the macroscale, serving as a stepping stone toward an ultimately more quantitative and predictive design approach.

Effect of Chain Polydispersity on the Elasticity of Disordered Polymer Networks

Due to their unique structural and mechanical properties, randomly cross-linked polymer networks play an important role in many different fields, ranging from cellular biology to industrial processes. In order to elucidate how these properties are controlled by the physical details of the network (e.g., chain-length and end-to-end distributions), we generate disordered phantom networks with different cross-linker concentrations C and initial densities ρinit and evaluate their elastic properties. We find that the shear modulus computed at the same strand concentration for networks with the same C, which determines the number of chains and the chain-length distribution, depends strongly on the preparation protocol of the network, here controlled by ρinit. We rationalize this dependence by employing a generic stress-strain relation for polymer networks that does not rely on the specific form of the polymer end-to-end distance distribution. We find that the shear modulus of the networks is a nonmonotonic function of the density of elastically active strands, and that this behavior has a purely entropic origin. Our results show that if short chains are abundant, as it is always the case for randomly cross-linked polymer networks, the knowledge of the exact chain conformation distribution is essential for correctly predicting the elastic properties. Finally, we apply our theoretical approach to literature experimental data, qualitatively confirming our interpretations.

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.

Engulfing Behavior of Nanoparticles into Thermoresponsive Microgels: A Mesoscopic Simulation Study

The engulfing of nanoparticles into microgels provides a versatile platform to design nano- and microstructured materials with various shape anisotropies and multifunctional properties. Manipulating the spontaneous engulfment process remains elusive. Herein, we report a mesoscopic simulation study on the engulfing behavior of nanoparticles into thermoresponsive microgels. The effects of the multiple parameters, including binding strength, temperature, and nanoparticle size, are examined systematically. Our simulation results disclose three engulfing states at different temperatures, namely full-engulfing, half-engulfing, and surface contact. The engulfing depth is determined by the complementary balance of interfacial elastocapillarity. Specifically, the van der Waals interaction of hybrid microgel-nanoparticle offers the capillary force while the internally networked structure of microgel reinforces the elasticity repulsion. Our study, validated by relevant experimental results, provides a mechanistic understanding of the interfacial elastocapillarity for nanoparticle-microgels.

Charge affinity and solvent effects in numerical simulations of ionic microgels

Ionic microgel particles are intriguing systems in which the properties of thermo-responsive polymeric colloids are enriched by the presence of charged groups. In order to rationalize their properties and predict the behaviour of microgel suspensions, it is necessary to develop a coarse-graining strategy that starts from the accurate modelling of single particles. Here, we provide a numerical advancement of a recently-introduced model for charged co-polymerized microgels by improving the treatment of ionic groups in the polymer network. We investigate the thermoresponsive properties of the particles, in particular their swelling behaviour and structure, finding that, when charged groups are considered to be hydrophilic at all temperatures, highly charged microgels do not achieve a fully collapsed state, in favorable comparison to experiments. In addition, we explicitly include the solvent in the description and put forward a mapping between the solvophobic potential in the absence of the solvent and the monomer-solvent interactions in its presence, which is found to work very accurately for any charge fraction of the microgel. Our work paves the way for comparing single-particle properties and swelling behaviour of ionic microgels to experiments and to tackle the study of these charged soft particles at a liquid-liquid interface.

Mechanochemical Cellular Membrane Internalization of Nanohydrogels: A Large-Scale Mesoscopic Simulation

By combining large-scale dissipative particle dynamics and steered molecular dynamics simulations, we investigate the mechanochemical cellular internalization pathways of homogeneous and heterogeneous nanohydrogels and demonstrate that membrane internalization is determined by the crosslink density and encapsulation ability of nanohydrogels. The homogeneous nanohydrogels with a high crosslink density and low encapsulation ability behave as soft nanoparticles partially wrapped by the membrane, while those with a low crosslink density and high encapsulation ability permeate into the membrane. Regardless of the crosslink density, the homogeneous nanohydrogels undergo typical dual morphological deformations. The local lipid nanodomains are identified at the contacting region between the membrane and nanohydrogels because of different diffusion behaviors between lipid and receptor molecules during the internalization process. The yolk@shell heterogeneous nanohydrogels present a different mechanochemical cellular internalization pathway. The yolk with strong affinity is directly in contact with the membrane, resulting in partial membrane wrapping, and the contacting area is much reduced when compared to homogenous nanohydrogels, leading to a smaller lipid nanodomain and thus avoiding related cellular toxicity. Our findings provide a critical mechanism understanding of the biological pathways of nanohydrogels and may guide the molecular design of the hydrogel-based materials for controlled release drug delivery, tissue engineering, and cell culture.

Nanostructures, Thermoresponsiveness, and Assembly Mechanism of Hydrogel Microspheres during Aqueous Free-Radical Precipitation Polymerization

Although techniques to produce uniformly sized hydrogel microspheres (microgels) by aqueous free-radical precipitation polymerization are well established, the details of the polymerization process remain mysterious. In the present study, the structural evolution and thermoresponsiveness of the developing microgels during the polymerization were evaluated by temperature-controlled high-speed atomic force microscopy. This analysis clarified that the swelling properties of the precursor microgels formed in the early stages of the polymerization are quite low due to the high incorporation of cross-linkers and that non-thermoresponsive deca-nanosized spherical domains are already present in the precursor microgels. Furthermore, we succeeded in tracking the formation of nuclei and their growth process, which has never been fully understood, in aqueous solution by real-time observations. These findings will help us to design functional microgels with the desired nanostructures via precipitation polymerization.

Controlling microgel deformation via deposition method and surface functionalization of solid supports

Soft matter at solid-liquid interfaces plays an important role in multiple scientific disciplines as well as in various technological fields. For microgels, representing highly interesting soft matter systems, we demonstrate that the preparation method, i.e. the way how the microgel is applied to the specific surface, plays a key role. Focusing on the three most common sample preparation methods (spin-coating, drop-casting and adsorption from solution), we performed a comparative study of the deformation behavior of microgels at the solid-liquid interface on three different surfaces with varying hydrophilicities. For in situ visualization of the deformation of pNIPMAM microgels, we conducted highly sensitive 3D super resolution fluorescence microscopy methods. We furthermore performed complementary molecular dynamics simulations to determine the driving force responsible for the deformation depending on the surface and the deposition method. The combination of experiments and simulations revealed that the simulated equilibrium structure obtained after simulation of the completely dry microgel after deposition is retained after rehydration and subsequent fluorescent imaging.

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.

Behavior and mechanics of dense microgel suspensions

Suspensions of soft and highly deformable microgels can be concentrated far more than suspensions of hard colloids, leading to their unusual mechanical properties. Microgels can accommodate compression in suspensions in a variety of ways such as interpenetration, deformation, and shrinking. Previous experiments have offered insightful, but somewhat conflicting, accounts of the behavior of individual microgels in compressed suspensions. We develop a mesoscale computational model to probe the behavior of compressed suspensions consisting of microgels with different architectures at a variety of packing fractions and solvent conditions. We find that microgels predominantly change shape and mildly shrink above random close packing. Interpenetration is only appreciable above space filling, remaining small relative to the mean distance between cross-links. At even higher packing fractions, microgels solely shrink. Remarkably, irrespective of the single-microgel properties, and whether the suspension concentration is changed via changing the particle number density or the swelling state of the particles, which can even result in colloidal gelation, the mechanics of the suspension can be quantified in terms of the single-microgel bulk modulus, which thus emerges as the correct mechanical measure for these type of soft-colloidal suspensions. Our results rationalize the many and varied experimental results, providing insights into the relative importance of effects defining the mechanics of suspensions comprising soft particles.

Pathways and challenges towards a complete characterization of microgels

Due to their controlled size, sensitivity to external stimuli, and ease-of-use, microgel colloids are unique building blocks for soft materials made by crosslinking polymers on the micrometer scale. Despite the plethora of work published, many questions about their internal structure, interactions, and phase behavior are still open. The reasons for this lack of understanding are the challenges arising from the small size of the microgel particles, complex pairwise interactions, and their solvent permeability. Here we describe pathways toward a complete understanding of microgel colloids based on recent experimental advances in nanoscale characterization, such as super-resolution microscopy, scattering methods, and modeling.