Time-Dependent Flow in Arrested States —Transient Behaviour

Time-Dependent Flow in Arrested States —
Transient Behaviour

K J Mutch M Laurati C P Amann M Fuchs and S U Egelhaaf
(a) Condensed Matter Physics Laboratory, Heinrich-Heine University,
40225 Düsseldorf, Germany

(b) Fachbereich Physik, University of Konstanz, 78457 Konstanz, Germany
* kevin.mutch@uni-duesseldorf.de

The transient behaviour of highly concentrated colloidal liquids and dynamically arrested states (glasses) under time-dependent shear is reviewed. This includes both theoretical and experimental studies and comprises the macroscopic rheological behaviour as well as changes in the structure and dynamics on a microscopic individual-particle level. The microscopic and macroscopic levels of the systems are linked by a comprehensive theoretical framework which is exploited to quantitatively describe these systems while they are subjected to an arbitrary flow history. Within this framework, theoretical predictions are compared to experimental data, which were gathered by rheology and confocal microscopy experiments, and display consistent results. Particular emphasis is given to (i) switch-on of shear flow during which the system can liquify, (ii) switch-off of shear flow which might still leave residual stresses in the system, and (iii) large amplitude oscillatory shearing. The competition between timescales and the dependence on flow history leads to novel features in both the rheological response and the microscopic structure and dynamics.

1 Introduction

The ubiquity of colloidal dispersions, coupled with their ability to act as suitable models for atomic and molecular systems [1], make them both academically and commercially very relevant. One of the primary advantages of studying colloidal systems is their size: they are often large enough to be seen with a microscope and slow enough to be followed. Moreover, the typical shear modulus of a colloidal dispersion is orders of magnitude lower than equivalent atomic and molecular systems [2, 3], rendering the flow properties of such systems experimentally accessible. This has relevance to industrial processes, where materials are often dispersed as colloidal particles in order to pour or transport them, dispersion paints being just one familiar example. It is thus not surprising that the rheology of dispersions has been the focus of longstanding research and a number of approaches have been developed in this field, some of which were recently reviewed [4].

For the case of hard sphere colloids, the thermodynamic control parameter is the volume fraction of particles, . In a monodisperse system, the colloids exhibit a disordered fluid phase until freezing at , whereby a first-order phase transition to a crystal phase of occurs [5, 6, 7]. An increase in the polydispersity of the colloidal dispersion leads to a frustration of the crystalline ordering and suppression of crystallisation for polydispersities greater than [8]. Instead, a non-equilibrium glassy state is observed for  [9, 6], at which point the dynamics become arrested.

Such concentrated colloidal systems can be conveniently characterised by dynamic light scattering (DLS) and confocal microscopy experiments. Observation of the intermediate scattering function (ISF) from DLS reveals a non-decaying component upon increasing beyond the glass transition, suggesting a dynamical slowing down [10, 11]. The mode-coupling theory (MCT) [12, 13] has been shown to accurately describe these experimentally determined ISF [14, 15]. Mean-squared displacements (MSD) can be derived from the ISF and can also be measured by confocal microscopy; the slowing down in the dynamics as increases is evidenced by an emerging plateau in the MSD, through which a localisation length related to particle caging can be determined. In the concentrated disordered fluid state, particles are temporarily caged by their neighbours but can in time escape these cages and display diffusion at long times. As the concentration increases towards and above the glass transition, the cages tighten and the particles become immobilised; their only motion is local movement within the cages. The restricted mobility has consequences for their response to mechanical stress, such as reversible elastic deformation and resistance to flow at small applied stresses followed by yielding and irreversible deformation at large applied stresses.

The low yield stress of colloidal glasses, i.e. their mechanical weakness (or softness), imparts interesting and varied behaviour, particularly under conditions of flow. A long observed phenomenon in the sheared steady state is that of shear thinning, i.e. the decrease of the shear viscosity with increasing rate of shear [16, 17, 18, 19, 20, 21]. In order to understand such macroscopic behaviour, it is necessary to draw a link to the microstructure of the system [22], which can be achieved through modern scattering and imaging techniques. For example, the application of light scattering echo [23, 24, 25, 26] and more recently confocal microscopy [27, 28] has allowed the origins of shear thinning to be better understood. From microscopy experiments, it was seen that the inverse relaxation time of a sheared glass was closely linked to the local rather than global shear rate. The converse case of shear thickening can also arise at high shear rates due to cluster formation induced by increased hydrodynamic interactions [29, 30, 31, 32, 33].

These studies of colloidal systems under flow, i.e. in the steady state of shear, effectively require “pre-yielding” and a certain strain to have been reached. In order to achieve this however, first the linear regime must be traversed, in which the stress is proportional to strain , followed by the transient regime. It is in this transient regime of intermediate strains that some particularly interesting phenomena can occur. The transport coefficients measured in quiescent equilibrium or linear response are known to differ greatly from those in steady shear states and hence the microscopic properties of a system exposed to a sudden change in externally applied shear are of particular interest. Moreover, the transient regime encompasses the shear melting of colloidal glasses and could potentially provide greater understanding of glassy states at rest.

First insight to the transient response can be gleaned by measuring the evolution of the stress, , under conditions of constant shear rate, . For certain colloidal systems, the stress evolves non-monotonically [34, 35, 36, 37]. Upon start-up of shear, the stress grows almost linearly with time before reaching a maximum value and subsequently decaying to the steady state asymptotic value. The intermediate maximum is referred to as the stress overshoot and has also been observed in other materials displaying a yield stress [38, 39, 40, 41, 42, 43, 44]. The ability of a system to yield is not a pre-requisite to display an overshoot however, as systems of entangled polymers and wormlike micelles also exhibit such behaviour [45, 46, 47, 48, 49, 50, 51, 52, 53, 54]; it has been proposed that these systems display yielding like behaviour due to the disentangling of polymer/micelle chains [54, 55]. Conversely, a microscopic understanding of yielding and the presence of a stress overshoot in dense colloidal systems is not so clear. Hence, a discussion of the rheology, structure and dynamics of colloidal suspensions during the transient regimes of both start-up and cessation of flow, along with the periodic switching on and off of shear in the form of large-amplitude oscillatory shear (LAOS), will form the focus of this review.

