Shear banding and flow-concentration coupling in colloidal glasses

Permanent shear localization in dense disordered materials due to microscopic inertia

In this work using computer simulations of 3D model of dense disordered solids we show, for the first time, the appearance of shear localization in the stationary flow under homogeneous driving conditions. To rationalize our simulation results we develop a continuum model, that couples the dynamics of the local flow to the evolution of a kinetic temperature field related to the local inertial dynamics. Our model predicts that the coupling of the flow field to this additional destabilizing field appears only as a necessary condition for shear localization, a minimum system size is necessary to accommodate the flow instability. Moreover we show that this size criterion resulting from our continuum description is in quantitative agreement with our particle-based simulation results.

Dense bidisperse suspensions under non-homogeneous shear

We study the rheological behaviour of bidisperse suspensions in three dimensions under a non-uniform shear flow, made by the superimposition of a linear shear and a sinusoidal disturbance. Our results show that (i) only a streamwise disturbance in the shear-plane alters the suspension dynamics by substantially reducing the relative viscosity, (ii) with the amplitude of the disturbance determining a threshold value for the effect to kick-in and its wavenumber controlling the amount of reduction and which of the two phases is affected. We show that, (iii) the rheological changes are caused by the effective separation of the two phases, with the large or small particles layering in separate regions. We provide a physical explanation of the phase separation process and of the conditions necessary to trigger it. We test the results in the whole flow curve, and we show that the mechanism remains substantially unaltered, with the only difference being the nature of the interactions between particles modified by the phase separation.

Mechanical theory of nonequilibrium coexistence and motility-induced phase separation

Nonequilibrium phase transitions are routinely observed in both natural and synthetic systems. The ubiquity of these transitions highlights the conspicuous absence of a general theory of phase coexistence that is broadly applicable to both nonequilibrium and equilibrium systems. Here, we present a general mechanical theory for phase separation rooted in ideas explored nearly a half-century ago in the study of inhomogeneous fluids. The core idea is that the mechanical forces within the interface separating two coexisting phases uniquely determine coexistence criteria, regardless of whether a system is in equilibrium or not. We demonstrate the power and utility of this theory by applying it to active Brownian particles, predicting a quantitative phase diagram for motility-induced phase separation in both two and three dimensions. This formulation additionally allows for the prediction of novel interfacial phenomena, such as an increasing interface width while moving deeper into the two-phase region, a uniquely nonequilibrium effect confirmed by computer simulations. The self-consistent determination of bulk phase behavior and interfacial phenomena offered by this mechanical perspective provide a concrete path forward toward a general theory for nonequilibrium phase transitions.

Directed colloidal assembly and banding via DC electrokinetics

Manipulating the transport and assembly of colloidal particles to form segregated bands or ordered supracolloidal structures plays an important role in many aspects of science and technology, from understanding the origin of life to synthesizing new materials for next-generation manufacturing, electronics, and therapeutics. One commonly used method to direct colloidal transport and assembly is the application of electric fields, either AC or DC, due to its feasibility. However, as colloidal segregation and assembly both require active redistribution of colloidal particles across multiple length scales, it is not apparent at first sight how a DC electric field, either externally applied or internally induced, can lead to colloidal structuring. In this Perspective, we briefly review and highlight recent advances and standing challenges in colloidal transport and assembly enabled by DC electrokinetics.

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.

Volume-shear coupling in a mesoscopic model of amorphous materials

We present a two-dimensional mesoscopic model of a yield stress material that includes the possibility of local volume fluctuations coupled to shear in such a way that the shear strength of the material decreases as the local density decreases. The model reproduces a number of effects well known in the phenomenology of this kind of material. In particular, we find that the volume of the sample increases as the deformation rate increases; shear bands are no longer oriented at 45^{∘} with respect to the principal axis of the applied stress (as in the absence of volume-shear coupling); and homogeneous deformation becomes unstable at low enough deformation rates if volume-shear coupling is strong enough. We also discuss the effect of this coupling on some out-of-equilibrium configurations, which can be relevant to the study of the shear bands observed in metallic glasses.

Elastic criterion for shear-banding instability in amorphous solids

In amorphous solids, plastic flow is prone to localization into shear bands via an avalanche of shear-transformation (ST) rearrangements of constituent atoms or particles. However, such banding instability still remains a lack of direct experimental evidence. Using a real 3D colloidal glass under shear as proof of principle, we study STs' avalanches into shear banding that is controlled by strain rates. We demonstrate that, accompanying the emergent shear banding, the elastic response fields of the system, typical of a quadrupole for shear and a centrosymmetry for dilatation, lose the Eshelby-type spatial symmetry; instead, a strong correlation appears preferentially along the banding direction. By quantifying the fields' spatial decay, we identify an elastic criterion for the shear-banding instability, that is, the strongly correlated length of dilatation is smaller than the full length of shear correlation. Specifically, ST-induced free volume has to be confined within the elastic shear domain of ST so that those STs can self-organize to trigger shear banding. This physical picture is directly visualized by tracing the real-space evolution of local dilatation and ST particles. The present work unites the two classical mechanisms: free volume and STs, for the fundamental understanding of shear banding in amorphous solids.

