Search

Thermoresponsive hydrogels exhibit reversible deswelling at the Volume Phase Transition (VPT), associated to a minimum of the Poisson's ratio $\nu$. Recent numerical investigations uncovered the occurrence of a Hyper-Auxetic Transition (HAT) ($\nu=-1$) for low-crosslinked hydrogels at low temperature, accompanied by a critical-like behavior. Here, we perform extensive numerical simulations to unveil the relation between these two transitions. We find that the HAT occurs at different temperatures $T$ up to a maximum value of the crosslinker concentration $c$, thus being clearly distinct from the VPT, taking place at fixed $T$ and for all $c$. Our results provide novel fundamental insights on the interplay between network collapse and mechanical properties of these fascinating materials.

Low-connectivity polymer networks were recently found to behave auxetically when subjected to small tensions, that is, their Poisson's ratio $\nu$ becomes negative. In addition, for specific state points, numerical simulations revealed that diamond-like networks reach the limit of mechanical stability, exhibiting values of $\nu = -1$, a condition that we define hyper-auxeticity. This behavior is interesting per se for its consequences in material science, but also appealing for fundamental physics because the mechanical instability is accompanied by evidences 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 a 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, that is found to be a relevant control parameter determining the elastic properties of the disordered system. The validity of the results in disordered conditions paves the way to an experimental investigation of this phenomenon in real systems, such as hydrogels.

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 crosslinked 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 crosslinks 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 crosslink localization length to the shear modulus of the system., Comment: 37 pages, 12 figures; fixed equation at pag. 24 and typos

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 $\nu$. The amount of deformation in response to the applied force can be at most equal to the imposed one, so that $\nu=-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.

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 [Gnan \& 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 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., Comment: 13 pages, 15 figures

Thermoresponsive microgels are one of the most investigated class 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. Here we fill this gap by combining experiments and 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 several microgels with different cross-linker concentrations, charge contents and charge distributions, we build a unifying master-curve able to predict the two-step deswelling. Our results have direct relevance on fundamental soft condensed matter science and on microgel applications ranging from materials to biomedical technologies., Comment: main article: 12 pages, 7 figures; supplementary materials: 11 pages, 12 figures, 2 tables

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., Comment: 13 pages, 9 figures

Due to their unique structural and mechanical properties, randomly-crosslinked 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 (\textit{e.g.} chain-length and end-to-end distributions), we generate disordered phantom networks with different crosslinker concentrations $C$ and initial density $\rho_{\rm 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 $\rho_{\rm init}$. We rationalise 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 non-monotonic function of the density of elastically-active strands, and that this behaviour has a purely entropic origin. Our results show that if short chains are abundant, as it is always the case for randomly-crosslinked polymer networks, the knowledge of the exact chain conformation distribution is essential for predicting correctly the elastic properties. Finally, we apply our theoretical approach to published experimental data, qualitatively confirming our interpretations., Comment: SM appended after the list of references

Soft colloids are increasingly used as model systems to address fundamental issues such as crystallisation and the glass and jamming transitions. Among the available classes of soft colloids, microgels are emerging as the gold standard. Since their great internal complexity makes their theoretical characterization very hard, microgels are commonly modelled, at least in the small-deformation regime, within the simple framework of linear elasticity theory. Here we show that there exist conditions where its range of validity can be greatly extended, providing strong numerical evidence that microgels adsorbed at an interface follow the two-dimensional Hertzian theory, and hence behave like 2D elastic particles, up to very large deformations, in stark contrast to what found in bulk conditions. We are also able to estimate the Young's modulus of the individual particles and, by comparing it with its counterpart in bulk conditions, we demonstrate a significant stiffening of the polymer network at the interface. Finally, by analyzing dynamical properties, we predict multiple reentrant phenomena: by a continuous increase of particle density, microgels first arrest and then re-fluidify due to the high penetrability of their extended coronas. We observe this anomalous behavior in a range of experimentally accessible conditions for small and loosely crosslinked microgels. The present work thus establishes microgels at interfaces as a new model system for fundamental investigations, paving the way for the experimental synthesis and research on unique high-density liquidlike states. In addition, these results can guide the development of novel assembly and patterning strategies on surfaces and the design of novel materials with desired interfacial behavior.

We perform stringent tests of thermodynamic theories of the glass transition over the experimentally relevant temperature regime for several simulated glass-formers. The swap Monte Carlo algorithm is used to estimate the configurational entropy and static point-to-set lengthscale, and careful extrapolations are used for the relaxation times. We first quantify the relation between configurational entropy and the point-to-set lengthscale in two and three dimensions. We then show that the Adam-Gibbs relation is generally violated in simulated models for the experimentally relevant time window. Collecting experimental data for several supercooled molecular liquids, we show that the same trends are observed experimentally. Deviations from the Adam-Gibbs relation remain compatible with random first order transition theory, and may account for the reported discrepancies between Kauzmann and Vogel-Fulcher-Tammann temperatures. Alternatively, they may also indicate that even near $T_g$ thermodynamics is not the only driving force for slow dynamics., Comment: 13 pages, 8 figures

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\"uckel 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.

Soft particles display highly versatile properties with respect to hard colloids, even more so at fluid-fluid interfaces. In particular, microgels, consisting of a cross-linked polymer network, are able to deform and flatten upon adsorption at the interface due to the balance between surface tension and internal elasticity. Despite the existence of experimental results, a detailed theoretical understanding of this phenomenon is still lacking due to the absence of appropriate microscopic models. In this work, we propose an advanced modelling of microgels at a flat water/oil interface. The model builds on a realistic description of the internal polymeric architecture and single-particle properties of the microgel and is able to reproduce its experimentally observed shape at the interface. Complementing molecular dynamics simulations with in-situ cryo-electron microscopy experiments and atomic force microscopy imaging after Langmuir-Blodgett deposition, we compare the morphology of the microgels for different values of the cross-linking ratios. Our model allows for a systematic microscopic investigation of soft particles at fluid interfaces, which is essential to develop predictive power for the use of microgels in a broad range of applications, including the stabilization of smart emulsions and the versatile patterning of surfaces.

Thermoresponsive microgels are soft colloids that find widespread use as model systems for soft matter physics. Their complex internal architecture, made of a disordered and heterogeneous polymer network, has been so far a major challenge for computer simulations. In this work we put forward a coarse-grained model of microgels whose structural properties are in quantitative agreement with results obtained with small-angle X-ray scattering experiments across a wide range of temperatures, encompassing the volume phase transition. These results bridge the gap between experiments and simulations of individual microgel particles, paving the way to theoretically address open questions about their bulk properties with unprecedented nano and microscale resolution.

We study the equilibrium statistical properties of the potential energy landscape of several glass models in a temperature regime so far inaccessible to computer simulations. We show that unstable modes of the stationary points undergo a localization transition in real space close to the mode-coupling crossover temperature determined from the dynamics. The concentration of localized unstable modes found at low temperature is a non-universal, finite dimensional feature not captured by mean-field glass theory. Our analysis reconciles, and considerably expands, previous conflicting numerical results and provides a characteristic temperature for glassy dynamics that unambiguously locates the mode-coupling crossover., Comment: 20 pages, 11 figures, submission to SciPost. Data and workflow relevant to this work are available at https://doi.org/10.5281/zenodo.1478600

An important open problem in materials science is whether a direct connection exists between single-particle elastic properties and the macroscopic bulk behavior. Here we address this question by focusing on the archetype of soft colloids: thermoresponsive microgels. These colloidal-sized polymer networks are often assumed to interact through a simple Hertzian potential, a classic model in linear elasticity theory. By developing an appropriate methodology, that can be generalized to any kind of soft particle, we are able to calculate all the elastic moduli of non-ionic microgels across their volume phase transition (VPT). Remarkably, we reproduce many features seen in experiments, including the appearance of a minimum of the Poisson's ratio close to the VPT. By calculating the particle-particle effective interactions and the resulting collective behavior, we find that the Hertzian model works well up to moderate values of the packing fraction.

The nature of the glass transition is theoretically understood in the mean-field limit of infinite spatial dimensions, but the problem remains totally open in physical dimensions. Nontrivial finite-dimensional fluctuations are hard to control analytically, and experiments fail to provide conclusive evidence regarding the nature of the glass transition. Here, we use Monte Carlo simulations that fully bypass the glassy slowdown, and access equilibrium states in two-dimensional glass-forming liquids at low enough temperatures to directly probe the transition. We find that the liquid state terminates at a thermodynamic glass transition at zero temperature, which is associated with an entropy crisis and a diverging static correlation length., Comment: 23 pages, 18 figures

We numerically study the evolution of the vibrational density of states $D(\omega)$ of zero-temperature glasses when their kinetic stability is varied over an extremely broad range, ranging from poorly annealed glasses obtained by instantaneous quenches from above the onset temperature, to ultrastable glasses obtained by quenching systems thermalised below the experimental glass temperature. The low-frequency part of the density of states splits between extended and quasi-localized modes. Extended modes exhibit a boson peak crossing over to Debye behaviour ($D(\omega) \sim \omega^2$) at low-frequency, with a strong correlation between the two regimes. Quasi-localized modes instead obey $D(\omega) \sim \omega^4$, irrespective of the glass stability. However, the prefactor of this quartic law becomes smaller in more stable glasses, and the corresponding modes become more localized and sparser. Our work is the first numerical observation of quasi-localized modes in a regime relevant to experiments, and it establishes a direct connection between glass stability and soft vibrational motion in amorphous solids., Comment: 8 pages, 6 figures

Successful computer studies of glass-forming materials need to overcome both the natural tendency to structural ordering and the dramatic increase of relaxation times at low temperatures. We present a comprehensive analysis of eleven glass-forming models to demonstrate that both challenges can be efficiently tackled using carefully designed models of size polydisperse supercooled liquids together with an efficient Monte Carlo algorithm where translational particle displacements are complemented by swaps of particle pairs. We study a broad range of size polydispersities, using both discrete and continuous mixtures, and we systematically investigate the role of particle softness, attractivity and non-additivity of the interactions. Each system is characterized by its robustness against structural ordering and by the efficiency of the swap Monte Carlo algorithm. We show that the combined optimisation of the potential's softness, polydispersity and non-additivity leads to novel computer models with excellent glass-forming ability. For such models, we achieve over ten orders of magnitude gain in the equilibration timescale using the swap Monte Carlo algorithm, thus paving the way to computational studies of static and thermodynamic properties under experimental conditions. In addition, we provide microscopic insights into the performance of the swap algorithm which should help optimizing models and algorithms even further., Comment: 22 pages, 15 figs

Liquids relax extremely slowly upon approaching the glass state. One explanation is that an entropy crisis, due to the rarefaction of available states, makes it increasingly arduous to reach equilibrium in that regime. Validating this scenario is challenging, because experiments offer limited resolution, while numerical studies lag more than eight orders of magnitude behind experimentally-relevant timescales. In this work we not only close the colossal gap between experiments and simulations but manage to create in-silico configurations that have no experimental analog yet. Deploying a range of computational tools, we obtain four estimates of their configurational entropy. These measurements consistently confirm that the steep entropy decrease observed in experiments is also found in simulations, even beyond the experimental glass transition. Our numerical results thus extend the new observational window into the physics of glasses and reinforce the relevance of an entropy crisis for understanding their formation., Comment: 4+23 pages, 3+12 figures; v2: final version, with various changes made. Data relevant to this work can be accessed at http://dx.doi.org/10.7924/G8ZG6Q9T

We implement and optimize a particle-swap Monte-Carlo algorithm that allows us to thermalize a polydisperse system of hard spheres up to unprecedentedly-large volume fractions, where \revise{previous} algorithms and experiments fail to equilibrate. We show that no glass singularity intervenes before the jamming density, which we independently determine through two distinct non-equilibrium protocols. We demonstrate that equilibrium fluid and non-equilibrium jammed states can have the same density, showing that the jamming transition cannot be the end-point of the fluid branch., Comment: 5 pages, 3 figs; To be published in Phys. Rev. Lett

We perform Monte-Carlo simulations to analyse the structure and microscopic dynamics of a viscous Lennard-Jones liquid coupled to a quenched reference configuration of the same liquid. The coupling between the two replicas is introduced via a field epsilon conjugate to the overlap Q between the two particle configurations. This allows us to study the evolution of various static and dynamic correlation functions across the (epsilon, T) equilibrium phase diagram. As the temperature is decreased, we identify increasingly marked precursors of a first-order phase transition between a low-Q and a high-Q phase induced by the field epsilon. We show in particular that both static and dynamic susceptibilities have a maximum at a temperature-dependent value of the coupling field, which defines a `Widom line'. We also show that, in the high-overlap regime, diffusion and structural relaxation are strongly decoupled because single particle motion mostly occurs via discrete hopping on the sites defined by the reference configuration. These results, obtained using conventional numerical tools, provide encouraging signs that an equilibrium phase transition exists in coupled viscous liquids, but also demonstrate that important numerical challenges must be overcome to obtain more conclusive numerical evidence., Comment: 14 pages, 8 figures. Accepted for publication in Molecular Physics (Special Issue in honour of J.-P. Hansen)

We discuss the observable-dependence of the effective temperature $T_{eff}$, defined via the fluctuation-dissipation relation, of an out-of-equilibrium system composed by homonuclear dumbbell molecules. $T_{eff}$ 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 calculation for different dumbbell elongations $\zeta$. At high elongations ($\zeta > 0.4 $), we find the same $T_{eff}$ for the two observables. At low elongations ($\zeta \leq 0.4$), only for very deep quenches $T_{eff}$ coincides. The observable-dependence of $T_{eff}$ 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.