2 Theory

2.1 Overview of theoretical approaches

Theoretical treatment of the rheology of dense colloidal suspensions generally falls into one of three main categories, each of which will be briefly summarised.

Firstly, phenomenological approaches are based on kinetic rules or dynamical equations for coarse-grained fields. These approaches traditionally start from continuum hydrodynamics by formulating constitutive equations for the stress-strain relations and closing the Navier-Stokes equations [56]. This led to, e.g., the ‘principle of material objectivity’, which constrains constitutive equations to be invariant under arbitrary (including time-dependent) solid body rotations. Nevertheless, a detailed microscopic picture of local structural deformations under flow has not been provided.

In the mesoscopic regime, other phenomenological approaches like shear transformation zone (STZ) theory have been applied to linear dispersion moduli [57]. STZ considers kinetic rate equations for density fields of flow defects which, when present in small numbers, could be observed in colloidal glasses using confocal microscopy [58]. The existence of flow defects at densities below the glass transition is, however, still unclear. Soft glassy rheology (SGR) is another phenomenological approach, whereby a set of elastic elements with distributed yield thresholds are considered [59]. They undergo activated dynamics governed by a noise temperature, which differs from the thermodynamic one. The microscopic basis of this model for colloids has not yet been established, and new results indicate that SGR is more adapted to jamming or athermal systems such as granular materials [60, 61]. They exhibit different shear-thinning behaviour than dense colloidal dispersions, which is discussed in the present review.

Other microscopic approaches start from the underlying stochastic processes, where the colloids display Brownian motion in an incompressible solvent. A multi-particle stochastic equation (Smoluchowski equation) captures the dynamics and rheology of colloidal dispersions for times longer than the diffusion time , i.e. the time it takes an isolated colloid to diffuse its own radius  [62]. Here, is the Stokes-Einstein-Sutherland diffusion coefficient, which can be determined from the solvent viscosity .

The second class of approaches specifically addresses the low density limit, where consideration of the two-particle distribution function suffices [32, 33, 63, 64]. The resulting theory for the relative density distributions, stress and flow fields can be solved rigorously, providing insight into shear-thinning and shear-thickening phenomena at low densities. The complete inclusion of hydrodynamic interactions has been achieved [63]; these interactions propagate through the solvent instantaneously, causing both lubrication forces in the near-field and long-ranged force-velocity correlations in the far-field. For large flow rates, the distortion of the probability of particles to approach each other varies on two very different length scales: for large separations, flow randomizes the correlations between the particle-pair, whilst for small separations, diffusion and flow compete in a narrow boundary layer at contact. For flows at a steady flow rate , the competition is characterized by the bare Péclet number , i.e. the ratio between the rates of external flow and internal diffusion. Asymptotically, the boundary layer obtains a width and height that increases with . In the strongly varying local density field, local stress fields arise and determine the macroscopic rheological properties [64, 63]. For colloidal hard spheres, the density increase is measured by the contact value of the pair-correlation function. This aspect is the basis for the extension to higher packing fractions [63]. The binary-correlation results obtained at low densities are scaled into the semi-dilute regime using the packing fraction dependent contact value. This scaling has also been suggested to hold close to random-close packing, where the contact value diverges [65]. Because this scaling focuses on the local packing, which is connected to motion characterised by the short-time self-diffusivity, the slow structural relaxation is neglected. The latter are connected to the long-time relaxation in the supercooled fluid. The approach [63] therefore does not apply to the region around the glass transition, where cooperative structural rearrangements reaching beyond the boundary layer dominate.

The third and final group of theories are those which directly address concentrated dispersions and their density fluctuations through a microscopic approach. For example, time dependent density functional theories (DDFT) have been intensively pursued [66]. Whilst some DDFT approaches contain an instability leading to glassy arrest, the incorporation of structural memory and slow dynamics remains unresolved. If the system at low enough density escapes the DDFT arrest-phenomenon, its long-time diffusivity returns to the value at infinite dilution. At long times, the supercooled fluid has no memory of the incipient glassy arrest in DDFT [67], which is clearly at odds with dynamic light scattering observations [68]. Another route to the nonlinear rheology of dense dispersions builds on the mode coupling theory (MCT) of the glass transition [12, 13] and extends it to nonequilibrium phenomena via the integration through transients (ITT) approach [69, 70]. The MCT-ITT approach is focused on the slow structural relaxation phenomena under weak flows which include shear-thinning as the main effect in complex fluids. Recent developments in the generalisation of MCT-ITT to incompressible and homogeneous but otherwise arbitrary weak flows in dense dispersions [71, 72], where weak entails , are reviewed in the following section.

2.2 Constitutive relations of microscopic MCT-ITT

The MCT-ITT approach was developed in order to describe steadily sheared dispersions where internal structural relaxation competes with an external driving force [73, 74]. It was recently generalized to arbitrary incompressible and homogeneous weak flows [71, 72], in order to incorporate transient phenomena and general flow geometries. Shear thinning remains the dominating phenomenon captured in the MCT-ITT equations. It results from the flow advection of density fluctuations, which are captured in wavevector space. Consideration of spatial Fourier-modes has proven to be an efficient method to describe structural relaxation, presumably because glassy arrest results from a nonlinear scattering process where plane wave density fluctuations at wavevectors and scatter and excite a fluctuation at . The wavevectors become time-dependent due to flow advection because the affine particle motion under flow tilts or distorts the plane wave modes. Describing the homogeneous flow with a velocity gradient tensor , the advected wavevector