Ferrihydrite nanoparticles entrapped in shear-induced multilamellar vesicles

HYPOTHESIS

Ferrihydrite (Fh) nanoparticles are receiving considerable scientific interest due to their large reactive surface areas, crystalline structures, and nanoparticle morphology. They are of great importance in biogeochemical processes and have the ability to sequester hazardous and toxic substances. Here, the working hypothesis was to entrap fractal-like Fh nanoparticles, with a radius of gyration of 6.2 nm and a primary building block of polydisperse spheres with a radius of 0.8 nm, in a shear-induced multilamellar vesicle (MLV) state using a 40 wt% polyethylene glycol dodecyl ether surfactant.

EXPERIMENTS

Small- and Wide- Angle X-ray scattering revealed the equilibrium state of the non-ionic planar lamellar phase, the Fh dispersion, and their mixture. The MLV state was induced by using a shear flow in a Taylor-Couette geometry of a rheometer.

FINDINGS

The nonionic surfactant initially exhibited a lamellar gel phase with two distinct d-spacings of 11.0 and 9.7 nm, which collapsed into the MLV state under shear flow. The Fh nanoparticles induced bilayer attraction by suppressing lamellar layer undulations, decreasing the d-spacing. These results are helpful in the understanding of the relationship between nanoparticle size and nanoparticle-bilayers interactions and provides insight on Fh encapsulations in a kinetically stable MLVs state.

Flow-Induced Concentration Nonuniformity and Shear Banding in Entangled Polymer Solutions

Recent models have predicted entangled polymer solutions could shear band due to unstable flow-induced demixing. This work provides the first experimental probe of the in situ concentration profile of entangled polymer solutions under shear. At shear rates above a critical value, we show that the concentration and velocity profiles can develop bands, in quantitative agreement with steady-state model predictions. These findings highlight the critical importance of flow-concentration coupling in entangled polymer solutions.

Shear driven vorticity aligned flocs in a suspension of attractive rigid rods

A combination of rheology, optical microscopy and computer simulations was used to investigate the microstructural changes of a semi-dilute suspension of attractive rigid rods in an imposed shear flow. The aim is to understand the relation of the microstructure with the viscoelastic response, and the yielding and flow behaviour in different shear regimes of gels built from rodlike colloids. A semi-dilute suspension of micron sized, rodlike silica particles suspended in 11 M CsCl salt solution was used as a model system for attractive rods' gel. Upon application of steady shear the gel microstructure rearranges in different states and exhibits flow instabilities depending on shear rate, attraction strength, volume fraction and geometrical confinement. At low rod volume fractions, the suspension forms large, vorticity aligned, particle rich flocs that roll in the flow-vorticity plane, an effect that is due to an interplay between hydrodynamic interactions and geometrical confinement as suggested by computer simulations. Experimental data allow the creation of a state diagram, as a function of volume fraction and shear rates, identifying regimes of stable (or unstable) floc formation and of homogeneous gel or broken clusters. The transition is related to dimensionless Mason number, defined as the ratio of shear forces to interparticle attractive force.

Interparticle attraction controls flow heterogeneity in calcite gels

We couple rheometry and ultrasonic velocimetry to study experimentally the flow behavior of gels of colloidal calcite particles dispersed in water, while tuning the strength of the interparticle attraction through physico-chemistry. We unveil, for the first time in a colloidal gel, a direct connection between attractive interactions and the occurrence of shear bands, as well as stress fluctuations.

Flow heterogeneities in supercooled liquids and glasses under shear

Using extensive nonequilibrium molecular dynamics simulations, we investigate a glass-forming binary Lennard-Jones mixture under shear. Both supercooled liquids and glasses are considered. Our focus is on the characterization of inhomogeneous flow patterns such as shear bands that appear as a transient response to the external shear. For the supercooled liquids, we analyze the crossover from Newtonian to non-Newtonian behavior with increasing shear rate γ[over ̇]. Above a critical shear rate γ[over ̇]_{c} where a non-Newtonian response sets in, the transient dynamics are associated with the occurrence of short-lived vertical shear bands, i.e., bands of high mobility that form perpendicular to the flow direction. In the glass states, long-lived horizontal shear bands, i.e., bands of high mobility parallel to the flow direction, are observed in addition to vertical ones. The systems with shear bands are characterized in terms of mobility maps, stress-strain relations, mean-squared displacements, and (local) potential energies. The initial formation of a horizontal shear band provides an efficient stress release, corresponds to a local minimum of the potential energy, and is followed by a slow broadening of the band towards the homogeneously flowing fluid in the steady state. Whether a horizontal or a vertical shear band forms cannot be predicted from the initial undeformed sample. Furthermore, we show that with increasing system size, the probability for the occurrence of horizontal shear bands increases.

Computational study of transient shear banding in soft jammed solids

We have designed three-dimensional numerical simulations of a soft spheres model, with size polidispersity and in athermal conditions, to study the transient shear banding that occurs during yielding of jammed soft solids. We analyze the effects of different types of drag coefficients used in the simulations and compare the results obtained using Lees-Edwards periodic boundary conditions with the case in which the same model solid is confined between two walls. The specific damping mechanism and the different boundary conditions indeed modify the load curves and the velocity profiles in the transient regime. Nevertheless, we find that the presence of a stress overshoot and of a related transient banding phenomenon, for large enough samples, is a robust feature for overdamped systems, where their presence do not depend on the specific drag used and on the different boundary conditions.

