Flow heterogeneities in supercooled liquids and glasses under shear

Brittle yielding in supercooled liquids below the critical temperature of mode coupling theory

Molecular dynamics computer simulations of a polydisperse soft-sphere model under shear are presented. The starting point for these simulations are deeply supercooled samples far below the critical temperature, T c, of mode coupling theory. These samples are fully equilibrated with the aid of the swap Monte Carlo technique. For states below T c, we identify a lifetime τlt that measures the time scale on which the system can be considered as an amorphous solid. The temperature dependence of τlt can be well described by an Arrhenius law. The existence of transient amorphous solid states below T c is associated with the possibility of brittle yielding, as manifested by a sharp stress drop in the stress–strain relation and shear banding. We show that brittle yielding requires, on the one hand, low shear rates and, on the other hand, the time scale corresponding to the inverse shear rate has to be smaller or of the order of τlt. Both conditions can only be met for a large lifetime τlt, i.e., for states far below T c.

Rheological response of a glass-forming liquid having large bidispersity

Using extensive numerical simulations, we investigate the flow behaviour of a model glass-forming binary mixture whose constituent particles have a large size ratio. The rheological response to applied shear is studied in the regime where the larger species are spatially predominant. We demonstrate that the macroscopic rigidity that emerges with increasing density occurs in the regime where the larger species undergo a glass transition while the smaller species continue to be highly diffusive. We analyse the interplay between the timescale imposed by the shear and the quiescent relaxation dynamics of the two species to provide a microscopic insight into the observed rheological response. Finally, by tuning the composition of the mixture, we illustrate that the systematic insertion of the smaller particles affects the rheology by lowering of viscosity of the system.

Continuum modeling of shear startup in soft glassy materials

Yield stress fluids (YSFs) display a dual nature highlighted by the existence of a critical stress σ_{y} such that YSFs are solid for stresses σ imposed below σ_{y}, whereas they flow like liquids for σ>σ_{y}. Under an applied shear rate γ[over ̇], the solid-to-liquid transition is associated with a complex spatiotemporal scenario that depends on the microscopic details of the system, on the boundary conditions, and on the system size. Still, the general phenomenology reported in the literature boils down to a simple sequence that can be divided into a short-time response characterized by the so-called "stress overshoot," followed by stress relaxation towards a steady state. Such relaxation can be either (1) long-lasting, which usually involves the growth of a shear band that can be only transient or that may persist at steady state or (2) abrupt, in which case the solid-to-liquid transition resembles the failure of a brittle material, involving avalanches. In the present paper, we use a continuum model based on a spatially resolved fluidity approach to rationalize the complete scenario associated with the shear-induced yielding of YSFs. A key feature of our model is to provide a scaling for the coordinates of the stress overshoot, i.e., stress σ_{M} and strain γ_{M} as a function of γ[over ̇], which shows good agreement with experimental and numerical data extracted from the literature. Moreover, our approach shows that the power-law scaling σ_{M}(γ[over ̇]) is intimately linked to the growth dynamics of a fluidized boundary layer in the vicinity of the moving boundary. Yet such scaling is independent of the fate of that layer, and of the long-term behavior of the YSF, i.e., whether the steady-state flow profile is homogeneous or shear-banded. Finally, when including the presence of "long-range" correlations, we show that our model displays a ductile to brittle transition, i.e., the stress overshoot reduces into a sharp stress drop associated with avalanches, which impacts the scaling σ_{M}(γ[over ̇]). This generalized model nicely captures subtle avalanche-like features of the transient shear banding dynamics reported in experiments. Our work offers a unified picture of shear-induced yielding in YSFs, whose complex spatiotemporal dynamics are deeply connected to nonlocal effects.

On the yielding of a point-defect-rich model crystal under shear: insights from molecular dynamics simulations

In real crystals and at finite temperatures point defects are inevitable. Under shear their dynamics severely influence the mechanical properties of these crystals, giving rise to non-linear effects, such as ductility. In an effort to elucidate the complex behavior of crystals under plastic deformation it is crucial to explore and to understand the interplay between the timescale related to the equilibrium point-defect diffusion and the shear-induced timescale. Based on extensive non-equilibrium molecular dynamics simulations we present a detailed investigation on the yielding behavior of cluster crystals, an archetypical model for a defect-rich crystal: in such a system clusters of overlapping particles occupy the lattice sites of a regular (FCC) structure. In equilibrium particles diffuse via site-to-site hopping while maintaining the crystalline structure intact. We investigate these cluster crystals at a fixed density and at different temperatures where the system remains in the FCC structure: temperature allows us to vary the diffusion timescale appropriately. We then expose the crystal to shear, thereby choosing shear rates which cover timescales that are both higher and lower than the equilibrium diffusion timescales. We investigate the macroscopic and microscopic response of our cluster crystal to shear and find that the yielding scenario of such a system does not rely on the diffusion of the particles - it is rather related to the plastic deformation of the underlying crystalline structure. The local bond order parameters and the measurement of local angles between neighboring clusters confirm the cooperative movement of the clusters close to the yield point. Performing complementary, related simulations for an FCC crystal formed by harshly repulsive particles reveals similarities in the yielding behavior between both systems. Still we find that the diffusion of particles does influence characteristic features in the cluster crystal, such as a less prominent increase of order parameters close to the yield point. Our simulations provide for the first time an insight into the role of the diffusion of defects in the yielding behavior of a defect-rich crystal under shear. These observations will thus be helpful in the development of theories for the plastic deformation of defect-rich crystals.