enters the theory. Here, the time dependence of the flow between the times and is given by the (time-ordered) exponential that also describes the deformation in real space:


The tensor is the deformation gradient tensor [56]. The principle of material objectivity, which requires that the nonlinear rheology be independent of solid-body rotations of the sample, is verified in MCT-ITT because at each appearance of it can be combined into the finger tensor , which is manifestly invariant [71, 72].

In order to establish a constitutive relation which links stress and strain, in MCT-ITT first a formally exact generalized Green-Kubo relation for the macroscopic stress tensor is derived. It can be considered a formally exact constitutive equation under homogeneous flow


Here, is the fluctuating stress tensor. The averaged stress tensor is a nonlinear functional of the velocity gradient tensor as appears in the adjoint Smoluchowski operator   [71, 72].

Building on the insights gained through MCT into the cage-effect in supercooled liquids [13], approximations are formulated which enable evaluation of the formally exact Green-Kubo relation. The first approximation expresses the stress fluctuations in Eq. (3) in terms of pair density fluctuations. Stresses are supposed to arise from density fluctuations and to become slow because of the slow structural rearrangements which are captured in the ensemble averaged transient density correlators


The forward-advected wavevector, Eq. (1), appears naturally in Eq. (4). The correlator arises in the approximated stress, Eq. (3), after a Wick-like decomposition of a four-density correlation function into the product of two density fluctuations that make up the transient correlator. The final result takes the form:


The MCT equations of motion are fully specified by the equilibrium structure factor and a single initial time scale. It is often connected to the short time diffusion coefficient. Eq. (5) thus predicts the nonlinear rheology from equilibrium structural correlations, and neglects hydrodynamic interactions; MCT presumes that they only shift the time scale and lead to a high-frequency viscosity . The transient density correlator in Eq. (5) depends only on the magnitude of the advected wavevector, rendering the MCT-ITT result dependent exclusively on the finger tensor and hence resulting in material objectivity [72]. Its validity in the exact Eq. (3) was shown by a discussion of the transient correlation function [72].

2.3 Schematic MCT-ITT equations

Structural correlations at the peak of the structure factor dominate the microscopic MCT-ITT equations. These are connected to the average particle separation via the inverse of the wave vector The spatially resolved equations for can be simplified to a single equation for a local density correlator , because the dynamics on all length scales are coupled strongly (factorization theorem of MCT [13]). In the resulting schematic model [75, 76], the modulus , which determines the shear stress via causal temporal integration


is expressed by the transient correlator as


The transient (shear) stress-stress correlation function depends on the full flow history arising from the time-dependent shear rate . In general, it depends on the two times corresponding to the underlying fluctuations separately. The strain-dependent function is an elastic coefficient that captures the coupling of stress to density fluctuations and was set to a constant in the original formulation of the model [75]. In the microscopic theory, Eq. (5), is given by an integration over wave vectors, including nontrivial weights that depend on time, since density fluctuations are advected by shear. To account for the dephasing of wavevector contributions, which render negative in the microscopic theory, a more refined model [76] allows to express the time-dependence of the prefactor through the accumulated strain :


The direction of the strain does not matter due to symmetry, leading to an even function . The two material parameters and denote characteristic strain values. If the accumulated strain exceeds , stress fluctuations become anticorrelated and the transient correlator in Eq. (7) becomes negative, corresponding to the stress maximum. Beyond stress fluctuations have decorrelated and vanishes, i.e. the steady state stress value has been reached. Eq. (8) was justified by a comparison with fully microscopic MCT-ITT calculations in two dimensions [76]. The single-mode density correlator (normalized to ) is obtained from an equation of motion simplified from the full MCT equations by again neglecting wave vector dependences.

3 Start-up of shear: switch on

3.1 Rheology

During a start-up shear experiment of fixed shear rate the evolution of the shear stress can be followed with the progression of time, or equivalently strain . Brownian dynamics simulations (Fig. 1(a,b) [76]) and rheological experiments (Fig. 1(c) [77]) display three distinct regimes: the linear, the transient and the steady state. In the linear regime, the stress grows almost linearly with increasing strain in an elastic-like fashion before, during the transient regime, reaching a maximum at an intermediate strain . Finally, there is a stress decay to an asymptotic value in the steady state, , where the system exhibits a fluidlike response.

Figure 1: (a,b) Stress vs. strain curves from Brownian dynamic simulations in two dimensions for various values of , as labelled, and packing fractions (a) and (b) [76]. (c) Stress vs. strain curves from PMMA particles of  nm and . (d) Switch-on MSD in the vorticity direction for waiting time  s for PMMA particles with  nm, and as indicated [77]. In all panels, solid black lines are fits from the model.

The stress overshoots observed in experiments of hard-sphere-like poly-methylmethacrylate (PMMA) particles, molecular dynamic (MD) simulations and Brownian dynamic (BD) simulations, were compared by Laurati et al. [77]. It was noted that both the relative magnitude of the overshoot and the strain at which it occurs increase with increasing . The diminishing ability of Brownian motion to fully relax the distorted structures as rises is thought to be the reason for this change. Therefore, systems exposed to higher shear rates are able to store more stress before yielding. Upon increasing the volume fraction, the stress overshoot becomes less prominent as the free volume around a particle decreases [78]. Upon approaching random close packing, where the system is fully jammed, the overshoot essentially disappears due to the inability of particle cages to deform and hence store stress.

