Plastic deformations in crystal, polycrystal, and glass in binary mixtures under shear: collective yielding

Reversible-to-irreversible transition of colloidal polycrystals under cyclic athermal quasistatic deformation

Cyclic loading on granular packings and amorphous media exhibits a transition from reversible elastic behavior to irreversible plasticity. The present study compares the irreversibility transition and microscopic details of colloidal polycrystals under oscillatory tensile-compressive and shear strain. Under both modes, the systems exhibit a reversible to irreversible transition. However, the strain amplitude at which the transition is observed is larger in the shear strain than in the tensile-compressive mode. The threshold strain amplitude is confirmed by analyzing the dynamical properties, such as mobility and atomic strain (von Mises shear strain and the volumetric strain). The structural changes are quantified using a hexatic order parameter. Under both modes of deformation, dislocations and grain boundaries in polycrystals disappear, and monocrystals are formed. We also recognize the dislocation motion through grains. The key difference is that strain accumulates diagonally in oscillatory tensile-compressive deformation, whereas in shear deformation, strain accumulation is along the x or y axis.

Acoustic resonance in periodically sheared glass: damping due to plastic events

Using molecular dynamics simulation, we study acoustic resonance in a low-temperature model glass by applying a small periodic shear at a boundary wall. Shear wave resonance occurs as the frequency ω approaches ωl = πc⊥l/L (l = 1, 2…). Here, c⊥ is the transverse sound speed and L is the cell width. At resonance, large-amplitude sound waves appear after many cycles even if the applied strain γ0 is very small. They then induce plastic events, which are heterogeneous on the mesoscopic scale and intermittent on timescales longer than the oscillation period tp = 2π/ω. We visualize them together with the extended elastic strains around them. These plastic events serve to damp sounds. We obtain the nonlinear damping Q-1 = tan δ due to the plastic events near the first resonance at ω ≅ ω1, which is linear in γ0 and independent of ω. After many resonant cycles, we observe an increase in the shear modulus (measured after switching-off the oscillation). We also observe catastrophic plastic events after a very long time (∼103tp), which induce system-size elastic strains and cause a transition from resonant to off-resonant states. At resonance, stroboscopic diffusion becomes detectable.

Scaling of relaxation and excess entropy in plastically deformed amorphous solids

When stressed sufficiently, solid materials yield and deform plastically via reorganization of microscopic constituents. Indeed, it is possible to alter the microstructure of materials by judicious application of stress, an empirical process utilized in practice to enhance the mechanical properties of metals. Understanding the interdependence of plastic flow and microscopic structure in these nonequilibrium states, however, remains a major challenge. Here, we experimentally investigate this relationship, between the relaxation dynamics and microscopic structure of disordered colloidal solids during plastic deformation. We apply oscillatory shear to solid colloidal monolayers and study their particle trajectories as a function of shear rate in the plastic regime. Under these circumstances, the strain rate, the relaxation rate associated with plastic flow, and the sample microscopic structure oscillate together, but with different phases. Interestingly, the experiments reveal that the relaxation rate associated with plastic flow at time t is correlated with the strain rate and sample microscopic structure measured at earlier and later times, respectively. The relaxation rate, in this nonstationary condition, exhibits power-law, shear-thinning behavior and scales exponentially with sample excess entropy. Thus, measurement of sample static structure (excess entropy) provides insight about both strain rate and constituent rearrangement dynamics in the sample at earlier times.

Poiseuille Flow of Soft Polycrystals in 2D Rough Channels

