A simple and accurate method to determine fluid-crystal phase boundaries from direct coexistence simulations.

One method for computationally determining phase boundaries is to explicitly simulate a direct coexistence between the two phases of interest. Although this approach works very well for fluid-fluid coexistences, it is often considered to be less useful for fluid-crystal transitions, as additional care must be taken to prevent the simulation boundaries from imposing unwanted strains on the crystal phase. Here, we present a simple adaptation to the direct coexistence method that nonetheless allows us to obtain highly accurate predictions of fluid-crystal coexistence conditions, assuming that a fluid-crystal interface can be readily simulated. We test our approach on hard spheres, the screened Coulomb potential, and a 2D patchy-particle model. In all cases, we find excellent agreement between the direct coexistence approach and (much more cumbersome) free-energy calculation methods. Moreover, the method is sufficiently accurate to resolve the (tiny) free-energy difference between the face-centered cubic and hexagonally close-packed crystal of hard spheres in the thermodynamic limit. The simplicity of this method also ensures that it can be trivially implemented in essentially any simulation method or package. Hence, this approach provides an excellent alternative to free-energy based methods for the precise determination of phase boundaries.

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.

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.

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.

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.

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.

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.