An alternative experimental system of colloids with a polystyrene core and a temperature sensitive shell of crosslinked poly(N-isopropylacrylamide) (PNIPAM), which shrinks upon increasing temperature, have also been studied by rheology [79, 80, 60, 76]. A similar increase in the magnitude of the stress overshoot with increasing was observed [76], in agreement with the behaviour of PMMA systems [77, 78]. Upon increasing however, the softer PNIPAM particles display an overshoot up to very high volume fractions, unlike the PMMA system. The origin of this difference is unclear but could be explained by particle softness; the PNIPAM shells have the ability to compress and hence the system can still store stress through an overall shell deformation [78]. Experimental data from both the PMMA and PNIPAM systems as well as simulation data can be well described by MCT [77, 76] (lines in Fig. 1).

3.2 Microscopic Structure

In order to further understand the mechanisms underlying stress storage and yielding, it is necessary to study the microscopic properties on the particle level. Whilst this has traditionally been possible through computer simulations, such as molecular dynamics (MD) [81, 82] or Brownian dynamics (BD) [83, 76], it required recent developments in confocal microscopy to provide experimental access to such quantities [84, 85, 86]. Through the use of either a purpose built shear cell [87, 88, 27, 89, 90], or the coupling of a rheometer with a microscope [91, 92, 93, 94], it is now possible to observe colloidal dispersions under shear. Use of tracking algorithms [95, 96] allow the locations and trajectories of particles to be determined and hence give access to local dynamical and structural information.

The relation between the stress-overshoot and yielding can be clarified by considering the local structure under shear. Experimentally, this requires the measurement of 3-dimensional image stacks, which implies moving either the sample or the objective in order to observe adjacent focal planes. For this reason, the time resolution is lower compared to the observation of a single 2D plane; determining the structure of suspensions under shear hence requires either very low shear rates [58] or very high scanning rates [27, 89]. Using an acousto-optic deflector to control the scanning of the laser beam of the confocal microscope [97], it is possible to quickly image colloidal suspensions under shear and directly observe the structural deformation which takes place in the transient regime.

Figure 2: 2-dimensional in the shear-gradient plane for PMMA particles of  nm and obtained in (a) the quiescent, unsheared state, (b) the steady state of shear with =0.0074 and (c) their difference, i.e. the unsheared case subtracted from the sheared case.

Fig. 2 shows the projection of the pair-distribution function, , in the quiescent, unsheared state and in the steady state of shear in the shear-gradient plane, along with the difference of the two [98]. These results are similar to those observed by BD and MD simulations, where the structure was also directly compared to the stress evolution under shear [99, 83]. The results show that at very low strains (), the structure displays only a slight distortion but this grows with increasing strain, showing a shift of the first peak along the extension axis, corresponding to a reduced probability of finding particles at contact in this direction, but an increased probability along the compression axis due to the squeezing of the particle cage. At strains corresponding to the peak stress, the maximum structural anisotropy is observed before displaying more diffuse lobes along the extension axis with further increase of the strain [99, 83]. In the steady state, the constant distortion of the nearest neighbour cage is caused both by the interplay of increased particle escape along the extension axis and crowding along the compression axis, and the continuous balance between cage destruction and reformation.

The value of the first maximum in gives a measure of the probability of particles contacting and displays deviating behaviour along the compression and extension axes as strain increases. While the probability of contact increases on the compression axis due to the squeezing of the particle cage, the stretching of the cage in the perpendicular direction results in a lower probability of contact and hence a smaller value of . The deviation is most prominent near the stress maximum. The change in the value of can be closely linked to the evolution of stresses in the system [65], despite these being transient states [83]. This indicates a direct correlation of the structural anisotropy with the transient macroscopic stress during yielding [83, 98].

3.3 Microscopic Dynamics

The transient dynamics of colloidal dispersions under shear, as measured by confocal microscopy and in MD simulations, were first described by Zausch et al. [82]. Experimentally, observation of a single plane of the sample through application of fast scanning confocal microscopy, allows particles to be followed with good time resolution and the MSDs to be calculated. Further development of this technique allowed a wider range of shear rates to be studied [77].

Fig. 1(d) shows the MSDs in the vorticity direction for varying shear rates as measured in the quiescent state, the steady-state and the transient state starting from the start-up of shear (waiting time  s, symbols) along with fits to the data from MCT (lines) [77]. The presence of a plateau in the quiescent MSD followed by an increase at long times signifies the caging of particles by their neighbours before diffusion at longer times takes place due to collective particle motion. In the steady state of shear, the MSDs exhibit faster dynamics than at rest, as seen by the earlier onset of diffusion and the correspondingly larger value of the long-time diffusion coefficient. The dynamics become faster upon increasing and long-time diffusion begins earlier, suggesting that the shear rate determines the charateristic timescale of the dynamics. For the transient case with  s, the MSDs follow the quiescent curve for low strains, whilst at long times displaying diffusive behaviour equivalent to the steady state curve, evidencing the shear-thinning nature of these dispersions. At intermediate times however, there is a rapid increase in the MSD with a growth exponent greater than the diffusive exponent of 1, i.e. super-diffusive motion is observed. The prominence of the super-diffusion is seen to increase with , which is similar to the increase of the stress overshoot and is in agreement with predictions from MCT. The range of strains over which super-diffusion is displayed is essentially independent of , being typically 1.5%–6%. That this range seems to be independent of shear rate and waiting time suggests the extent of the transient regime is related to a static lengthscale of the system such as average neighbour separation [77].

MCT relates the stress overshoot to the super diffusive regime of the MSD via the generalized shear modulus; in the quiescent state, the shear modulus is always positive, whilst under shear it exhibits a small negative portion which gives rise to both the stress overshoot and the super-diffusive regime of the MSD. The dip in the generalized shear modulus arises from the secondary -relaxation, which is attributed to the breaking up of nearest neighbour cages [82].