Rheology of colloidal and metallic glass formers

Colloidal hard-sphere suspensions are convenient experimental models to understand soft matter, and also by analogy the structural-relaxation behavior of atomic or small-molecular fluids. We discuss this analogy for the flow and deformation behavior close to the glass transition. Based on a mapping of temperature to effective hard-sphere packing, the stress–strain curves of typical bulk metallic glass formers can be quantitatively compared with those of hard-sphere suspensions. Experiments on colloids give access to the microscopic structure under deformation on a single-particle level, providing insight into the yielding mechanisms that are likely also relevant for metallic glasses. We discuss the influence of higher-order angular signals in connection with non-affine particle rearrangements close to yielding. The results are qualitatively explained on the basis of the mode-coupling theory. We further illustrate the analogy of pre-strain dependence of the linear-elastic moduli using data on PS-PNiPAM suspensions.

Emergence and persistence of flow inhomogeneities in the yielding and fluidization of dense soft solids

In three-dimensional computer simulations of model non-Brownian jammed suspensions, we compute the time required to reach homogeneous flow upon yielding, by analyzing stresses and particle packing at different shear rates, with and without confinement. We show that the stress overshoot and persistent shear banding preceding the complete fluidization are controlled by the presence of overconstrained microscopic domains in the initial solids. Such domains, identifiable with icosahedrally packed regions in the model used, allow for stress accumulation during the shear startup. Their structural reorganization under deformation controls the emergence and the persistence of the shear banding.

How clay particulates affect flow cessation and the coiling stability of yield stress-matched cementing suspensions

The remarkable increase in the flow resistance of dense suspensions can hinder 3D-printing processes on account of flow cessation in the extruder, and filament fragility/rupture following deposition. Understanding the nature of rheological changes that occur is critical to manipulate flow conditions or to dose flow modifiers for 3D-printing. Therefore, this paper elucidates the influences of clay particulates on controlling flow cessation and the shape stability of dense cementing suspensions that typically feature poor printability. A rope coiling method was implemented with varying stand-off distances to probe the buckling stability and tendency to fracture of dense suspensions that undergo stretching and bending during deposition. The contributions of flocculation and short-term percolation due to the kinetics of structure formation to deformation rate were deconvoluted using a stepped isostress method. It is shown that the shear stress indicates a divergence with a power-law scaling when the particle volume fraction approaches the jamming limit; φ → φj ≈ φmax. Such a power-law divergence of the shear stress decreases by a factor of 10 with increasing clay dosage. Such behavior in clay-containing suspensions arises from a decrease in the relative packing fraction (φ/φmax) and the formation of fractally-architected aggregates with stronger interparticle interactions, whose uniform arrangement controls flow cessation in the extruder and suspension homogeneity, thereby imparting greater buckling stability. The outcomes offer new insights for assessing/improving the extrudability and printability behavior during slurry-based 3D-printing process.

i-Rheo: determining the linear viscoelastic moduli of colloidal dispersions from step-stress measurements

We report on the application of a Fourier transform-based method, 'i-Rheo', to evaluate the linear viscoelastic moduli of hard-sphere colloidal dispersions, both in the fluid and glass states, from a direct analysis of raw step-stress (creep) experimental data. We corroborate the efficacy of i-Rheo by comparing the outputs of creep tests performed on homogenous complex fluids to conventional dynamic frequency sweeps. A similar approach is adopted for a number of colloidal suspensions over a broad range of volume fractions. For these systems, we test the limits of the method by varying the applied stress across the materials' linear and non-linear viscoelastic regimes, and we show that the best results are achieved for stress values close to the upper limit of the materials' linear viscoelastic regime, where the signal-to-noise ratio is at its highest and the non-linear phenomena have not appeared yet. We record that, the range of accessible frequencies is controlled at the higher end by the relative weight between the inertia of the instrument and the elasticity of the complex material under investigation; whereas, the lowest accessible frequency is dictated by the extent of the materials' linear viscoelastic regime. Nonetheless, despite these constrains, we confirm the effectiveness of i-Rheo for gaining valuable information on the materials' linear viscoelastic properties even from 'creep ringing' data, confirming its potency and general validity as an accurate method for determining the material's rheological behaviour for a variety of complex systems.

Mechanical oscillation accelerating nucleation and nuclei growth in hard-sphere colloidal glass

Crystallization from amorphous solids is generally caused by activating phonons in a wide frequency range during heat treatment. In contrast, the activation of phonons in a narrow frequency range using ultrasonic treatment also causes crystallization below the glass transition temperature. These behaviors indicate that crystallization is related to the atomic motion in the glass state, and it is suggested that the activation of specific atomic motion can cause crystallization without increasing temperature. In this study, we observe nucleation and nuclei growth caused by mechanical oscillation in a hard-sphere colloidal glass and evaluate the effect of mechanical oscillation on the structural evolution in the early stage of the crystallization. Oscillation between 5 and 100 Hz is applied to the colloidal glass, and it is observed that the nucleation rate increases under the 70 Hz oscillation, resulting in formation of stable nuclei in a short amount of time. The nuclei growth is also accelerated by the 70 Hz oscillation, whereas increases in the nucleation rate and nuclei growth were not observed at other frequencies. Finally, activation of the diffusion-based rattling of particles by caging is considered as a possible mechanism of the observations.

