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

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.

Critical behavior of quorum-sensing active particles.

It is still a debated issue whether all critical active particles belong to the same universality class. Here we numerically study the critical behavior of quorum sensing active particles that represents the archetypal model for interpreting motility-induced phase separation. Mean-field theory predicts that this model should undergo a full phase separation if particles slow-down enough when sensing the presence of their neighbors and that the coexistence line terminates in a critical point. By performing large-scale numerical simulations, we confirm this scenario, locate the critical point and use finite-size scaling analysis to show that the static and dynamic critical exponents of this active system substantially agree with those of the Ising universality class.

Universality class of the motility-induced critical point in large scale off-lattice simulations of active particles.

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.

Effective potentials induced by mixtures of patchy and hard co-solutes.

The addition of co-solutes to colloidal suspensions is often employed to induce tunable depletion interactions. In this work, we investigate effective colloidal interactions arising from binary co-solute mixtures of hard spheres and patchy particles. By changing the relative concentration of the two species, we show that the resulting effective potential Veff continuously changes from the one obtained for a single-component hard sphere co-solute to that mediated by the single-component patchy particle co-solute. Interestingly, we find that, independent of the relative concentration of the two components, the resulting Veff is additive, i.e., it is well-described by the linear combination of the effective interactions mediated by respective pure co-solutes. However, a breakdown of the additivity occurs when the co-solute mixture is close to the onset of a demixing transition. These results represent a step forward in understanding and predicting colloidal behavior in complex and crowded environments and for exploiting this knowledge to design targeted colloidal superstructures.

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.

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.

Numerical insights on ionic microgels: structure and swelling behaviour.

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.

Rheological investigation of gels formed by competing interactions: A numerical study.

A transition from solid-like to liquid-like behavior occurs when colloidal gels are subjected to a prolonged exposure to a steady shear. This phenomenon, which is characterized by a yielding point, is found to be strongly dependent on the packing fraction. However, it is not yet known how the effective inter-particle potential affects this transition. To this aim, we present a numerical investigation of the rheology of equilibrium gels in which a short-range depletion is complemented by a long-range electrostatic interaction. We observe a single yielding event in the stress-strain curve, occurring at a fixed strain. The stress overshoot is found to follow a power-law dependence on the Péclet number, with an exponent larger than that found in depletion gels, suggesting that its value may depend systematically on the underlying colloid-colloid interactions. We also establish a mapping between equilibrium states and steady states under shear, which allows us to identify the structural modifications induced by the presence of the shear. Remarkably, we find that steady states corresponding to the same Péclet number, obtained by different combinations of shear rate and solvent viscosity, show identical structural and rheological properties. Our results highlight the importance of understanding the coupling between colloidal interactions, solvent effects, and flow to be able to describe the microscopic organization of colloidal particles under shear.

Effective potentials induced by self-assembly of patchy particles.

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 Veff 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, Veff 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.

Discontinous change from thermally- to geometrically-dominated effective interactions in colloidal solutions.

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.

How soft repulsion enhances the depletion mechanism.

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.

Multiple glass singularities and isodynamics in a core-softened model for glass-forming systems.

We investigate the slow dynamics of a simple glass former whose interaction potential is the sum of a hard core and a square shoulder repulsion. According to mode coupling theory, the competition between the two repulsive length scales gives rise to a complex dynamic scenario: besides the fluid-glass line, the theory predicts a glass-glass line in the temperature-packing fraction plane with two end points. Interestingly, for critical values of the square-shoulder parameters, such end points can be accessed from the liquid phase. We verify, via extensive numerical simulations, the existence of both points through the observation of an unconventional subdiffusive/logarithmic dynamical behavior. Unexpectedly, we also discover that the simultaneous presence of two end points generates special loci in the state diagram along which the dynamics is identical at all length and time scales.

Observable-dependence of the effective temperature in off-equilibrium diatomic molecular liquids.

We discuss the observable-dependence of the effective temperature Teff, defined via the fluctuation-dissipation relation, of an out-of-equilibrium system composed by homonuclear dumbbell molecules. Teff is calculated by evaluating the fluctuation and the response for two observables associated, respectively, to translational and to rotational degrees of freedom, following a sudden temperature quench. We repeat our calculations for different dumbbell elongations ζ. At high elongations (ζ > 0.4), we find the same Teff for the two observables. At low elongations (ζ ⩽ 0.4), only for very deep quenches Teff coincides. The observable-dependence of Teff for low elongations and shallow quenches stresses the importance of a strong coupling between orientational and translational variables for a consistent definition of the effective temperature in glassy systems.

Generalized fluctuation-dissipation relation and effective temperature upon heating a deeply supercooled liquid.

We show that a generalized fluctuation-dissipation relation applies upon instantaneously increasing the temperature of a deeply supercooled liquid. This has the same two-step shape of the relation found upon cooling the liquid, but with opposite violation, indicating an effective temperature that is lower than bath temperature. We show that the effective temperature exhibits some sensible time dependence and that it retains its connection with the partitioned phase space visited in aging. We underline the potential relevance of our numerical results for experimental studies of the fluctuation-dissipation relation in glassy systems.

Unveiling the complex glassy dynamics of square shoulder systems: simulations and theory.