3.4 Summary of transient properties under switch-on of shear

A fuller picture of yielding can now be conjectured: at the start-up of shear, stresses build in the system due to the stretching of particle cages. The structural anisotropy reaches its maximum at strains similar to the maximum in the stress and corresponding to the point at which particle cages start to break up, essentially the yielding point. During cage yielding, super-diffusive dynamics are observed as particles move from their cages due to the shear. In the steady state, the system is essentially fluidlike, displaying diffusive dynamics and a less pronounced though still significant structural anisotropy induced by shear and due to the constant breaking and reformation of particle cages. In the limit of low shear rates, no stress overshoot is observed since stress is completely relaxed through Brownian diffusion before the cage is considerably distorted. At higher shear rates however, cage deformation is rapid enough that stress can be stored before yielding and relaxes partly at the steady state.

4 Cessation of shear: switch off

4.1 Rheology

Upon sudden cessation of shear, the relaxation of dense colloidal suspensions displays behaviour that is heavily dependent on both the properties of the system, such as , and on the previously sheared state [99, 100]. Through a combination of rheology on hard and soft spheres (systems as in §3.1), MD simulations and MCT, it was recently shown that previously sheared glasses can display history dependent residual stresses [100].

For dense liquids, the shear stress decays to zero on the structural-relaxation timescale, , which is typically much larger than the timescale for free diffusion, (Fig. 3(c)). However, as increases towards , grows beyond experimentally accessible times and a plateau in the decay of the shear stress develops. Experimentally measured stress plateaus decay very slowly, possibly as a result of creep, whereas within MCT, pre-sheared glassy states display a nonrelaxing persistent residual stress. This is due to the divergence of the structural relaxation time as full relaxation of density fluctuations is impeded by permanent local caging of particles [100].

Whilst the yield stress of a concentrated colloidal system typically displays an increase with , the residual stresses show a decrease, which suggests that the stress relaxation is more effective for strongly fluidised samples. Scaling the stresses by their flow-curve values, , and the time with , a separation into two distinct classes of liquids and glasses is observed, Fig. 3(c). The ability to rescale by is particularly remarkable given that the shear has already been stopped and the equations of motion contain no reference to the previous perturbation, i.e. the relaxing state has some memory of the shear. The final stress decay for the observed liquids does not scale with as this is governed by the equilibrium relaxation time [100].

Figure 3: (a) MSDs in the vorticity direction after cessation of shear from confocal microscopy on PMMA colloids at . The equilibrium (eq), sheared steady state (NESS) for (red) and (blue) and the transient data for  s (solid lines) and (open symbols) are shown along with a fit to the data from MCT-ITT (green line). (b) MSDs from schematic MCT-ITT model for two state points (glass shown by red lines, liquid by blue lines) with and , for . (Labels as in (a).) (c) Stress relaxation after shear cessation on PMMA colloids at as indicated [100].

4.2 Microscopic Dynamics

Similar to the case of start-up shear, further understanding can be gained by measuring the dynamics of the relaxing state by microscopy and simulations. The MSDs for dense fluid samples after cessation of shear as experimentally determined and calculated from ITT-MCT are shown in Fig. 3(a,b). The short time dynamics follows the steady state curve, which is diffusive, for a time of , i.e. even after cessation of shear the system shows shear-rate dependent dynamics. The transient regime displays sub-diffusion leading to long-time diffusive behaviour, resembling the dynamics of the quiescent state. This intermediate plateau is reminiscent of the quiescent case where it can be linked to particle caging.

The evolution of the transient MSD for short waiting times can be predicted by ITT-MCT based upon the MSDs of the quiescent state and the steady sheared state (green line, Fig. 3(a)). Within MCT, the intermediate plateau represents a second length scale beyond the quiescent localisation length, arising from the competition of shear-induced fluidisation and arrest after cessation, suggesting that this is linked to the residual stresses of the system. It is an effective stopping distance caused by the progression of the distribution function from the perturbed state to the quiescent state.

For glassy states, the evolution of the dynamics after cessation of shear has been studied by MD simulations [100], but has yet to be measured experimentally due to limitations in the stability of the experimental setup. Moreover, highly concentrated samples can exhibit shear banding under low shear rates [93]; the effect of this on the switch off structure can be investigated once sufficiently high shear rates are experimentally accessible. Overcoming these issues and studying both the dynamics and the structural origin of residual stresses after shear cessation will form the focus of future research.

5 Large amplitude oscillatory shear (LAOS)

Key aspects of the nonlinear rheology of viscoelastic dispersions and yielding soft solids is also displayed upon application of an oscillatory shear strain


introducing the parameters strain amplitude and frequency . The response of a dispersion at long times, such that transient effects from switching on the strain have decayed, is dicussed [60].

Figure 4: Stress response measured in LAOS experiments for strain amplitudes (black lines) and associated Lissajous figures illustrating the nonlinear character of the response. The experiments on PNIPAM microgels were performed on a glassy sample at C and a frequency of 1 Hz (corresponding to ). The yield stress is indicated by the green dashed lines, theoretical results by red dashed lines [101].

The LAOS stress signals measured in PNIPAM microgel dispersions at a frequency in the -window (0.01 Hz, corresponding to ) are shown in the right-hand panels of Fig. 4. The deformation varies from in the linear regime with an almost entirely elastic response to deep in the nonlinear regime with an almost purely viscous response. The increase in dissipation with increasing is reflected in the increasing area enclosed by the Lissajous curves. For small strain amplitudes, the linear viscoelastic behaviour is indicated by the nearly perfect sinusoidal stress . The time dependent signal becomes distorted away from a pure sinusoid when the peak of encounters the dynamic yield stress (green dashed lines, Fig.4). Theory (red dashed lines) predicts that for even lower frequencies, the peak maximum should approach the yield stress value, but this could not be tested experimentally. For the intermediate chosen frequency, the stress signal displays flattened asymmetric peaks at intermediate values of , consistent with a regime of cage breaking around . The schematic model calculations included in the figure are in good agreement with the experimental data. Data at higher frequencies show somewhat larger deviations from the model calculations. Similar information can be obtained from the ‘Lissajous curves’ (left panels, Fig. 4).