Local elasticity in nonlinear rheology of interacting colloidal glasses revealed by neutron scattering and rheometry

The flow of colloidal suspensions is ubiquitous in nature and industry. Colloidal suspensions exhibit a wide range of rheological behavior, which should be closely related to the microscopic structure of the systems. With in situ small-angle neutron scattering complemented by rheological measurements, we investigated the deformation behavior of a charge-stabilized colloidal glass at particle level undergoing steady shear. A short-lived, localized elastic response at particle level, termed as the transient elasticity zone (TEZ), was identified from the neutron spectra. The existence of the TEZ, which could be promoted by the electrostatic interparticle potential, is a signature of deformation heterogeneity: the body of fluids under shear behaves like an elastic solid within the spatial range of the TEZ but like fluid outside the TEZ. The size of the TEZ shrinks as the shear rate increases in the shear thinning region, which shows that the shear thinning is accompanied by a diminishing deformation heterogeneity. More interestingly, the TEZ is found to be the structural unit that provides the resistance to the imposed shear, as evidenced by the quantitative agreement between the local elastic stress sustained by the TEZ and the macroscopic stress from rheological measurements at low and moderate shear rates. Our findings provide an understanding on the nonlinear rheology of interacting colloidal glasses from a micro-mechanical view.

Geological implication of grain-size segregation in dense granular matter

To the current common belief, grain size segregation in granular matter requires sufficient porosity. Therefore, grain size segregation found in a natural fault gouge could imply elevated fluid pressure and the reduced normal stress on fault, possibly caused by the frictional heat during an earthquake. To clarify whether fluidization is essential to grain size segregation, we conduct numerical simulation on a simple model of fault gouge in a plane shear geometry under constant volume condition: the volume fraction is fixed at 0.6, at which the granular system possesses yield stress. We observe apparent grain size segregation at this volume fraction, meaning that grain size segregation alone does not imply fluidization of granular matter. We also show that segregation is driven by the nonlinear velocity profile, and that the gravity is not essential to segregation. The physical condition tested here may be relevant to earthquake faults: the normal stress of 1 MPa, the sliding velocity of 1 m s^-1, and the duration of 0.1 s.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Dynamic Vorticity Banding in Discontinuously Shear Thickening Suspensions

It has recently been argued that steady-state vorticity bands cannot arise in shear thickening suspensions because the normal stress imbalance across the interface between the bands will set up particle migrations. In this Letter, we develop a simple continuum model that couples shear thickening to particle migration. We show by linear stability analysis that homogeneous flow is unstable towards vorticity banding, as expected, in the regime of negative constitutive slope. In full nonlinear computations, we show, however, that the resulting vorticity bands are unsteady, with spatiotemporal patterns governed by stress-concentration coupling. We furthermore show that these dynamical bands also arise in direct particle simulations, in good agreement with the continuum model.

Non-local stresses in highly non-uniformly flowing suspensions: The shear-curvature viscosity

For highly non-uniformly flowing fluids, there are contributions to the stress related to spatial variations of the shear rate, which are commonly referred to as non-local stresses. The standard expression for the shear stress, which states that the shear stress is proportional to the shear rate, is based on a formal expansion of the stress tensor with respect to spatial gradients in the flow velocity up to leading order. Such a leading order expansion is not able to describe fluids with very rapid spatial variations of the shear rate, like in micro-fluidics devices and in shear-banding suspensions. Spatial derivatives of the shear rate then significantly contribute to the stress. Such non-local stresses have so far been introduced on a phenomenological level. In particular, a formal gradient expansion of the stress tensor beyond the above mentioned leading order contribution leads to a phenomenological formulation of non-local stresses in terms of the so-called "shear-curvature viscosity". We derive an expression for the shear-curvature viscosity for dilute suspensions of spherical colloids and propose an effective-medium approach to extend this result to concentrated suspensions. The validity of the effective-medium prediction is confirmed by Brownian dynamics simulations on highly non-uniformly flowing fluids.

Shear-density coupling for a compressible single-component yield-stress fluid

Flow behavior of a single-component yield stress fluid is addressed on the hydrodynamic level. A basic ingredient of the model is a coupling between fluctuations of density and velocity gradient via a Herschel-Bulkley-type constitutive model. Focusing on the limit of low shear rates and high densities, the model approximates well-but is not limited to-gently sheared hard sphere colloidal glasses, where solvent effects are negligible. A detailed analysis of the linearized hydrodynamic equations for fluctuations and the resulting cubic dispersion relation reveals the existence of a range of densities and shear rates with growing flow heterogeneity. In this regime, after an initial transient, the velocity and density fields monotonically reach a spatially inhomogeneous stationary profile, where regions of high shear rate and low density coexist with regions of low shear rate and high density. The steady state is thus maintained by a competition between shear-induced enhancement of density inhomogeneities and relaxation via overdamped sound waves. An analysis of the mechanical equilibrium condition provides a criterion for the existence of steady state solutions. The dynamical evolution of the system is discussed in detail for various boundary conditions, imposing either a constant velocity, shear rate, or stress at the walls.