Polycrystals are partially ordered solids where crystalline order extends over mesoscopic length scales, namely, the grain size. We study the Poisuielle flow of such materials in a rough channel. In general, similar to yield stress fluids, three distinct dynamical states, namely, flowing, stick-slip, and jammed can be observed, with a yield threshold dependent on channel width. Importantly, the interplay between the finite channel width, and the intrinsic ordering scale (the grain size) leads to a new type of spatiotemporal heterogeneity. In wide channels, although the average flow profile remains pluglike, at the underlying granular level, there is vigorous grain remodeling activity resulting from the velocity heterogeneity among the grains. As the channel width approaches typical grain size, the flowing polycrystalline state breaks up into a spatially heterogeneous mixture of flowing liquid like patches and chunks of nearly static grains. Despite these static grains, the average velocity still shows a parabolic profile, dominated by the moving liquidlike patches. However, the solid-liquid front moves at nearly constant speed in the opposite direction of the external drive.

Discrete fluidization of dense monodisperse emulsions in neutral wetting microchannels

The rheology of pressure-driven flows of two-dimensional dense monodisperse emulsions in neutral wetting microchannels is investigated by means of mesoscopic lattice Boltzmann simulations, capable of handling large collections of droplets, in the order of several hundreds. The simulations reveal that the fluidization of the emulsion proceeds through a sequence of discrete steps, characterized by yielding events whereby layers of droplets start rolling over each other, thus leading to sudden drops of the relative effective viscosity. It is shown that such discrete fluidization is robust against loss of confinement, namely it persists also in the regime of small ratios of the droplet diameter over the microchannel width. We also develop a simple phenomenological model which predicts a linear relation between the relative effective viscosity of the emulsion and the product of the confinement parameter (global size of the device over droplet radius) and the viscosity ratio between the disperse and continuous phases. The model shows excellent agreement with the numerical simulations. The present work offers new insights to enable the design of microfluidic scaffolds for tissue engineering applications and paves the way to detailed rheological studies of soft-glassy materials in complex geometries.

Theory of applying heat flow from thermostatted boundary walls: Dissipative and local-equilibrium responses and fluctuation theorems

We construct a microscopic theory of applying a heat flow from thermostatted boundary walls in the film geometry. We treat a classical one-component fluid, but our method is applicable to any fluids and solids. We express linear response of any variable B in terms of the time-correlation functions between B and the heat flows J_K from the thermostats to the particles. Furthermore, the surface variables J_K can be written in the form of space integrals of bulk quantities from the equations of motion. Owing to this surface-to-bulk relation, the steady-state response functions consist of dissipative and local-equilibrium parts, where the former gives rise to Fourier's law with Green's expression for the thermal conductivity. In the nonlinear regime, we derive the steady-state distribution in the phase space in the McLennan-Zubarev form from the first principles. Some fluctuation theorems are also presented.

Theory of applying shear strains from boundary walls: Linear response in glasses

We construct a linear response theory of applying shear deformations from boundary walls in the film geometry in Kubo's theoretical scheme. Our method is applicable to any solids and fluids. For glasses, we assume quasi-equilibrium around a fixed inherent state. Then, we obtain linear-response expressions for any variables including the stress and the particle displacements, even though the glass interior is elastically inhomogeneous. In particular, the shear modulus can be expressed in terms of the correlations between the interior stress and the forces from the walls. It can also be expressed in terms of the inter-particle correlations, as has been shown in the previous literature. Our stress relaxation function includes the effect of the boundary walls and can be used for inhomogeneous flow response. We show the presence of long-ranged, long-lived correlations among the fluctuations of the forces from the walls and the displacements of all the particles in the cell. We confirm these theoretical results numerically in a two-dimensional model glass. As an application, we describe emission and propagation of transverse sounds after boundary wall motions using these time-correlation functions. We also find resonant sound amplification when the frequency of an oscillatory shear approaches that of the first transverse sound mode.

Irreversibility transition of colloidal polycrystals under cyclic deformation