We performed extensive molecular dynamics (MD) simulations, supplemented by Mode Coupling Theory (MCT) calculations, for the square shoulder model, a purely repulsive potential where the hardcore is complemented by a finite shoulder. For the one-component version of this model, MCT predicted [Sperl et al., Phys. Rev. Lett. 104, 145701 (2010)] the presence of diffusion anomalies both upon cooling and upon compression and the occurrence of glass-glass transitions. In the simulations, we focus on a non-crystallising binary mixture, which, at the investigated shoulder width, shows a non-monotonic behaviour of the diffusion upon cooling but not upon isothermal compression. In addition, we find the presence of a disconnected glass-glass line in the phase diagram, ending in two higher order singularities. These points generate a logarithmic dependence of the density correlators as well as a subdiffusive behaviour of the mean squared displacement, although with the interference of the nearby liquid-glass transition. We also perform novel MCT calculations using as input the partial structure factors obtained within MD, confirming the simulation results. The presence of two hard sphere glasses, differing only in their hardcore length, is revealed, showing that the simple competition between the two is sufficient for creating a rather complex dynamical behaviour.

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.

Tuning effective interactions close to the critical point in colloidal suspensions.

We report a numerical investigation of two colloids immersed in a critical solvent, with the aim of quantifying the effective colloid-colloid interaction potential. By turning on an attraction between the colloid and the solvent particles we follow the evolution from the case in which the solvent density close to the colloids changes from values smaller than the bulk to values larger than the bulk. We thus effectively implement the so-called (+, +) and (-, -) boundary conditions defined in field theoretical approaches focused on the description of critical Casimir forces. We find that the effective potential at large distances decays exponentially, with a characteristic decay length compatible with the bulk critical correlation length, in full agreement with theoretical predictions. We also investigate the case of (+, -) boundary condition, where the effective potential becomes repulsive. Our study provides a guidance for a design of the interaction potential which can be exploited to control the stability of colloidal systems.

Properties of patchy colloidal particles close to a surface: a Monte Carlo and density functional study.

We investigate the behavior of a patchy particle model close to a hard-wall via Monte Carlo simulation and density functional theory (DFT). Two DFT approaches, based on the homogeneous and inhomogeneous versions of Wertheim's first order perturbation theory for the association free energy are used. We evaluate, by simulation and theory, the equilibrium bulk phase diagram of the fluid and analyze the surface properties for two isochores, one of which is close to the liquid side of the gas-liquid coexistence curve. We find that the density profile near the wall crosses over from a typical high-temperature adsorption profile to a low-temperature desorption one, for the isochore close to coexistence. We relate this behavior to the properties of the bulk network liquid and find that the theoretical descriptions are reasonably accurate in this regime. At very low temperatures, however, an almost fully bonded network is formed, and the simulations reveal a second adsorption regime which is not captured by DFT. We trace this failure to the neglect of orientational correlations of the particles, which are found to exhibit surface induced orientational order in this regime.

Pressure-energy correlations in liquids. V. Isomorphs in generalized Lennard-Jones systems.

This series of papers is devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids with more than 90% correlation between their virial W and potential energy U fluctuations in the NVT ensemble. Paper IV [N. Gnan et al., J. Chem. Phys. 131, 234504 (2009)] showed that strongly correlating liquids have "isomorphs," which are curves in the phase diagram along which structure, dynamics, and some thermodynamic properties are invariant in reduced units. In the present paper, using the fact that reduced-unit radial distribution functions are isomorph invariant, we derive an expression for the shapes of isomorphs in the WU phase diagram of generalized Lennard-Jones systems of one or more types of particles. The isomorph shape depends only on the Lennard-Jones exponents; thus all isomorphs of standard Lennard-Jones systems (with exponents 12 and 6) can be scaled onto a single curve. Two applications are given. One tests the prediction that the solid-liquid coexistence curve follows an isomorph by comparing to recent simulations by Ahmed and Sadus [J. Chem. Phys. 131, 174504 (2009)]. Excellent agreement is found on the liquid side of the coexistence curve, whereas the agreement is less convincing on the solid side. A second application is the derivation of an approximate equation of state for generalized Lennard-Jones systems by combining the isomorph theory with the Rosenfeld-Tarazona expression for the temperature dependence of the potential energy on isochores. It is shown that the new equation of state agrees well with simulations.

A single-agent extension of the SIR model describes the impact of mobility restrictions on the COVID-19 epidemic.

Mobility restrictions are successfully used to contain the diffusion of epidemics. In this work we explore their effect on the epidemic growth by investigating an extension of the Susceptible-Infected-Removed (SIR) model in which individual mobility is taken into account. In the model individual agents move on a chessboard with a Lévy walk and, within each square, epidemic spreading follows the standard SIR model. These simple rules allow to reproduce the sub-exponential growth of the epidemic evolution observed during the Covid-19 epidemic waves in several countries and which cannot be captured by the standard SIR model. We show that we can tune the slowing-down of the epidemic spreading by changing the dynamics of the agents from Lévy to Brownian and we investigate how the interplay among different containment strategies mitigate the epidemic spreading. Finally we demonstrate that we can reproduce the epidemic evolution of the first and second COVID-19 waves in Italy using only 3 parameters, i.e , the infection rate, the removing rate, and the mobility in the country. We provide an estimate of the peak reduction due to imposed mobility restrictions, i. e., the so-called flattening the curve effect. Although based on few ingredients, the model captures the kinetic of the epidemic waves, returning mobility values that are consistent with a lock-down intervention during the first wave and milder limitations, associated to a weaker peak reduction, during the second wave.