Yielding of a model glass former: An interpretation with an effective system of icosahedra

We consider the yielding under simple shear of a binary Lennard-Jones glass former whose super-Arrhenius dynamics are correlated with the formation of icosahedral structures. We recast this glass former as an effective system of icosahedra [Pinney et al., J. Chem. Phys. 143, 244507 (2015)JCPSA60021-960610.1063/1.4938424]. Looking at the small-strain region of sheared simulations, we observe that shear rates affect the shear localization behavior particularly at temperatures below the glass transition as defined with a fit to the Vogel-Fulcher-Tamman equation. At higher temperature, shear localization starts immediately on shearing for all shear rates. At lower temperatures, faster shear rates can result in a delayed start in shear localization, which begins close to the yield stress. Building from a previous work which considered steady-state shear [Pinney et al., J. Chem. Phys. 143, 244507 (2015)JCPSA60021-960610.1063/1.4938424], we interpret the response to shear and the shear localization in terms of a local effective temperature with our system of icosahedra. We find that the effective temperatures of the regions undergoing shear localization increase significantly with increasing strain (before reaching a steady-state plateau).

Shear-Induced Heterogeneity in Associating Polymer Gels: Role of Network Structure and Dilatancy

We study associating polymer gels under steady shear using Brownian dynamics simulation to explore the interplay between the network structure, dynamics, and rheology. For a wide range of flow rates, we observe the formation of shear bands with a pronounced difference in shear rate, concentration, and structure. A striking increase in the polymer pressure in the gradient direction with shear, along with the inherently large compressibility of the gels, is shown to be a crucial factor in destabilizing homogeneous flow through shear-gradient concentration coupling. We find that shear has only a modest influence on the degree of association, but induces marked spatial heterogeneity in the network connectivity. We attribute the increase in the polymer pressure (and polymer mobility) to this structural reorganization.

Accelerated crystallization of colloidal glass by mechanical oscillation

Crystallization of a hard-sphere colloidal glass by mechanical oscillation is investigated, and accelerated crystallization is found at a specific frequency. The crystallization frequency increases as attractive force between particles increases, indicating that interparticle interaction affects the crystallization frequency. Time scale of the mechanical oscillation is different from that of the slow relaxation, and notable relationship with the low-frequency mode is not observed. The experimental results are not explained by the previously proposed model for crystallization by oscillatory shear. Conversely, we speculate that activations of the fast relaxation and particle motion in crystal-like clusters are possible causes of the observations.

Shear banding, discontinuous shear thickening, and rheological phase transitions in athermally sheared frictionless disks

We report on numerical simulations of simple models of athermal, bidisperse, soft-core, massive disks in two dimensions, as a function of packing fraction ϕ, inelasticity of collisions as measured by a parameter Q, and applied uniform shear strain rate γ[over ̇]. Our particles have contact interactions consisting of normally directed elastic repulsion and viscous dissipation, as well as tangentially directed viscous dissipation, but no interparticle Coulombic friction. Mapping the phase diagram in the (ϕ,Q) plane for small γ[over ̇], we find a sharp first-order rheological phase transition from a region with Bagnoldian rheology to a region with Newtonian rheology, and show that the system is always Newtonian at jamming. We consider the rotational motion of particles and demonstrate the crucial importance that the coupling between rotational and translational degrees of freedom has on the phase structure at small Q (strongly inelastic collisions). At small Q, we show that, upon increasing γ[over ̇], the sharp Bagnoldian-to-Newtonian transition becomes a coexistence region of finite width in the (ϕ,γ[over ̇]) plane, with coexisting Bagnoldian and Newtonian shear bands. Crossing this coexistence region by increasing γ[over ̇] at fixed ϕ, we find that discontinuous shear thickening can result if γ[over ̇] is varied too rapidly for the system to relax to the shear-banded steady state corresponding to the instantaneous value of γ[over ̇].

Onset of shear thinning in glassy liquids: Shear-induced small reduction of effective density

We propose a simple mechanism for describing the onset of shear thinning in a high-density glassy liquid. In a shear flow, along the compression axis, the overlap between neighboring particles is more enhanced than that at equilibrium, meaning that the "effective" size is reduced along this axis. On the other hand, along the extension axis perpendicular to the compression axis, the average structural configurations are stretched, but it does not indicate the expansion of the "effective" size itself. This asymmetric shear flow effect for particles results in a small reduction of the "effective" density. Because, in glass-forming liquids, the structural relaxation time τ_{α} strongly depends on the density ρ, even a very small reduction of the effective density should lead to a significant decrease of the relaxation time under shear flow. We predict that the crossover shear rate from Newtonian to non-Newtonian flow behaviors is given by γ[over ̇]_{c}=[ρ(∂τ_{α}/∂ρ)]^{-1}, which can be much smaller than 1/τ_{α} near the glass transition point. It is shown that this prediction is consistent with the results of molecular dynamics simulations.