6 Conclusions

The transient response of dense colloidal suspensions both to the application of an external shear and to its cessation, have been studied with the aim of elucidating the shear-induced structural and dynamical changes in these systems. These insights have been made possible through advancements in mode coupling theory and experimental techniques, namely confocal microscopy and the associated data analysis.

Application of external shear results in cage deformations on the microscopic particle level as well as an increased (reduced) probability to find particles at contact in the compression (extension) direction, respectively, which reflects the build up of stresses. At the same time, the transient dynamics exhibit super-diffusive behaviour. Thus the applied shear enforces the decay of correlations even in a state that would otherwise be dynamically arrested. MCT links the super-diffusion to the overshoot observed in the rheological response. The external shear rate determines the relaxation rate, which is much faster than the internal structural relaxation rate, and allows for local particle motion leading to cage distortion and shear-thinning. With increasing strain and beyond the yield strain, a steady-state is reached that is characterised by a reduced though still significant structural anisotropy, purely diffusive dynamics and a fluid-like rheological response. After cessation of shear, the shear stress does not fully decay but history-dependent residual stresses remain, which are smaller in previously strongly fluidised samples. The dynamics are sub-diffusive with a plateau which represents a second length scale in addition to the quiescent localisation length. The localisation length is caused by the competition of shear-induced fluidisation and cessation-induced arrest, suggesting a link to the residual stress.