Cyclically loaded disordered particle systems, such as granular packings and amorphous media, display a non-equilibrium phase transition towards irreversibility. Here, we investigate numerically the cyclic deformation of a colloidal polycrystal with impurities and reveal a transition to irreversible behavior driven by the displacement of dislocations. At the phase transition we observe enhanced particle diffusion, system size effects and broadly distributed strain bursts. In addition to provide an analogy between the deformation of amorphous and polycrystalline materials, our results allow to reinterpret Zener pinning of grain boundaries as a way to prevent the onset of irreversible crystal ordering.

Role of disorder in finite-amplitude shear of a 2D jammed material

A material's response to small but finite deformations can reveal the roots of its response to much larger deformations. Here, we identify commonalities in the responses of 2D soft jammed solids with different amounts of disorder. We cyclically shear the materials while tracking their constituent particles, in experiments that feature a stable population of repeated structural relaxations. Using bidisperse particle sizes creates a more amorphous material, while monodisperse sizes yield a more polycrystalline one. We find that the materials' responses are very similar, both at the macroscopic, mechanical level and in the microscopic motions of individual particles. However, both locally and in bulk, crystalline arrangements of particles are stiffer (greater elastic modulus) and less likely to rearrange. Our work supports the idea of a common description for the responses of a wide array of materials.

Plasticity of a colloidal polycrystal under cyclic shear

We use confocal microscopy and time-resolved light scattering to investigate plasticity in a colloidal polycrystal, following the evolution of the network of grain boundaries as the sample is submitted to thousands of shear deformation cycles. The grain boundary motion is found to be ballistic, with a velocity distribution function exhibiting nontrivial power law tails. The shear-induced dynamics initially slow down, similarly to the aging of the spontaneous dynamics in glassy materials, but eventually reach a steady state. Surprisingly, the crossover time between the initial aging regime and the steady state decreases with increasing probed length scale, hinting at a hierarchical organization of the grain boundary dynamics.

Slow relaxations and stringlike jump motions in fragile glass-forming liquids: breakdown of the Stokes-Einstein relation

We perform molecular dynamics simulation on a glass-forming liquid binary mixture with the soft-core potential in three dimensions. We investigate crossover of the configuration changes caused by stringlike jump motions. With lowering the temperature T, the motions of the particles composing strings become larger in sizes and displacements, while those of the particles surrounding strings become smaller. Then the contribution of the latter to time-correlation functions tends to be long-lived as T is lowered. As a result, the relaxation time τ(α) and the viscosity η grow more steeply than the inverse diffusion constant D(-1) at low T, leading to breakdown of the Stokes-Einstein relation. At low T, the diffusion occurs as activation processes and may well be described by short-time analysis of rare jump motions with broken bonds and large displacements. Some characteristic features of the Van Hove self-correlation function arise from escape jumps over high potential barriers. We also visualize the particle motions at string formation taking place in a very short time.

Relationship between bond-breakage correlations and four-point correlations in heterogeneous glassy dynamics: configuration changes and vibration modes

We investigate the dynamic heterogeneities of glassy particle systems in the theoretical schemes of bond breakage and four-point correlation functions. In the bond-breakage scheme, we introduce the structure factor S(b)(q,t) and the susceptibility χ(b)(t) to detect the spatial correlations of configuration changes. Here χ(b)(t) attains a maximum at t=t(b)(max) as a function of time t, where the fraction of the particles with broken bonds φ(b)(t) is about 1/2. In the four-point scheme, treating the structure factor S(4)(q,t) and the susceptibility χ(4)(t), we detect superpositions of the heterogeneity of bond breakage and that of thermal low-frequency vibration modes. While the former grows slowly, the latter emerges quickly to exhibit complex space-time behavior. In two dimensions, the vibration modes extending over the system yield significant contributions to the four-point correlations, which depend on the system size logarithmically. A maximum of χ(4)(t) is attained at t=t(4)(max), where these two contributions become of the same order. As a result, t(4)(max) is considerably shorter than t(b)(max).

The physics of the colloidal glass transition