Structural and cooperative length scales in polymer gels

Understanding the relationship between the material structural details, the geometrical confining constraints, the local dynamical events and the global rheological response is at the core of present investigations on complex fluid properties. In the present article, this problem is addressed on a model yield stress fluid made of highly entangled polymer gels of Carbopol which follows at the macroscopic scale the well-known Herschel-Bulkley rheological law. First, performing local rheology measurements up to high shear rates ([Formula: see text] s^-1)and under confinement, we evidence unambiguously the breakdown of bulk rheology associated with cooperative processes under flow. Moreover, we show that these behaviors are fully captured with a unique cooperativity length [Formula: see text] over the whole range of experimental conditions. Second, we introduce an original optical microscopy method to access structural properties of the entangled polymer gel in the direct space. Performing image correlation spectroscopy of fluorophore-loaded gels, the characteristic size D of carbopol gels microstructure is determined as a function of preparation protocol. Combining both dynamical and structural information shows that the measured cooperative length [Formula: see text] corresponds to 2-5 times the underlying structural size D, thus providing a strong grounding to the "Shear Transformation Zones" modeling approach.

A stress-controlled shear cell for small-angle light scattering and microscopy

We develop and test a stress-controlled, parallel plates shear cell that can be coupled to an optical microscope or a small angle light scattering setup, for simultaneous investigation of the rheological response and the microscopic structure of soft materials under an imposed shear stress. In order to minimize friction, the cell is based on an air bearing linear stage, the stress is applied through a contactless magnetic actuator, and the strain is measured through optical sensors. We discuss the contributions of inertia and of the small residual friction to the measured signal and demonstrate the performance of our device in both oscillating and step stress experiments on a variety of viscoelastic materials.

Yielding of glass under shear: A directed percolation transition precedes shear-band formation

Under external mechanical loading, glassy materials, ranging from soft matter systems to metallic alloys, often respond via formation of inhomogeneous flow patterns, during yielding. These inhomogeneities can be precursors to catastrophic failure, implying that a better understanding of their underlying mechanisms could lead to the design of smarter materials. Here, extensive molecular dynamics simulations are used to reveal the emergence of heterogeneous dynamics in a binary Lennard-Jones glass, subjected to a constant strain rate. At a critical strain, this system exhibits for all considered strain rates a transition towards the formation of a percolating cluster of mobile regions. We give evidence that this transition belongs to the universality class of directed percolation. Only at low shear rates, the percolating cluster subsequently evolves into a transient (but long-lived) shear band with a diffusive growth of its width. Finally, the steady state with a homogeneous flow pattern is reached. In the steady state, percolation transitions also do occur constantly, albeit over smaller strain intervals, to maintain the stationary plastic flow in the system.

Two-dimensional magnetic colloids under shear

Complex rheological properties of soft disordered solids, such as colloidal gels or glasses, inspire a range of novel applications. However, the microscopic mechanisms of their response to mechanical loading are not well understood. Here, we elucidate some aspects of these mechanisms by studying a versatile model system, i.e. two-dimensional superparamagnetic colloids in a precessing magnetic field, whose structure can be tuned from a hexagonal crystal to a disordered gel network by varying the external field opening angle θ. We perform Langevin dynamics simulations subjecting these structures to a constant shear rate and observe three qualitatively different types of material response. In hexagonal crystals (θ = 0°), at a sufficiently low shear rate, plastic flow occurs via successive stress drops at which the stress releases due to the formation of dislocation defects. The gel network at θ = 48°, on the contrary, via bond rearrangement and transient shear banding evolves into a homogeneously stretched network at large strains. The latter structure remains metastable after switching off of the shear. At θ = 50°, the external shear makes the system unstable against phase separation and causes a failure of the network structure leading to the formation of hexagonal close packed clusters interconnected by particle chains. At a microcopic level, our simulations provide insight into some of the mechanisms by which strain localization as well as material failure occur in a simple gel-like network. Furthermore, we demonstrate that new stretched network structures can be generated by the application of shear.

Wall slip across the jamming transition of soft thermoresponsive particles

Flows of suspensions are often affected by wall slip, that is, the fluid velocity v(f) in the vicinity of a boundary differs from the wall velocity v(w) due to the presence of a lubrication layer. While the slip velocity v(s)=|v(f)-v(w)| robustly scales linearly with the stress σ at the wall in dilute suspensions, there is no consensus regarding denser suspensions that are sheared in the bulk, for which slip velocities have been reported to scale as a v(s)∝σ(p) with exponents p inconsistently ranging between 0 and 2. Here we focus on a suspension of soft thermoresponsive particles and show that v(s)) actually scales as a power law of the viscous stress σ-σ(c), where σ(c) denotes the yield stress of the bulk material. By tuning the temperature across the jamming transition, we further demonstrate that this scaling holds true over a large range of packing fractions ϕ on both sides of the jamming point and that the exponent p increases continuously with ϕ, from p=1 in the case of dilute suspensions to p=2 for jammed assemblies. These results allow us to successfully revisit inconsistent data from the literature and pave the way for a continuous description of wall slip above and below jamming.

Convective Cage Release in Model Colloidal Glasses