We acknowledge support by the Deutsche Forschungsgemeinschaft (DFG) through the SFB-Transregio TR6 “The Physics of Colloidal Dispersions in External Fields” (project A6) during the period 2005–2013. We thank J.M. Brader, M. Krüger, O. Henrich, Th. Voigtmann and F. Weysser for joint work on this project. We also thank M.E. Cates, N. Koumakis, G. Petekidis, M. Siebenbürger, and M. Ballauff for the fruitful cooperation and discussions, and J. Horbach (project A5) for collaboration.


  • [1] P. N. Pusey, Liquids, Freezing and the Glass Transition, In: Proceedings of the Les Houches Summer School, J. P. Hansen, D. Levesque, and J. Zinn-Justin (eds.), Elsevier, Amsterdam, 1991.
  • [2] W. C. K. Poon, A day in the life of a hard-shpere suspension, In: Soft and Fragile Matter: Nonequilibrium Dynamics, Metastability and Flow, M. E. Cates and M. R. Evans (eds.), Institute of Physics Pub., copublished by Scottish Universities Summer School in Physics, Oxford, 2000.
  • [3] T. McLeish, Rheology of linear and branched polymers, In: Soft and Fragile Matter: Nonequilibrium Dynamics, Metastability and Flow, M. E. Cates and M. R. Evans (eds.), Institute of Physics Pub., copublished by Scottish Universities Summer School in Physics, Oxford, 2000.
  • [4] J. M. Brader, J. Phys.: Condens. Matter 22, 363101 (2010).
  • [5] W. G. Hoover and F. H. Ree, J. Chem. Phys. 49, 3609 (1968).
  • [6] P. N. Pusey and W. van Megen, Nature 320, 340 (1986).
  • [7] U. Gasser, J. Phys.: Condens. Matter 21, 203101 (2009).
  • [8] P. N. Pusey, E. Zaccarelli, C. Valeriani, E. Sanz, W. C. K. Poon, and M. E. Cates, Phil. Trans. R. Soc. A 367, 4993 (2009).
  • [9] L. V. Woodcock, Ann. N. Y. Acad. Sci. 37, 274 (1981).
  • [10] P. N. Pusey and W. van Megen, Phys. Rev. Lett. 59, 2083 (1987).
  • [11] W. van Megen and P. N. Pusey, Phys. Rev. A 43, 5429 (1991).
  • [12] W. Götze and L. Sjögren, Rep. Prog. Phys. 55, 241 (1992).
  • [13] W. Götze, Complex Dynamics of Glass-Forming Liquids, A Mode-Coupling Theory, Oxford University Press, 2009.
  • [14] W. van Megen and S. M. Underwood, Phys. Rev. Lett. 70, 2766 (1993).
  • [15] W. van Megen and S. M. Underwood, Phys. Rev. E 49, 4206 (1994).
  • [16] I. M. Krieger, Adv. Colloid Interface Sci. 3, 111 (1972).
  • [17] I. M. Krieger and M. Eguiluz, Trans. Soc. Rheol. 20, 29 (1976).
  • [18] B. J. Ackerson, Physica A 174, 15 (1991).
  • [19] K. N. Pham, G. Petekidis, D. Vlassopoulos, S. U. Egelhaaf, P. N. Pusey, and W. C. K. Poon, Europhys. Lett. 75, 624 (2006).
  • [20] K. N. Pham, G. Petekidis, D. Vlassopoulos, S. U. Egelhaaf, W. C. K. Poon, and P. N. Pusey, J. Rheol. 52, 649 (2008).
  • [21] T. Sentjabrskaja, E. Babaliari, J. Hendricks, M. Laurati, G. Petekidis, and S. U. Egelhaaf, Soft Matter 9, 4524 (2013).
  • [22] G. K. Batchelor, J. Fluid Mech. 83, 97 (1977).
  • [23] G. Petekidis, P. N. Pusey, A. Moussaïd, S. U. Egelhaaf, and W. C. K. Poon, Phys. A 306, 334 (2002).
  • [24] G. Petekidis, A. Moussaïd, and P. N. Pusey, Phys. Rev. E 66, 051402 (2002).
  • [25] G. Petekidis, D. Vlassopoulos, and P. N. Pusey, Faraday Discuss. 123, 287 (2003).
  • [26] K. N. Pham, S. U. Egelhaaf, A. Moussaïd, and P. N. Pusey, Rev. Sci. Instrum. 75, 2419 (2004).
  • [27] R. Besseling, E. R. Weeks, A. B. Schofield, and W. C. K. Poon, Phys. Rev. Lett. 99, 028301 (2007).
  • [28] K. J. Mutch, N. Koumakis, M. Laurati, and S. U. Egelhaaf, to be published .
  • [29] P. D’Haene, J. Mewis, and G. G. Fuller, J. Colloid Interface Sci. 156, 350 (1993).
  • [30] J. W. Bender and N. J. Wagner, J. Rheol. 40, 899 (1996).
  • [31] B. J. Maranzano and N. J. Wagner, J. Chem. Phys. 117, 10291 (2002).
  • [32] D. R. Foss and J. F. Brady, J. Fluid Mech. 407, 167 (2000).
  • [33] J. F. Brady, Chem. Eng. Sci. 56, 2921 (2001).
  • [34] Q. D. Nguyen and D. V. Boger, J. Rheol. 27, 321 (1983).
  • [35] J. Persello, J. Chang, B. Cabane, A. Magnin, and J. M. Piau, J. Rheol. 38, 1845 (1994).
  • [36] C. Derec, G. Ducouret, A. Adjari, and F. Lequeux, Phys. Rev. E 67, 061403 (2003).
  • [37] F. Mahaut, X. Chateau, P. Coussot, and G. Ovarlez, J. Rheol. 52, 287 (2008).
  • [38] S. A. Khan, C. A. Schnepper, and R. C. Armstrong, J. Rheol. 32, 69 (1988).
  • [39] M. T. Islam, N. Rodriguez-Hornedo, S. Ciotti, and C. Ackermann, Pharm. Res. 21, 1192 (2004).
  • [40] A. Mohraz and M. J. Solomon, J. Rheol. 49, 657 (2005).
  • [41] L. Bécu, S. Manneville, and A. Colin, Phys. Rev. Lett. 96, 138302 (2006).
  • [42] P. Akcora, H. Liu, S. K. Kumar, J. Moll, Y. Li, B. C. Benicewicz, L. S. Schadler, D. Acehan, A. Z. Panagiotopoulos, V. Pryamitsyn, V. Ganesan, J. Ilavsky, P. Thiyagarajan, R. H. Colby, and J. F. Douglas, Nature Mater. 8, 354 (2009).
  • [43] V. Carrier and G. Petekidis, J. Rheol. 53, 245 (2009).
  • [44] T. Divoux, C. Barentin, and S. Manneville, Soft Matter 7, 9335 (2011).
  • [45] J.-F. Berret, Langmuir 13, 2227 (1997).
  • [46] S. Lerouge, J.-P. Decruppe, and C. Humbert, Phys. Rev. Lett. 81, 5457 (1998).
  • [47] J. F. A. Soltero, F. Bautista, J. E. Puig, and O. Manero, Langmuir 15, 1604 (1999).
  • [48] K. Osaki, T. Inoue, and T. Isomura, J. Polym. Sci. B 38, 2043 (200).
  • [49] M. T. Islam and L. A. Archer, J. Polym. Sci. B 39, 2275 (2001).
  • [50] J.-P. Decruppe, S. Lerouge, and J.-F. Berret, Phys. Rev. E 63, 022501 (2001).
  • [51] P. E. Boukany, O. Hemminger, S.-Q. Wang, and L. J. Lee, Phys. Rev. Lett. 105, 027802 (2010).
  • [52] R. Ganapathy and A. K. Sood, J. Non-Newtonian Fluid Mech. 149, 78 (2008).
  • [53] S. Ravindranath and S.-Q. Wang, J. Rheol. 52, 681 (2008).
  • [54] P. Tapadia and S.-Q. Wang, Macromolecules 37, 9083 (2004).
  • [55] P. Tapadia and S.-Q. Wang, Phys. Rev. Lett. 96, 016001 (2006).
  • [56] R. G. Larson, The structure and rheology of complex fluids, Oxford University Press, New York, 1999.
  • [57] E. Bouchbinder and J. S. Langer, Phys. Rev. Lett. 106, 148301 (2011).
  • [58] P. Schall, D. A. Weitz, and F. Spaepen, Science 318, 1895 (2007).
  • [59] P. Sollich, F. Lequeux, P. Hébraud, and M. E. Cates, Phys. Rev. Lett. 78, 2020 (1997).
  • [60] M. Siebenbürger, M. Fuchs, and M. Ballauff, Soft Matter 8, 4025 (2012).
  • [61] A. Ikeda, L. Berthier, and P. Sollich, Phys. Rev. Lett. 109, 018301 (2012).
  • [62] J. K. G. Dhont, An Introduction to Dynamics of Colloids, Elsevier Science, Amsterdam, 1996.
  • [63] J. Bergenholtz, J. Brady, and M. Vicic, J. Fluid. Mech. 456, 239 (2002).
  • [64] W. Russel, D. A. Saville, and W. R. Schowalter, Colloidal dispersions, Cambridge University Press, 1989.
  • [65] J. F. Brady, J. Chem. Phys. 99, 567 (1993).
  • [66] J. M. Brader and M. Krüger, Mol. Phys. 109, 1029 (2011).
  • [67] P. Hopkins, A. Fortini, A. J. Archer, and M. Schmidt, J. Chem. Phys. 133, 224505 (2010).
  • [68] W. van Megen, T. C. Mortensen, S. Williams, and J. Müller, Phys. Rev. E 58, 6073 (1998).
  • [69] K. Kawasaki and J. D. Gunton, Phys. Rev. A 8, 2048 (1973).
  • [70] D. J. Evans and G. P. Morriss, Statistical Mechanics of Nonequilibrium Liquids, Cambridge University Press, 2008.
  • [71] J. M. Brader, M. Cates, and M. Fuchs, Phys. Rev. Lett. 101, 138301 (2008).
  • [72] J. M. Brader, M. E. Cates, and M. Fuchs, Phys. Rev. E 86, 021403 (2012).
  • [73] M. Fuchs and M. E. Cates, Phys. Rev. Lett. 89, 248304 (2002).
  • [74] M. Fuchs and M. E. Cates, J. Rheol. 53, 957 (2009).
  • [75] J. M. Brader, Th. Voigtmann, M. Fuchs, R. G. Larson, and M. E. Cates, Proc. Natl. Acad. Sci. USA 106, 15186 (2009).
  • [76] C. P. Amann, M. Siebenbürger, M. Krüger, F. Weysser, M. Ballauff, and M. Fuchs, J. Rheol. 57, 149 (2013).
  • [77] M. Laurati, K. J. Mutch, N. Koumakis, J. Zausch, C. P. Amann, A. B. Schofield, G. Petekidis, J. F. Brady, J. Horbach, M. Fuchs, and S. U. Egelhaaf, J. Phys.: Condens. Matter 24, 464104 (2012).
  • [78] N. Koumakis, A. Pamvouxoglou, A. S. Poulos, and G. Petekidis, Soft Matter 8, 4271 (2012).
  • [79] J. J. Crassous, M. Siebenbürger, M. Ballauff, M. Drechsler, D. Hajnal, O. Henrich, and M. Fuchs, J. Chem. Phys. 128, 204902 (2008).
  • [80] M. Siebenbürger, M. Fuchs, H. Winter, and M. Ballauff, J. Rheol. 53, 707 (2009).
  • [81] F. Varnik, L. Bocquet, and J.-L. Barrat, J. Chem. Phys. 120, 2788 (2004).
  • [82] J. Zausch, J. Horbach, M. Laurati, S. U. Egelhaaf, J. M. Brader, T. Voigtmann, and M. Fuchs, J. Phys.: Condens. Matter 20, 404210 (2008).
  • [83] N. Koumakis, M. Laurati, S. U. Egelhaaf, J. F. Brady, and G. Petekidis, Phys. Rev. Lett. 108, 098303 (2012).
  • [84] A. van Blaaderen and P. Wiltzius, Science 270, 1177 (1995).
  • [85] E. R. Weeks, J. C. Crocker, A. C. Levitt, A. Schofield, and D. A. Weitz, Science 287, 627 (1999).
  • [86] A. D. Dinsmore, E. R. Weeks, V. Prasad, A. C. Levitt, and D. A. Weitz, Appl. Opt. 40, 4152 (2001).
  • [87] D. Derks, H. Wisman, A. van Blaaderen, and A. Imhof, J. Phys.: Condens. Matter 16, S3917 (2004).
  • [88] P. A. Smith, G. Petekidis, S. U. Egelhaaf, and W. C. K. Poon, Phys. Rev. E 76, 041402 (2007).
  • [89] D. Chen, D. Semwogerere, J. Sato, V. Breedveld, and E. R. Weeks, Phys. Rev. E 81, 011403 (2010).
  • [90] X. Cheng, J. H. McCoy, J. N. Israelachvili, and I. Cohen, Science 333, 1276 (2011).
  • [91] P. Ballesta, R. Besseling, L. Isa, G. Petekidis, and W. C. K. Poon, Phys. Rev. Lett. 101, 258301 (2008).
  • [92] R. Besseling, L. Isa, E. R. Weeks, and W. C. K. Poon, Adv. Colloid Interface Sci. 146, 1 (2009).
  • [93] R. Besseling, L. Isa, P. Ballesta, G. Petekidis, M. E. Cates, and W. C. K. Poon, Phys. Rev. Lett. 105, 268301 (2010).
  • [94] S. K. Dutta, A. Mbi, R. C. Arevalo, and D. L. Blair, Rev. Sci. Instrum. 84, 063702 (2013).
  • [95] J. C. Crocker and D. G. Grier, J. Colloid Interface Sci. 179, 298 (1996).
  • [96] M. C. Jenkins and S. U. Egelhaaf, Adv. Colloid Interface Sci. 136, 65 (2008).
  • [97] D. Semwogerere and E. R. Weeks, Confocal microscopy, In: Encyclopedia of Biomaterials and Biomedical Engineering, G. Wnek and G. Bowlin (eds.), Taylor & Francis, 2005.
  • [98] N. Koumakis, M. Laurati, K. J. Mutch, S. U. Egelhaaf, J. F. Brady, and G. Petekidis, to be published .
  • [99] J. Zausch and J. Horbach, Europhys. Lett. 88, 60001 (2009).
  • [100] M. Ballauff, J. M. Brader, S. U. Egelhaaf, M. Fuchs, J. Horbach, N. Koumakis, M. Krüger, M. Laurati, K. J. Mutch, G. Petekidis, M. Siebenbürger, T. Voigtmann, and J. Zausch, Phys. Rev. Lett. 110, 215701 (2013).
  • [101] J. M. Brader, M. Siebenbürger, M. Ballauff, K. Reinheimer, M. Wilhelm, S. J. Frey, F. Weysser, and M. Fuchs, Phys. Rev. E 82, 061401 (2010).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description