No Arabic abstract
Colloidal gel networks are disordered elastic solids that can form even in extremely dilute particle suspensions. With interaction strengths comparable to the thermal energy, their stress-bearing network can locally restructure via breaking and reforming interparticle bonds. This allows for yielding, self-healing, and adaptive mechanics under deformation. Designing such features requires controlling stress transmission through the complex structure of the gel and this is challenging because the link between local restructuring and overall response of the network is still missing. Here, we use a space resolved analysis of dynamical processes and numerical simulations of a model gel to gain insight into this link. We show that consequences of local bond breaking propagate along the gel network over distances larger than the average mesh size. This provides the missing microscopic explanation for why nonlocal constitutive relations are necessary to rationalize the nontrivial mechanical response of colloidal gels.
Using molecular dynamics computer simulations we investigate the aging dynamics of a gel. We start from a fractal structure generated by the DLCA-DEF algorithm, onto which we then impose an interaction potential consisting of a short-range attraction as well as a long-range repulsion. After relaxing the system at T=0, we let it evolve at a fixed finite temperature. Depending on the temperature T we find different scenarios for the aging behavior. For T>0.2 the fractal structure is unstable and breaks up into small clusters which relax to equilibrium. For T<0.2 the structure is stable and the dynamics slows down with increasing waiting time. At intermediate and low T the mean squared displacement scales as t^{2/3} and we discuss several mechanisms for this anomalous time dependence. For intermediate T, the self-intermediate scattering function is given by a compressed exponential at small wave-vectors and by a stretched exponential at large wave-vectors. In contrast, for low T it is a stretched exponential for all wave-vectors. This behavior can be traced back to a subtle interplay between elastic rearrangements, fluctuations of chain-like filaments, and heterogeneity.
We use numerical simulations and an athermal quasi-static shear protocol to investigate the yielding of a model colloidal gel. Under increasing deformation, the elastic regime is followed by a significant stiffening before yielding takes place. A space-resolved analysis of deformations and stresses unravel how the complex load curve observed is the result of stress localization and that the yielding can take place by breaking a very small fraction of the network connections. The stiffening corresponds to the stretching of the network chains, unbent and aligned along the direction of maximum extension. It is characterized by a strong localization of tensile stresses, that triggers the breaking of a few network nodes at around 30% of strain. Increasing deformation favors further breaking but also shear-induced bonding, eventually leading to a large-scale reorganization of the gel structure at the yielding. At low enough shear rates, density and velocity profiles display significant spatial inhomogeneity during yielding in agreement with experimental observations.
The atomic theory of elasticity of amorphous solids, based on the nonaffine response formalism, is extended into the nonlinear stress-strain regime by coupling with the underlying irreversible many-body dynamics. The latter is implemented in compact analytical form using a qualitative method for the many-body Smoluchowski equation. The resulting nonlinear stress-strain (constitutive) relation is very simple, with few fitting parameters, yet contains all the microscopic physics. The theory is successfully tested against experimental data on metallic glasses, and it is able to reproduce the ubiquitous feature of stress-strain overshoot upon varying temperature and shear rate. A clear atomic-level interpretation is provided for the stress overshoot, in terms of the competition between the elastic instability caused by nonaffine deformation of the glassy cage and the stress buildup due to viscous dissipation.
The specific heat capacity $c_v$ of glass formers undergoes a hysteresis when subjected to a cooling-heating cycle, with a larger $c_v$ and a more pronounced hysteresis for fragile glasses than for strong ones. Here, we show that these experimental features, including the unusually large magnitude of $c_v$ of fragile glasses, are well reproduced by kinetic Monte Carlo and equilibrium study of a distinguishable particle lattice model (DPLM) incorporating a two-state picture of particle interactions. The large $c_v$ in fragile glasses is caused by a dramatic transfer of probabilistic weight from high-energy particle interactions to low-energy ones as temperature decreases.
We have made experimental observations of the force networks within a two-dimensional granular silo similar to the classical system of Janssen. Models like that of Janssen predict that pressure within a silo saturates with depth as the result of vertical forces being redirected to the walls of the silo where they can then be carried by friction. By averaging ensembles of experimentally-obtained force networks in different ways, we compare the observed behavior with various predictions for granular silos. We identify several differences between the mean behavior in our system and that predicted by Janssen-like models: We find that the redirection parameter describing how the force network transfers vertical forces to the walls varies with depth. We find that changes in the preparation of the material can cause the pressure within the silo to either saturate or to continue building with depth. Most strikingly, we observe a non-linear response to overloads applied to the top of the material in the silo. For larger overloads we observe the previously reported giant overshoot effect where overload pressure decays only after an initial increase [G. Ovarlez et al., Phys. Rev. E 67, 060302(R) (2003)]. For smaller overloads we find that additional pressure propagates to great depth. This effect depends on the particle stiffness, as given for instance by the Youngs modulus, E, of the material from which the particles are made. Important measures include E, the unscreened hydrostatic pressure, and the applied load. These experiments suggest that when the load and the particle weight are comparable, particle elasticity acts to stabilize the force network, allowing non-linear network effects to be seen in the mean behavior.