The mechanism of flow in glassy materials is interrogated using mechanical spectroscopy applied to model nearly hard sphere colloidal glasses during flow. Superimposing a small amplitude oscillatory motion orthogonal onto steady shear flow makes it possible to directly evaluate the effect of a steady state flow on the out-of-cage (α) relaxation as well as the in-cage motions. To this end, the crossover frequency deduced from the viscoelastic spectra is used as a direct measure of the inverse microstructural relaxation time, during flow. The latter is found to scale linearly with the rate of deformation. The microscopic mechanism of flow can then be identified as a convective cage release. Further insights are provided when the viscoelastic spectra at different shear rates are shifted to scale the alpha relaxation and produce a strain rate-orthogonal frequency superposition, the colloidal analogue of time temperature superposition in polymers with the flow strength playing the role of temperature. Whereas the scaling works well for the α relaxation, deviations are observed both at low and high frequencies. Brownian dynamics simulations point to the origins of these deviations; at high frequencies these are due to the deformation of the cages which slows down the short-time diffusion, while at low frequency, deviations are most probably caused by some mild hydroclustering.

Colloidal binary mixtures at fluid-fluid interfaces under steady shear: structural, dynamical and mechanical response

We experimentally study the link between structure, dynamics and mechanical response of two-dimensional (2D) binary mixtures of colloidal microparticles spread at water/oil interfaces. The particles are driven into steady shear by a microdisk forced to rotate at a controlled angular velocity. The flow causes particles to layer into alternating concentric rings of small and big colloids. The formation of such layers is linked to the local, position-dependent shear rate, which triggers two distinct dynamical regimes: particles either move continuously ("Flowing") close to the microdisk, or exhibit intermittent "Hopping" between local energy minima farther away. The shear-rate-dependent surface viscosity of the monolayers can be extracted from a local interfacial stress balance, giving "macroscopic" flow curves whose behavior corresponds to the distinct microscopic regimes of particle motion. Hopping regions reveal a higher resistance to flow compared to the flowing regions, where spatial organization into layers reduces dissipation.

Understanding the approximations of mode-coupling theory for sheared steady states of colloids

The lack of clarity of various mode-coupling theory (MCT) approximations, even in equilibrium, makes it hard to understand the relation between various MCT approaches for sheared steady states as well as their regime of validity. Here we try to understand these approximations indirectly by deriving the MCT equations through two different approaches for a colloidal system under shear, first through a microscopic approach, as suggested by Zaccarelli et al., and second through fluctuating hydrodynamics, where the approximations used in the derivation are quite clear. The qualitative similarity of our theory with a number of existing theories show that linear response theory might play a role in various approximations employed in deriving those theories and one needs to be careful while applying them for systems arbitrarily far away from equilibrium, such as a granular system or when shear is very strong. As a by-product of our calculation, we obtain the extension of the Yvon-Born-Green (YBG) equation for a sheared system and under the assumption of random-phase approximation, the YBG equation yields the distorted structure factor that was earlier obtained through different approaches.

Shear-induced breaking of cages in colloidal glasses: Scattering experiments and mode coupling theory

We employ x-ray scattering on sheared colloidal suspensions and mode coupling theory to study structure factor distortions of glass-forming systems under shear. We find a transition from quadrupolar elastic distortion at small strains to quadrupolar and hexadecupolar modes in the stationary state. The latter are interpreted as signatures of plastic rearrangements in homogeneous, thermalized systems. From their transient evolution with strain, we identify characteristic strain and length-scale values where these plastic rearrangements dominate. This characteristic strain coincides with the maximum of the shear stress versus strain curve, indicating the proliferation of plastic flow. The hexadecupolar modes dominate at the wavevector of the principal peak of the equilibrium structure factor that is related to the cage-effect in mode coupling theory. We hence identify the structural signature of plastic flow of glasses.

Fracture of Jammed Colloidal Suspensions

Concentrated colloidal suspensions display dramatic rises in viscosity, leading to jamming and granulation, with increasing shear rate. It has been proposed that these effects result from inter particle friction, as lubrication forces are overcome. This suggests the jamming of concentrated colloidal suspensions should exhibit some shared phenomenology with macroscopic granular systems where friction leads to two different types of jammed state. Here we show that transient rheological measurements can be used to probe the processes of granulation in concentrated colloidal suspensions. Our results support the idea that frictional contacts are created between jammed particles. The jamming behaviour displays two qualitatively different regimes separated by a critical strain rate with qualitatively different types of fracture/break up behaviour. In the lower strain rate regime, it is found that vibrations can be used to control jamming and granulation, resulting in a flowable fluid.

S-shaped flow curves of shear thickening suspensions: direct observation of frictional rheology

We study the rheological behavior of concentrated granular suspensions of simple spherical particles. Under controlled stress, the system exhibits an S-shaped flow curve (stress vs shear rate) with a negative slope in between the low-viscosity Newtonian regime and the shear thickened regime. Under controlled shear rate, a discontinuous transition between the two states is observed. Stress visualization experiments with a fluorescent probe suggest that friction is at the origin of shear thickening. Stress visualization shows that the stress in the system remains homogeneous (no shear banding) if a stress is imposed that is intermediate between the high- and low-stress branches. The S-shaped shear thickening is then due to the discontinuous formation of a frictional force network between particles upon increasing the stress.