As one increases the concentration of a colloidal suspension, the system exhibits a dramatic increase in viscosity. Beyond a certain concentration, the system is said to be a colloidal glass; structurally, the system resembles a liquid, yet motions within the suspension are slow enough that it can be considered essentially frozen. For several decades, colloids have served as a valuable model system for understanding the glass transition in molecular systems. The spatial and temporal scales involved allow these systems to be studied by a wide variety of experimental techniques. The focus of this review is the current state of understanding of the colloidal glass transition, with an emphasis on experimental observations. A brief introduction is given to important experimental techniques used to study the glass transition in colloids. We describe features of colloidal systems near and in glassy states, including increases in viscosity and relaxation times, dynamical heterogeneity and ageing, among others. We also compare and contrast the glass transition in colloids to that in molecular liquids. Other glassy systems are briefly discussed, as well as recently developed synthesis techniques that will keep these systems rich with interesting physics for years to come.

Inhomogeneous shear flows in soft jammed materials with tunable attractive forces

We perform molecular dynamics simulations to characterize the occurrence of inhomogeneous shear flows in soft jammed materials. We use rough walls to impose a simple shear flow and study the athermal motion of jammed assemblies of soft particles in two spatial dimensions, both for purely repulsive interactions and in the presence of an additional short-range attraction of varying strength. In steady state, pronounced flow inhomogeneities emerge for all systems when the shear rate becomes small. Deviations from linear flow are stronger in magnitude and become very long lived when the strength of the attraction increases, but differ from permanent shear bands. Flow inhomogeneities occur in a stress window bounded by the dynamic and static yield stress values. Attractive forces enhance the flow heterogeneities because they accelerate stress relaxation, thus effectively moving the system closer to the yield stress regime where inhomogeneities are most pronounced. The present scenario for understanding the effect of particle adhesion on shear localization, which is based on detailed molecular dynamics simulations with realistic particle interactions, differs qualitatively from previous qualitative explanations and ad hoc theoretical modeling.

Plasticity and dynamical heterogeneity in driven glassy materials

Many amorphous glassy materials exhibit complex spatio-temporal mechanical response and rheology, characterized by an intermittent stress strain response and a fluctuating velocity profile. Under quasistatic and athermal deformation protocols this heterogeneous plastic flow was shown to be composed of plastic events of various sizes, ranging from local quadrupolar plastic rearrangements to system spanning shear bands. In this paper, through numerical study of a 2D Lennard-Jones amorphous solid, we generalize the study of the heterogeneous dynamics of glassy materials to the finite shear rate (gamma not equal to 0) and temperature case (T not equal to 0). In practice, we choose an effectively athermal limit (T approximately 0) and focus on the influence of shear rate on the rheology of the glass. In line with previous works we find that the model Lennard-Jones glass follows the rheological behavior of a yield stress fluid with a Herschel-Bulkley response of the form, sigma = sigmaY + c1gamma(beta). The global mechanical response obtained through the use of Molecular Dynamics is shown to converge in the limit gamma --> 0 to the quasistatic limit obtained with an energy minimization protocol. The detailed analysis of the plastic deformation at different shear rates shows that the glass follows different flow regimes. At sufficiently low shear rates the mechanical response reaches a shear-rate-independent regime that exhibits all the characteristics of the quasistatic response (finite-size effects, cascades of plastic rearrangements, yield stress, ...). At intermediate shear rates the rheological properties are determined by the externally applied shear rate and the response deviates from the quasistatic limit. Finally at higher shear the system reaches a shear-rate-independent homogeneous regime. The existence of these three regimes is also confirmed by the detailed analysis of the atomic motion. The computation of the four-point correlation function shows that the transition from the shear-rate-dominated to the quasistatic regime is accompanied by the growth of a dynamical cooperativity length scale xi that is shown to diverge with shear rate as xi is proportional to gamma(-nu), with nu approximately 0.2 -0.3. This scaling is compared with the prediction of a simple model that assumes the diffusive propagation of plastic events.