Search

In this article, we examine the collective particle dynamics, as expressed by the time correlation function of the longitudinal particle current density, of several different fluids in the vicinity of their freezing points/lines. We consider and compare results obtained by dynamic light scattering for a suspension of hard spheres and by molecular dynamics for fluids with hard sphere and Lennard-Jones interactions. The latter are performed along both an isotherm and an isochore. In all cases, we find a qualitative change in the collective dynamics, within the resolution of the data, when their respective freezing lines are crossed. We associate this change with the onset of caging. The new results for the Lennard-Jones fluid reported here confirm that the occurrence of caging, found previously for systems of hard spheres, is a more general feature that distinguishes a metastable fluid from one in thermodynamic equilibrium. [ABSTRACT FROM AUTHOR]

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., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

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.

The growing interest in the dynamical properties of colloidal suspensions, both in equilibrium and under an external drive such as shear or pressure flow, requires the development of accurate methods to correctly include hydrodynamic effects due to the suspension in a solvent. In the present work, we generalize Multiparticle Collision Dynamics (MPCD) to be able to deal with soft, polymeric colloids. Our methods build on the knowledge of the monomer density profile that can be obtained from monomer-resolved simulations without hydrodynamics or from theoretical arguments. We hereby propose two different approaches. The first one simply extends the MPCD method by including in the simulations effective monomers with a given density profile, thus neglecting monomer-monomer interactions. The second one considers the macromolecule as a single penetrable soft colloid (PSC), which is permeated by an inhomogeneous distribution of solvent particles. By defining an appropriate set of rules to control the collision events between the solvent and the soft colloid, both linear and angular momenta are exchanged. We apply these methods to the case of linear chains and star polymers for varying monomer lengths and arm number, respectively, and compare the results for the dynamical properties with those obtained within monomer-resolved simulations. We find that the effective monomer method works well for linear chains, while the PSC method provides very good results for stars. These methods pave the way to extend MPCD treatments to complex macromolecular objects such as microgels or dendrimers and to work with soft colloids at finite concentrations.

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.