Tuning colloidal gels by shear

Using a powerful combination of experiments and simulations we demonstrate how the microstructure and its time evolution are linked with mechanical properties in a frustrated, out-of-equilibrium, particle gel under shear. An intermediate volume fraction colloid-polymer gel is used as a model system, allowing quantification of the interplay between interparticle attractions and shear forces. Rheometry, confocal microscopy and Brownian dynamics reveal that high shear rates, fully breaking the structure, lead after shear cessation to more homogeneous and stronger gels, whereas preshear at low rates creates largely heterogeneous weaker gels with reduced elasticity. We find that in comparison, thermal quenching cannot produce structural inhomogeneities under shear. We argue that external shear has strong implications on routes towards metastable equilibrium, and therefore gelation scenarios. Moreover, these results have strong implications for material design and industrial applications, such as mixing, processing and transport protocols coupled to the properties of the final material.

Shear banding in drying films of colloidal nanoparticles

Drying suspensions of colloidal nanoparticles exhibit a variety of interesting strain release mechanisms during film formation. These result in the selection of characteristic length scales during failure processes such as cracking and subsequent delamination. A wide range of materials (e.g., bulk metallic glasses) release strain through plastic deformations which occur in a narrow band of material known as a shear band. Here we show that drying colloidal films also exhibit shear banding. Bands are observed to form a small distance behind the drying front and then to propagate rapidly at ∼45° to the direction of drying. It is shown that the spacing of the bands depends on salt concentration and the evaporation rate of the colloidal suspension. These combined observations suggest that there is a critical shear rate (related to the film yield stress) which controls the ratio of bandwidth to band spacing. Local deformations were measured in the early stages of drying using fluorescent tracer particles. The measurements were used to show that the existence of shear bands is linked to the compaction of particles perpendicular to the drying front. The spacing of shear bands was also found to be strongly correlated with the characteristic length scale of the compaction process. These combined studies elucidate the role of plastic deformation during pattern formation in drying films of colloidal nanoparticles.

Classical XY model with conserved angular momentum is an archetypal non-Newtonian fluid

We find that the classical one-dimensional XY model, with angular-momentum-conserving Langevin dynamics, mimics the non-Newtonian flow regimes characteristic of soft matter when subjected to counterrotating boundaries. An elaborate steady-state phase diagram has continuous and first-order transitions between states of uniform flow, shear-banding, solid-fluid coexistence and slip planes. Results of numerical studies and a concise mean-field constitutive relation offer a paradigm for diverse nonequilibrium complex fluids.

Manipulating shear-induced non-equilibrium transitions in colloidal films by feedback control

Using Brownian Dynamics (BD) simulations we investigate non-equilibrium transitions of sheared colloidal films under controlled shear stress σxz. In our approach the shear rate [small gamma, Greek, dot above] is a dynamical variable, which relaxes on a time scale τc such that the instantaneous, configuration-dependent stress σxz(t) approaches a pre-imposed value. Investigating the dynamics under this "feedback-control" scheme we find unique behavior in regions where the flow curve σxz([small gamma, Greek, dot above]) of the uncontrolled system is monotonic. However, in non-monotonic regions our method allows to select between dynamical states characterized by different in-plane structure and viscosities. Indeed, the final state strongly depends on τc relative to an intrinsic relaxation time of the uncontrolled system. The critical values of τc are estimated on the basis of a simple model.

Flow instability due to coupling of shear-gradients with concentration: non-uniform flow of (hard-sphere) glasses

Flow-induced instabilities that lead to non-uniform stationary flow profiles have been observed in many different soft-matter systems. Two types of instabilities that lead to banded stationary states have been identified, which are commonly referred to as gradient- and vorticity-banding. The molecular origin of these instabilities is reasonably well understood. A third type of instability that has been proposed phenomenologically [Europhys. Lett., 1986, 2, 129 and Phys. Rev. E, 1995, 52, 4009] is largely unexplored. Essential to this "Shear-gradient Concentration Coupling" (SCC-) instability is a mass flux that is induced by spatial gradients of the shear rate. A possible reason that this instability has essentially been ignored is that the molecular origin of the postulated mass flux is not clear, and no explicit expressions for the shear-rate and concentration dependence of the corresponding transport coefficient exist. It is therefore not yet known what types of flow velocity- and concentration-profiles this instability gives rise to. In this paper, an expression for the transport coefficient corresponding to the shear-gradient induced mass flux is derived in terms of the shear-rate dependent pair-correlation function, and Brownian dynamics simulations for hard-spheres are presented that specify the shear-rate and concentration dependence of the pair-correlation function. This allows to explicitly formulate the coupled advection-diffusion equation and an equation of motion for the suspension flow velocity. The inclusion of a non-local contribution to the stress turns out to be essential to describe the SCC-banding transition. The coupled equations of motion are solved numerically, and flow- and concentration-profiles are discussed. It is shown that the SCC-instability occurs within the glass state at sufficiently small shear rates, leading to a banded flow-profile where one of the bands is non-flowing.