Deformation profiles and microscopic dynamics of complex fluids during oscillatory shear experiments

Oscillatory shear tests are widely used in rheology to characterize the linear and non-linear mechanical response of complex fluids, including the yielding transition. There is an increasing urge to acquire detailed knowledge of the deformation field that is effectively present across the sample during these tests; at the same time, there is mounting evidence that the macroscopic rheological response depends on the elusive microscopic behavior of the material constituents. Here we employ a strain-controlled shear-cell with transparent walls to visualize and quantify the dynamics of tracers embedded in various cyclically sheared complex fluids, ranging from almost-ideal elastic to yield stress fluids. For each sample, we use image correlation processing to measure the macroscopic deformation field, and echo-Differential Dynamic Microscopy to probe the microscopic irreversible sample dynamics in reciprocal space; finally, we devise a simple scheme to spatially map the rearrangements in the sheared sample, once again without tracking the tracers. For the yield stress sample, we obtain a wave-vector dependent characterization of shear-induced diffusion across the yielding transition, which is accompanied by a three-order-of-magnitude speed-up of the dynamics and by a transition from localized, intermittent rearrangements to a more spatially homogeneous and temporally uniform activity. Our tracking free approach is intrinsically multi-scale, can successfully discriminate between different types of dynamics, and can be automated to minimize user intervention. Applications are many, as it can be translated to other imaging modes, including fluorescence, and can be used with sub-resolution tracers and even without tracers, for samples that provide intrinsic optical contrast.


Introduction
As compared with molecular hard solids -such as steel, wood or bones -soft materials made of micron-sized constituents promptly deform with relatively small mechanical loads: we easily spread butter on a slice of bread, mindlessly distribute toothpaste on our teeth, effortlessly rub our faces with cream, smoothly tickle our palate with a chocolate mousse, and readily pull our cheeks with a little pinch. Beyond having practically ubiquitous consequences in our lives, this softness opens up unique possibilities when doing experimental research on soft matter. A paradigmatic example is passive microrheology, in which the linear moduli of a soft material are probed by embedded microscopic tracer particles that operate as micro-scale rheometers powered by tiny thermal fluctuations 1,2 . Despite being a technique that has had enormous success, passive microrheology is unable to provide complete rheological information in all cases of interest. For instance, if the storage modulus of the material is too large or if one is interested in the non-linear rheological properties, it is often preferred to recur to active microrheology, in which at least one of the embedded particles is driven by an external force obtained, for instance, via optical or magnetic trapping 1,3-5 . Passive and active microrheology techniques are very powerful tools for the characterization of soft materials, as they can access high frequencies and enable a space-resolved study of the rheological properties of complex, heterogeneous systems 3 . However, the fact that the energy is injected at a local, particle scale makes them unsuitable to study More insightful information on the microscopic mechanisms underlying plasticity in soft materials can be obtained by imposing a controlled macroscopic deformation to the sample -for instance by displacing one or more of the confining surfaces of a shear rheometer -while monitoring over time the structural rearrangements occurring at the microscale. Rearrangements are typically monitored by accurate tracking of embedded tracer particles, performed in direct (rheo-microscopy or rheo-imaging) or in reciprocal space (rheo-scattering), very often but not necessarily by optical means [6][7][8][9] . These rearrangements may be partly or entirely reversible, with the irreversible component becoming dominant for increasing shear strain across the transition to nonlinearity [10][11][12][13][14] . This kind of experiments provides thus very useful insights on how the microscopic structure and dynamics determine the material rheology, which can be used to design enhanced materials, e.g. with a longer lifetime, controlled relaxation, improved ductility or emerging memory 15 .
A variety of shear rheology tests is compatible with the simultaneous characterization of the tracer particles' motion. The latter can be generally split into the sum of affine and non-affine displacements 16 ; the relative importance of these two contributions depends on the type of test, and on the scheme used to sample the particles' configuration in space and time. For instance, during Small/Medium/Large Amplitude Oscillatory Shear (SAOS/MAOS/LAOS) rheology experiments, the so-called echo protocol enables one sampling the microscopic configuration at each deformation cycle, by obtaining a stroboscopic movie of the particles configurations either in direct or in reciprocal space. In this way, the affine contribution is suppressed and one can focus on the non-affine irreversible particle rearrangements.
Most of the initial experiments with the echo protocol were performed with Diffusing Wave Spectroscopy (DWS) 10 , and Dynamic Light Scattering (DLS) 17 . Being based on the collection of multiply scattered light, DWS is very sensitive to a single, small length scale, whereas DLS is intrinsically a multi-scale technique, at the expense of an increased experimental complexity 18 . More recently, microscopy has also been employed to track particle motion in oscillatory shear experiments with the echo protocol 12 .
If one is interested in going beyond the study of irreversible non-affine particle displacements typically obtained with the echo protocol, imaging and microscopy offer an important advantage over scattering experiments: the affine particle motion that characterizes the intra-cycle particle dynamics can be readily removed from the particle trajectories, whereas more caution is needed in scattering experiments 19 . However, the information extracted from scattering experiments is characterized by an enviable statistical robustness of the obtained information. For these reasons, it would be interesting to have both techniques available during oscillatory shear experiments.
A possible compromise is represented by Digital Fourier Microscopy (DFM) 20 , a family of techniques based on collecting in direct space temporal sequences of images of the sample, which are subsequently spatially Fourier transformed to obtain scattering information in the reciprocal space as a function of the wavevector q. The most popular DFM technique is Differential Dynamic Microscopy (DDM) [21][22][23][24] , which relies on subtracting images acquired at different times, a feature that appears well suited to echo experiments: when two images acquired at an integer multiple of the shear oscillation period are subtracted, only the irreversible particle displacements are captured and contribute to the final result. At the same time, one can take advantage of the availability of the direct space movies to remove the rigid translations associated to the affine intra-cycle particle dynamics 25 , a step that may benefit from the large amount of image registration algorithms available 26 .
The potential of DDM in combination with shear rheology experiments was recently suggested in two companion articles by Aime and Cipelletti 6,19 , in which the authors performed a critical comparison with DLS for probing shear-induced rearrangements in Fourier space. In both cases, the intermediate scattering functions for q parallel to the direction of the applied shear (shear direction) are mostly dominated by the relaxation due to the affine particle displacements, whereas in the perpendicular direction (vorticity direction) only the non-affine particle rearrangements are probed. The main advantage of DDM over DLS is that the limited spatial coherence of illumination typical of brightfield DDM experiments results in a finite, q-dependent depth-offocus, which can be used to probe slices of the sample perpendicular to the optical axis 22 . This capability is surely appealing for those cases where the deformation field deviates from a simple homogeneous profile, such as in the presence of wall slip, shear banding, or edge fracture 27 . All these different phenomena have the common effect of making the actual deformation field inhomogeneous and difficult to predict, causing problems in the interpretation of the results of both rheological and DLS experiments. Indeed, while in a rheological experiment, the characterization of the material constitutive equation for the shear stress σ (γ), requires the accurate knowledge of the actual shear strain γ of the sample, in DLS experiments, where the signal is the result of an average over the entire scattering volume, it is almost impossible to discern the contributions coming from different portions of a macroscopically inhomogeneous sample. Aime and Cipelletti have clearly and convincingly showed that a proper DDM characterization of the affine displacements is possible, providing that an instrument calibration of the q-dependent depth-of-focus is performed 19 . However, the full potential of combining DDM analysis with an echo protocol, where the effect of affine displacements is minimized, remains yet unexplored.
In this work, we perform oscillatory shear experiments with a strain-controlled cell, characterized by a parallel-sliding-plates geometry. During the sample oscillation, we monitor the dynamics of tracer particles by using bright-field DDM in conjunction with the echo protocol (echo-DDM). The depth-of-focus is tuned and optimized by varying the numerical aperture of the microscope condenser so that, even if the deformation field is not homogeneous across the whole gap, it is possible to selectively image narrow "slices" of the sample, where the shear can be assumed to be approximately constant (see Fig.1). We apply this optimization procedure to the imaging of various samples, some of which presenting heterogeneity of the deformation field in the direction of the velocity gradient (gradient direction). We perform a preliminary accurate characterization of the deformation field across the whole gap in standard samples with pure elastic solid or pure viscous liquid response, retrieving the expected linear deformation profile. By applying the same protocol to more complex soft matter samples (a viscoelastic fluid and a yield stress material), we demonstrate that our approach can also effectively probe deviations from such ideal profile, due for example to wall slip or shear banding. Quantification of these effects is key to correctly relate the microscopic irreversible dynamics to the local macroscopic deformation in echo experiments.
Once obtained a reliable characterization of the local strain for all our samples, we use echo-DDM to probe the presence of a shear-induced non-affine rearrangements. The elastic solid, characterized by the absence of a true irreversible non-affine dynamics, serves as a reference to estimate the instrumental limits; as expected, Newtonian and Maxwell liquids do not exhibit any appreciable variation in non-affine activity as a function of the shear amplitude; by contrast, in the yield stress sample, we find intermittent, small scale displacements of the tracers for small strain amplitudes, and shear-induced diffusion 28-31 for amplitudes larger than the one corresponding to the crossover point of the viscoelastic moduli. All the results obtained with echo-DDM are compared and found in agreement with particle tracking (PT) analysis of the tracer motion. Finally, we propose a simple yet robust scheme to draw activity maps highlighting the rearrangements that occur in the sheared sample.
The manuscript is organized as follows: in the Materials and Methods section we describe the tested samples (Sect 2.1), the experimental set-up (Sect.2.2), and the experimental protocols enabling the characterization of the macroscopic deformation field (Sect.2.3) and of the irreversible non-affine dynamics (Sect.2.4). The results obtained in the characterization of the deformation profile and of the shear-induced rearrangements are presented and discussed in Sect.3.1 and Sect.3.2, respectively, before pre-senting some concluding comments in Sect.4.

Samples
We test our experimental approach on samples that are representative of different classes of materials: an elastomer (Sylgard), a Newtonian fluid (silicon Oil), a Maxwell fluid (worm-like micelles, WLM), and a simple yield stress fluid (Carbopol). All the samples in our study are fairly transparent and are thus seeded with tracer particles of known size. However, seeding is not needed if the sample already provides a sufficient optical contrast, because DDM can capture the dynamics of a variety of samples including molecular mixtures 32,33 and liquid crystals 34 .
Sylgard 184 (Dow Corning) is as a two-component system made of a silicon oil and a curing agent. A small fraction (φ ∼ 0.1%) of sterically stabilized PMMA particles of diameter 260 nm are dispersed in the silicon oil base before adding the curing agent. The sample is then mixed thoroughly, centrifuged to get rid of bubbles, and finally loaded in the shear-cell before the completion of the curing process, which takes about 48 h at room temperature. The so-obtained sample is expected to behave like a perfectly homogeneous elastic material over a wide range of shear strains 35 . The silicon oil that we use is the viscous component described above in the preparation of Sylgard. This sample behaves as a Newtonian viscous fluid with a shear viscosity η 5 Pa s 36 . Worm-like micelles are obtained by dispersing Hexadecylpyridinium chloride monohydrate (CpyCl) surfactant in water with Sodium Salicylate (NaSal) as a counter-ion 37,38 . Both components have been purchased from Sigma-Aldrich as powders. CpyCl is dispersed in Milliq water at a concentration 200 mM and taken to 40 • C for 30 minutes, until it turns to a clear solution. NaSal is dispersed in Milliq water at a concentration of 120 mM. The two dispersions are slowly mixed in a 1:1 ratio and left on a rotating mixer overnight. The final molarities are 100 mM and 60 mM for CpyCl and NaSal, respectively. Polystyrene particles of 2 µm diameter (PS2, Microparticles GmbH) are dispersed at a volume fraction φ = 0.05%. The dispersion is initially mixed manually with a spatula and left on a rotating mixer for one night. Carbopol 971 P NF has been purchased from Lubrizol as a single component powder. The powder is dispersed in MilliQ water to obtain a concentrated starting solution with concentration c = 5% wt. The obtained dispersion has an acidic pH, that is neutralized by adding few drops of NaOH 10 M to 10 ml of concentrated dispersion, while gently stirring for several days by keeping the pH controlled. At the end of the mixing phase, the neutralized dispersion is then diluted to obtain the sample at the desired concentration c = 0.5% wt. We found that this preparation procedure ensures the reproducibility of the rheological properties of the Carbopol samples. Eventually, PS2 tracer particles are dispersed at a volume concentration φ = 0.05%. Before the experiment, all samples are centrifuged to get rid of the air bubbles. We note that the choice of the tracer particle size for the Carbopol sample was guided by preparatory experiments in which we investigated the tracer dynamics in Carbopol at rest for different particle diameters: 0.5, 0.8, 2, and 5 microns. For the 0.5 and 0.8 micron particles, we found that roughly half of the particles appeared stuck and the other half underwent caged dynamics. By contrast, for the 2 and 5 micron tracers, all the particles were substantially trapped. These observations, in agreement with seminal microrheological studies of Carbopol EDT 2050 39 , point to the existence of a characteristic microstructural length scale, which in our Carpobol sample is smaller than 2 micron.

Shear-cell
The samples are loaded between the two horizontal parallel plates of a commercial strain-controlled shear-cell (RheOptiCAD, CAD Instruments, France). With this shear-cell, a modified version of the one described in Ref. 40 to enable bright-field imaging, we can impose an oscillatory displacement of the top plate while the bottom plate remains stationary. In our experiments, the oscillation frequency is kept constant (f = 1 Hz), while the amplitude A 0 is varied in the range [0.01 − 1] mm. The displacement of the top plate identifies the shear direction, corresponding to the orientation of the x-axis in Fig.1. The height of gap between the plates is 1 mm for Sylgard, 250 µm for Worm-like micelles, and 500 µm for the other samples. The plates are standard glass slides whose surface was sandblasted to reduce wall slip, except for small a circular window of diameter 4 mm, which is left clear to enable optical access to the sample. Even though we have not conducted an independent characterization of the roughness of the sandblasted surfaces, we expect it to be of a few micrometers, slightly larger than both the microstructure length-scale of the Carbopol and the particle size. The orientation of the plates is controlled by micrometer screws. The parallelism of the two plates and the alignment of the moving plate along the shear direction are essential to produce a simple shear deformation. The proper alignment of the top plate is assessed by observing the plate surface under the microscope while it is displaced along the shear direction and checking whether it remains in sharp focus. As described in more detail in ESI, the angle between the two plates is measured by means of a custom interferometric system, coupled to the rear port of the microscope. The alignment accuracy than can be achieved is within δ θ ∼ 10 −3 rad. Due to a memory buffer register limitation, with the current version of the shear-cell control software the length of an oscillatory experiment is limited to 500 cycles. Moreover, the shear-cell does not provide a real-time output signal proportional to the position of the moving plate, which could have been used to trigger a synchronized camera acquisition. Synchronization is a crucial issue in our set-up and it will be discussed in detail in Sec. 2.4.3.

Microscope
The shear-cell described in the previous section is mounted on an inverted microscope (Nikon Eclipse Ti-S), equipped with a 20x 0.45 NA objective. The working distance of the objective (WD= 8 mm) is large enough to enable scanning the whole gap between the two plates confining the sample. As discussed also in the the ESI, the vertical resolution in the reconstructed deformation field is determined by the depth-of-focus of the imaging system. In our setup, we minimize the depth-of-focus (to about 20 µm, see ESI for details) by keeping the condenser diaphragm almost completely open (LWD, NA=0.52).

Deformation field
It is known that in several conditions of interest, the deformation field within a sheared sample may be non-homogeneous even when a uniform stress is applied 41 . Two examples of deviation from homogeneity are shear banding and wall slip. These phenomena affect a wide variety of samples 42 43 44 and while they have an indisputable effect on the mechanical response and thus on the interpretation of rheological data, it is not easy to pinpoint their presence from purely mechanical measurements. In the following, we provide a simplified description of these effects in the case of the parallel sliding plates geometry, which is the one used in our experiments and is schematically represented in Fig.1. A key advantage of this geometry, compared to other common configurations like rotating plate-plate or rotating cone-plate, is that it allows imposing a uniform strain at the confining surfaces while ensuring easy optical access to the sample 6,12,13 .
Ideally, when a uniform displacement ∆X 0 = A 0ûx is imposed on the top plate, while the bottom one is fixed, an affine displacement field is produced in the sample where z is the vertical coordinate (corresponding to the velocity gradient direction), and has the origin of the axis at the bottom plate;û x in the unit vector oriented along the x (shear) direction and the ratio γ 0 = A 0 /h between the imposed displacement at the moving plate and the gap width h is the imposed shear strain. The displacement is uniform over each horizontal plane, with no components along other directions. The corresponding strain, which is given by the spatial gradient of the displacement field ∇ ∆X , is perfectly homogeneous and the only non-zero component of the strain tensor corresponds to the shear strain Deviations from this ideal case occur for example in case of wall slip, when the adhesion between the sample and the confining plates is not sufficient to prevent their relative motion. When wall slip occurs, in the absence of other effects, like shear-banding, Eq.2 modifies to where z 0 is the coordinate of the horizontal plane (outside the sample) on which the extrapolated displacement profile would vanish, and c < 1. In this scenario, the strain field is still homogeneous, even if with the actual shear strain γ is smaller than the imposed one γ = cγ 0 < γ 0 . Since the actual deformation is obtained multiplying the imposed one by c, we define c as the deformationloss-factor. Eq. 4 can reproduce any (symmetric or asymmetric) linear profile with "perfect slip" at one or two surfaces, by suitably choosing the parameters 0 ≤ c ≤ 1 and −h ≤ z 0 ≤ 0.
Another important example of deviation from homogeneity of the deformation field is shear banding, corresponding to the separation of the sample into macroscopic regions (bands) with different local shear strains 44,45 . As an example, one can consider the simplistic case of a system separated in two bands, in the absence of wall slip, which is described by the displacement profile where γ 1 and γ 2 are the local shear strains in the two bands. Even in this simple case where just two bands are formed, many profiles are compatible with a given imposed deformation γ 0 , corresponding to all possible choices of the parameters z 1 and γ 1 satisfying the constraint z 1 γ 1 ≤ hγ 0 . Under these conditions, even if the average shear strain is bound to be γ 0 , the local deformation can be very different, assuming in principle any value from 0 to infinity, which reinforces the importance of mapping the local strain. Eqs. 4 and 5 can be easily combined to obtain a description of the deformation field with coexisting wall slip and shear banding. The simple model described by Eq. 5 involves a sharp discontinuity of the local strain that may be avoided with more complex modelling in which the local strain changes continuously 46,47 . For all the experiments described in this work, Eq. 4 provides a satisfactory description.
In our bright-field microscopy experiments, we exploit the optical sectioning capability resulting from the partial coherence of illumination to obtain an accurate characterization of the displacement ∆X(z) at different heights z, and thus a faithful reconstruction of the actual shear strain profile γ(z) within the sample.

Strain field under oscillatory shear
We impose an oscillatory displacement of the top plate described by Eq.1, with imposed amplitude A 0 and shear strain γ 0 = A 0 /h. We then perform observations of the sheared sample at different vertical positions z i ∈ [0, h]. For each height z i , a stack of 200 images is acquired with a frame rate 1/δt = 100 fps (see ESI).
We calculate the spatial cross-correlation function c(r,t) of each pair of consecutive frames collected at a given height The positionr(t) of the peak of c(r,t) provides an estimate for the rigid displacement occurred between the two images, taken at times t and t + δt, respectively. By summing up all contributions from consecutive image pairs, we reconstruct the full displacement profile ∆X(z i ,t) = ∑ t ≤tr (t ) at the considered vertical position, as schematically shown in Fig.2(a). Fitting of this profile to a sinusoidal function enables extracting the corresponding amplitude A(z i ) (see Fig.2(b)). We then estimate the local shear strain by numerically evaluating the partial derivative of the displacement field along the vertical direction γ(z) = ∂ A 0 (z)/∂ z. During an oscillatory shear experiment, we perform repeated scans at different times to monitor a possible time evolution of the deformation profile (further detains in the ESI). . The local shear strain can be then estimated as γ(z) = ∂ z A(z), which becomes z-independent for the example considered in a).

Microscopic non-affine displacements
The emergence of both reversible and irreversible microscopic non-affine displacements has been reported in simulations and experiments on sheared soft solids, and it has been related to the transition from linear to non-linear mechanical response, as well as to their yielding or flow behavior 17,31,48,49 . In such cases, even if the macroscopic displacement field ∆X remains homogeneous (see Eq.2), the actual displacement field: includes also a non-affine contribution ∆x N.A. due to rearrangements occurring on the scale of the elementary building blocks of the material. In general, ∆x N.A. is inhomogeneous in the shearvorticity plane, and it is not necessarily oriented along the shear direction. In the rest of Sect.2, after briefly reviewing the working principle of DDM, we discuss how it can be adapted to capture and quantify irreversible non-affine displacements in a periodically sheared sample.

Differential Dynamic Microscopy
The first step of the DDM analysis is the calculation of the difference between two images acquired at times t and t + ∆t, respectively (see Fig.3a-c). As proposed in Ref. 50 , we then multiply the so-obtained image difference with a windowing function before calculating its Fourier power spectrum and successively average over all the results obtained with the same time delay ∆t but different initial times t. The result of this average is the image structure function (Fig.3d) whereÎ(q,t) = FFT{I(x,t)} is the Fast Fourier Transform of the image at time t, q is the wavevector, and the symbol · indicates an average over the initial time t. This calculation is iterated for all the time delays ∆t that are needed to properly follow the sample relaxation. The interpretation of the results is based on the fact that D(q, ∆t) can be written in terms of the real part f (q, ∆t) of the intermediate scattering function (ISF), which is the quantity encoding information on the relaxation of the density Fourier modes, and whose squared amplitude is typically measured with DLS experiments: 21 The amplitude term a(q) is determined by the scattering properties of the sample and by the microscope transfer function, whereas the additive term b(q) accounts for the noise in the de-tection chain 20 . To extract the ISF from the image structure function, a(q) and b(q) need to be determined. Since f (q, ∆t → 0) = 1, b(q) could be in principle estimated as the limit for small ∆t of D(q, ∆t). In practice, we fit a quadratic function to the image structure function over a small-time interval close to the origin (typically including the first five delay times) and we estimate b(q) as the intercept of the best fitting curve 51 . Similarly, since f (q, ∆t → ∞) = 0, one could think to extract a(q) from the longtime plateau of D(q, ∆t). Unfortunately, because of the finite acquisition time set by the experimental constraints, in many cases we are not able to capture the full relaxation of f (q, ∆t) for all wavevectors. Consequently, we decided to adopt a different strategy 51 . We identify a background image I 0 (x) incorporating all the static, sample-independent contributions to the images (see ESI for detail) and we estimate the amplitude a(q) directly from the power spectrum of the raw images, corrected for the background Once a(q) and b(q) are known, Eq.9 can be inverted to obtain the ISF. The sample dynamics can be then extracted, for example, by assuming a specific model for the statistics of the particle displacements and fitting it to the ISF. For a dispersion of identical, non-interacting Brownian particles in a Newtonian fluid, the ISF is azimuthally symmetric and is given by f (q, ∆t) = e −D 0 q 2 ∆t , where D 0 is the diffusion coefficient of the particles. In this case, D 0 can be estimated as a coefficient of the q-dependent relaxation rate Γ(q) = D 0 q 2 , obtained by fitting an exponential function e −Γ(q)∆t to the ISF for different wavevectors (see Fig.3(e-f)). More generally, at least in the case of identical, non-interacting particles, the IFS can be expressed in terms of the Fourier transform of the probability distribution function of particle displacements P(r, ∆t), i.e. the self part of the van Hove function 52 If the dynamics is anisotropic, it can be worth considering separately the behavior of the ISF along specific directions in Fourier space. In this case, the azimuthal average of f (q, ∆t) in the qplane can be restricted to narrow angular sectors of aperture ∆θ oriented along the direction of interest 34 . Under the same conditions under which Eq.10 holds, we can write which enables estimating the mean square displacement along a selected direction (say, the x-axis) directly from the low-q behavior of the ISF along the corresponding direction in the Fourier space.

Echo-Differential Dynamic Microscopy (Echo-DDM)
In the presence of an imposed oscillatory deformation, an effective strategy to single out shear-induced irreversible rearrangements is to adopt a stroboscopic or echo approach, where the microscopic configuration of the system is probed at a sampling frequency that matches the oscillation frequency. This approach to shear-induced dynamics has been previously exploited in combination with light scattering, microscopy or ultrasonic echog-raphy 10,12,17,53,54 . In the case of microscopy, one considers sequences of echo images I n (x) taken at times t n = t 0 + nT , where t 0 is a time offset, T is the oscillation period and n is an integer number comprised between 1 and the total number of cycles N. This choice ensures that each frame is captured in correspondence to the same realization of the macroscopic displacement field and, as a consequence, any difference between the images can be attributed to the occurrence of non-affine, irreversible displacements within the sample. Since DDM is a very sensitive probe of the microscopic dynamics, one could think to use it directly to extract the effective, shear-induced dynamics present in the echo images. However, testing this idea requires first minimizing a number of potential artifacts which, if not adequately accounted for, can introduce spurious contributions to the echo dynamics. These artifacts can arise for example from drifts or flows due to an imperfect confinement of the sample or from a mismatch between the oscillation frequency and the sampling frequency. In the following section, we present a pre-processing procedure aimed at minimizing such effects and enabling robust DDM analysis.

Correction of real and apparent drifts
As observed in previous work 12 , the imperfect synchronization between deformation and acquisition can lead to an apparent, phase-dependent displacement in the echo images, possibly impairing the result of the experiment (see Fig.4). A relatively sim- Fig. 4 Apparent displacement due to imperfect synchronization between deformation and acquisition. In each panel, the continuous blue line represents the displacement along the shear direction of the observed plane as a function of time ∆X(t) = sin (2πft), the circles indicate the position at the sampling times t n = (ϕ i /2π + n)/ν, where ν = f + ∆ν is the sampling frequency and n is a positive integer number. Circles of different colors correspond to different initial phases ϕ i . In each panel, a different value of the frequency mismatch ∆ν is considered ∆ν/f = 0, 1%, and 10%, respectively. If ∆ν = 0 each sampling series always captures the system in the same position. When ∆ν = 0 instead, an apparent displacement is observed in each sampling series, the velocity of the apparent displacement depending on both ∆ν and on the initial phase ϕ i .
ple solution for obtaining a good synchronization would be to trigger the camera acquisition with a signal obtained from the shear-cell, proportional to the displacement imposed on the upper plate. Unfortunately, our shear-cell does not provide such a real-time output signal 40 , which forced us to act both at the instrumental level and at the image processing level. The proposed solution encompasses the following steps: accurate characterization of the oscillation frequency of the shear-cell, adoption of the best-matching image acquisition frequency, measurement of the residual drift in the collected images, and correction of the residual drift via image registration. We now describe these corrective steps.
Characterization of the oscillation frequency. In our setup, the shear-cell and the camera are controlled by two different computers. The first step is to characterize the frequency of the shear-cell in terms of the clock of the camera's computer. We observed that the actual frequency of oscillation is systematically different from the nominal one in an amplitude-dependent fashion, so that the actual frequency f turns out to be a function of the two input parameters f(f 0 , A 0 ). This function for f 0 = 1 Hz, A 0 ∈ (10 µm, 1 mm) is reported in ESI (Fig.S2). The relative difference between the nominal and actual oscillation frequency (f − f 0 )/f 0 is typically or the order of 0.1%. This difference turns out to be significant, especially for long experiments lasting some hundreds of oscillations. The experiment-to-experiment variability is estimated with the root mean square deviation measured for fixed input parameters (f 0 , A 0 ),δ f ∼ 1 · 10 −5 Hz (details of the frequency characterization are given in ESI).
Frequency matching: To increase the statistics while keeping the duration of the experiments fixed, it is convenient to acquire more than one image per period. Thus the acquisition frequency ν has to meet the condition ν = m ϕ f. The frequency ν is set indirectly by setting the delay time between two consecutive images ∆t 0 , that can be determined with a precision of 0.01 ms. This introduces a rounding error δ ν round /f ∼ 10 −5 of the same order of magnitude as the intrinsic shear frequency variability (see ESI for details). The parameter m ϕ has been chosen in order to minimize the rounding error. Since different shear amplitudes lead to a different actual frequency, the number of images per period may vary between experiments at different amplitude, taking values in the interval m ϕ ∈ [4,8].
Residual drift and registration: The residual frequency mismatch results in a phase-dependent apparent drift (Fig.4) that is quite evident in long acquisitions. By means of the same image cross-correlation algorithm used for the determination of the displacement field (see Sect. 2.3.1), we characterize the residual drift, and we remove its effect with a custom sub-pixel registration algorithm implemented in the Fourier space and described in detail in ESI. In the case of fluid samples, actual drifts, corresponding to macroscopic flows, could also occur as a consequence of the loose lateral confinement. With the aforementioned registration algorithm, we can effectively remove these global contributions as well.

Summary of the echo strategy
In summary, our optimized strategy consists in acquiring a timed image sequence I(x,t) including exactly m ϕ images per period. For each experiment at a given amplitude, we have m ϕ different stacks of echo images, that we label with their initial phase ϕ m = 2πm/m ϕ , m = 0, 1, ..., (m ϕ − 1) I n (x|ϕ m ) = I(x, (n + ϕ m /2π)T ). (12) Registration and computation of the dynamic structure function are then sequentially performed on each stack separately by working at constant m and calculating the image structure function for various delay cycles indexed by n. Finally, the dynamic structure functions obtained for different initial phases are averaged to give a single dynamic structure function D(q, ∆t), representative of the average response of the sample within the cycle.

Echo-Particle Tracking
As a cross-check of the results obtained with echo-DDM, we also perform a direct space analysis of the same image sequences, at least in those cases where the size of the tracer particles is large enough to enable them to be tracked individually. For each acquisition, echo image sequences are extracted and pre-processed as described in Sect. 2.4.4. Images in each echo sequence are convoluted with a Wiener kernel (with standard deviation of 2µm) to reduce noise and then analyzed with a customized version of the MATLAB particle-tracking code developed by the Kilfoil group at the University of Massachusetts 55 and available at https://github.com/dsseara/microrheology. We extract the trajectories [x (i) (t), y (i) (t)] of all particles (typically ∼ 100) in the field of view, and we evaluate the corresponding displacement probability distribution functions (PDF) P s (∆x, ∆t) and P v (∆y, ∆t), along the shear (x) and the vorticity (y) direction, respectively In the above expressions, the average is performed over all the initial times t, all the tracked particles i and all the echo sequences corresponding to different initial phases. From the same data, we also calculate the mean square particle displacement (MSD) along the shear and the vorticity direction, respectively

Measurement of the deformation profiles
For each of the four samples, we perform a detailed characterization of the deformation profile which is produced when an oscillatory shear deformation of different amplitude is imposed.

Carbopol
In Fig.5(a), the results of vertical scans of the Carbopol sample for different shear amplitudes in the interval γ 0 ∈ [10, 100]% are shown. For each amplitude (labeled by a different color) two successive full scans are performed. The excellent overlap between these pairs of measurements indicates that the observed deformation profiles do not significantly evolve over time and do not depend on the shear history. For all the imposed shear ampli-tudes, the displacement profile is almost perfectly linear. Fitting of each profile with Eq.4 enables estimating the corresponding deformation-loss-factor c = γ/γ 0 , which quantifies the discrepancy between the imposed strain amplitude γ 0 and the measured one γ. In the case of Carbopol, the deformation-loss-factor is almost constant c 0.9 (see Fig.5(c)), as can also be seen from the collapse of all the data onto a single master curve (panel (b)), which is obtained by normalizing each displacement profile by the displacement A 0 imposed on the top plate of the shear-cell. The fact that the displacement profiles are linear for all the investigated strain amplitudes indicates that the material response remains homogeneous even in the transition region. This is a non-trivial result, as the occurrence of inhomogeneous weakening, accompanied by the formation of shear bands, is often observed for both hard materials and for other yield stress fluids when driven out of the linear regime 41,44,56 .

Sylgard
For this sample, representative of an ideal elastomer, with nearly Hookean solid response, we performed three successive scans for each amplitude in the range γ 0 ∈ [1%, 70%]. Compared to Carbopol, the measured shear strain is in even better agreement with the imposed one (c 0.98), as shown in Fig.5(c). The absence of wall slip is indicative of the strong adhesion between the gel sample and the confining plates, which was obtained by maintaining the sample in contact with the glass slides during the curing process.

Silicon oil
Very similar results are also obtained for the Newtonian sample. The reconstructed deformation profile is homogeneous, and c 0.98 (see Fig.5(c)).

Worm-like micelles
Experiments performed with this sample reveal the presence of significant wall slip, which becomes increasingly stronger at large shear amplitudes, as indicated by the marked decrease of the deformation loss-factor as a function of the imposed deformation amplitude (yellow line in Fig.5(c)). Interestingly, for each amplitude, the reconstructed deformation profile is still homogeneous, characterized by a well-defined shear strain γ < γ 0 (see ESI for more details).

Echo-DDM
Once a detailed characterization of the sample deformation field is available, we proceed to quantify the shear-induced irreversible dynamics with the echo-DDM approach introduced in Sect. 2.4.2.
As pointed out in Sect. 2.4.3, a key requirement for the accurate determination of the non-affine, irreversible tracer displacements in the sample is the absence of spurious contributions (such as drifts or external vibrations) that could mask the genuine signal. Small imperfections or asymmetries in the setup may lead to a stress unbalance which, in the absence of an effective lateral confinement, can induce a net, albeit typically very small, sample flow. The relevance of this effect is found to crucially depend on the rheological properties of the material and it is particularly severe for low or no yield stress samples (see ESI for more details). As shear-induced flows are often not perfectly homogeneous, they cannot be completely removed with the registration procedure described in Sect. 2.4.3. As a matter of fact, for both the silicon oil and the WLM samples, the residual effect of shear-induced flows are found be large enough to dominate the microscopic dynamics (see ESI ), suggesting that our shear geometry is not optimal to measure irreversible tracer displacements in fluid-like samples. We stress that the effect under discussion has no impact on the results for the deformation profiles reported in Sect. 3.1, as the displacement per period associated with shearinduced flows is typically at least two orders of magnitude smaller than the amplitude of the affine displacement (see Fig.S5).
We thus focus on the two samples that display a solid-like behavior at rest: Sylgard and Carbopol. For these samples, all the echo acquisitions are performed at the center of the cell: z = h/2. This allows to neglect boundary effects (including the small wall slip) since the distance between the focal plane and the sample edges is way larger than the depth-of-focus h/2 L f . In principle, it would be also possible to perform measurements in the vicinity of the moving wall to study in more detail the plastic activity triggered by the interaction of the sample with the boundaries and, more generally, the role of surface boundary conditions for the flow of soft materials 57 .

Sylgard
This sample displays an elastic, reversible mechanical response over the whole explored range of shear amplitudes (up to γ = 80%) 35 and, as demonstrated above, it is firmly attached to the confining plates. Based on these premises, we do not expect to observe any genuine echo-DDM signal associated to irreversible non-affine displacements. As a consequence, these experiments can be thought of as a sort of calibration of the instrumental "background noise". For Sylgard, the ISFs obtained from an echo-DDM experiment with imposed shear amplitude γ 0 = 75% are shown in Fig.6. As expected, the dynamics, along both the shear and the vorticity direction, is extremely slow and only the very first part of the decay of the ISFs can be observed. Nevertheless, when the ISFs are plotted as function of the rescaled delay time q∆t, a nice collapse of all the curves is observed (insets of Fig.6), suggesting the presence of a ballistic-like relaxation mechanism 20 . Assuming an exponential-like decay f (q, ∆t) e −vq∆t , we can get a rough, order-of-magnitude estimate of the characteristic associated velocity along the shear direction v s ∼ 3 · 10 −4 µm/s.
We attribute this ultra-slow dynamics to the presence of an apparent relative motion of the different planes contributing to the image, which is due to the combination of imperfect synchronization and axial velocity gradient (see ESI for details). Consistently with this interpretation, we found that v b compares fairly well with the characteristic velocity spread 2πδ νγ 0 L f ∼ 5 · 10 −4 expected in the presence of a synchronization error δ ν ∼ 10 −5 Hz (see Sect.2.4.3) and a finite depth-of-focus L f ∼ 10 µm. Moreover, we observe that the characteristic velocity related with the ballistic decay of the ISF along the vorticity direction is markedly slower (by about a factor of three, see inset of Fig.6(b)) compared to v s . This is consistent with the presence in the echo images of an apparent ballistic motion mainly directed along the shear direction and marginally contributing to the dynamics in the perpendicular direction due to the finite angular aperture of the bow-tie regions used to calculate f v (q, ∆t) (see Sect.2.4.1).
These results show that the residual effects of imperfect syn-chronization of the echo dynamics are rather small, corresponding to characteristic displacements of less than 1 nm per cycle. As it will be explicitly demonstrated in the next subsection, as soon as a genuine contribution to the dynamics, due to structural rearrangements, is present, the ballistic term discussed in this subsection can be safely neglected.

Carbopol
This sample is a yield stress fluid. As such, it exhibits a transition from solid to liquid, as the amplitude of the applied shear strain exceeds a certain threshold. This threshold can be broadly identified with the cross-over point of G and G in a dynamic strain sweep oscillatory rheological test, which we take here as an estimate of the sample yield strain 58 (see ESI for details). We thus expect the Carbopol sample to exhibit a dramatically different dynamics below and above the yield strain γ C 60%. While a systematic study of this transition is surely interesting, to assess the potential of our methodology we selected two different shear amplitudes γ, respectively below (γ 0 = 20%) and above (γ 0 = 80%) the yield strain. In both cases, after the sudden application of the oscillatory shear, the sample exhibited a transient, whose effects were made negligible by waiting 120 oscillation cycles before starting the image acquisition. Below the crossover point When the imposed amplitude is small (γ 0 = 20%), visual inspection of the echo image sequences reveals that most of the tracers do not move appreciably, while a few of them undergo rare, intermittent displacements (see supplementary movie M1). The occurrence of these isolated rearrangement events can be highlighted by considering instantaneous (Fig.7b) or integrated echo activity maps (Fig.7c). Instantaneous echo activity map are obtained as the square modulus of the difference between two consecutive images within an echo sequence |I n (x) − I n+1 (x)| 2 , while the integrated map corresponds to the time average of all the instantaneous maps in the sequence, and is a proxy of the cumulative rearrangements occurred in the sample during the whole experiment.
To be more quantitative, we apply echo-DDM by considering separately the dynamics along the shear and the vorticity direction. As it can be seen from Fig.8(a,b), for this amplitude we cannot capture the full relaxation of the ISFs within our observation window. Nevertheless, as shown in the inset of panels (a,b), all the ISFs obtained for q ∈ [0.15 − 2.2] µm −1 nicely collapse onto a single master curve when plotted as a function of the rescaled time delay q 2 ∆t. This indicates that even below yielding the system undergoes an (extremely slow) diffusivelike relaxation, which results from irreversible rearrangements. This is in agreement with the seminal observations made with DWS in Ref. 10 . Here, we have the advantage that we can characterize this irreversible dynamics as a function of the wavevector q on a well-defined plane in the sample, while being able to simultaneously measure the local strain. In order to measure the q-dependent relaxation rates, we thus fitted the phenomenological model f (q, ∆t) = [1 + ∆tΓ(q)] −1 , to the experimental ISFs. According to Eq.11, this model provides the Fickean scaling ∆x 2 (∆t) = 2D 0 ∆t for the 1D MSD when the DDM relax- ation rates exhibit a quadratic dependence Γ(q) = D 0 q 2 on the wavevector. This is exactly the scaling that we observe in both directions, with a moderate deviation only for the smallest values of q (Fig.8(c)). By fitting the relaxation rates we obtain the estimates D s = (2.2±0.2)·10 −4 µm 2 /s and D v = (2.6±0.2)·10 −4 µm 2 /s for the diffusion coefficient associated with the shear-induced effective dynamics of the tracers along the shear and the vorticity direction, respectively. Within the experimental error, the shearinduced dynamics is isotropic D s D v with values for the diffusion coefficients about one order of magnitude larger than the one associated with the thermal motion of the same tracers in Carbopol under stationary conditions (see ESI for details). Beyond being slower, the tracer dynamics in the sample at rest is also non-intermittent, which contrasts with the rare, intermittent rearrangements that can be spotted in the sample sheared below the yield strain. Interestingly, this small number of intermittently rearranging particles coexist with a larger number of slowly, yet persistently rearranging ones, with the result that the overall tracer dynamics remains diffusive. In these conditions, the determination of the diffusion coefficient is sensitive to the choice of the field of view, and it would be interesting to study the spatial correlation between rearrangements. A detailed study of these effects goes beyond the scope of the present study.
The results obtained with DDM on the echo image sequences are confirmed by PT analysis of the same sequences. The estimated PDF of the tracer displacements along the shear and the vorticity direction, is shown in Fig.9(a) and (b), respectively. In agreement with DDM analysis, we do not find any statistically significant difference between the shear and the vorticity direction, confirming that the dynamics is substantially isotropic. Interestingly, for all investigated time delays, the PDFs is not Gaussian, as one would expect for ordinary diffusion, but displays exponentiallike tails. Nevertheless, the MSDs in the shear and vorticity directions, reported in panel (c), display a clean linear dependence on the delay time, MSD s (∆t) 2D PT,v ∆t, MSD v (∆t) 2D PT,v ∆t, indicating Fickean diffusion. The value of the associated diffusion coefficient is found to be D PT,s = (1.14 ± 0.08) · 10 −4 µm 2 /s and D PT,v = (1.24 ± 0.1) · 10 −4 µm 2 /s respectively. The discrepancy (about a factor of 2) between the estimates for the diffusion coefficient obtained from DDM and PT can be probably attributed to (c) q-dependent relaxation rates Γ s (q) (orange circles) and Γ v (q) (blue circles), along the shear and vorticity directions, respectively. Large circles correspond to q-values whose ISFs are shown in panels (a-b), with the same colors. Dashed lines represent the best fits to the relaxation rates with a quadratic model Γ(q) = Dq 2 . (d-f) same of (a-c) for imposed shear amplitude γ 0 = 80%. In (f) the horizontal black dashed line marks the maximum relaxation rate that can be reliably measured, corresponding to the oscillation frequency ν = 1Hz . differences in the populations whose dynamics is captured by the two techniques (DDM sees all the particles and PT only a userdependent selection).

Above the crossover point
When the sample is subjected to a large amplitude shear deformation γ 0 = 80%, exceeding the crossover threshold γ C , a dramatic change in the echo dynamics is observed. The amplitude of tracer displacements increases by orders of magnitude and the shearinduced activity shows a much larger degree of spatial homogeneity and temporal uniformity (see supplementary movie M2). This can be also appreciated by inspecting the echo activity maps reported in Fig.7(e-f), outlining the presence of plastic, irreversible rearrangements distributed across the whole field of view. The speed-up of the dynamics is well captured by echo-DDM analysis. In contrast with the small deformation case, the obtained ISFs display an almost complete temporal relaxation for all the investigated q-values, in both the shear and the vorticity direction (see Fig.8(d-e)). By plotting the ISFs as a function of the rescaled time delay q 2 ∆t, an excellent collapse of the ISFs on a single master curve is observed, highlighting once again a diffusivelike dynamics (inset of panels (d-e)). As it can be appreciated from Fig.8(d-e), over the whole considered q-range, the decay of ISFs is very well captured by the same one-parameter fitting model f (q, ∆t) = [1 + ∆tΓ(q)] −1 used for the small deformation case. The obtained q-dependent relaxation rates show a quadratic scaling with q along both directions (Fig.8f). Remarkably, in this case the dynamics is no longer isotropic, displaying a faster relaxation along the shear direction. Indeed, fitting the experimental relaxation rates to Γ v (q) = D v q 2 and Γ s (q) = D s q 2 provides D v = (1.9 ± 0.3) · 10 −1 µm 2 s −1 and D s = (2.9 ± 0.3) · 10 −1 µm 2 s −1 , respectively. Comparing the diffusion coefficients for γ 0 = 20% and γ 0 = 80%, we observe a ∼ 1000-fold speed-up of the dynamics for large amplitude deformation, which is a strong signature of shear-induced diffusion. PT-based analysis of the same data fully corroborates the results obtained with echo-DDM, and in this case the agreement with DDM is even better than for small deformation. In particular, they exhibit a linear dependence on ∆t. From a linear fit MSD = 2D∆t of the MSDs along both the investigated directions (Fig.9f) we obtain D PT,v = (1.7 ± 0.2) · 10 −1 µm 2 s −1 and D PT,s = (2.8 ± 0.2) · 10 −1 µm 2 s −1 , respectively, in fairly good agreement with the results obtained from echo-DDM.

Conclusion
In this paper, we combined oscillatory deformation in a simple shear geometry with bright-field optical microscopy. For four different samples, we characterized the deformation field, finding stationary and homogeneous deformation fields. We observed the presence of partial slip, and we measured the actual local deformation. In echo measurements, we observed drifts of different origin: "apparent" drifts due to an imperfect synchronization between the shear and the acquisition frequencies, and "real" drifts due to unavoidable instrumental imperfections and the lack of horizontal confinement. The effect of drift on the microscopic dynamics is corrected with a custom sub-pixel registration algorithm. Eventually, the DDM algorithm is applied to echo movies. In the case of an "ideal" elastic solid, the applied procedure shows, as expected, the absence of shear-induced-dynamics. In the fluidlike samples, the dynamics is dominated by flows generated by small asymmetry in the set-up and by the lack of lateral confinement of the samples. Although our registration algorithm can significantly reduce the impact of these drifts, their residual effect is still large enough to prevent the observation of any additional shear-induced dynamics.
With Carbopol, we perform experiments by using two different strain amplitudes, below and above the yield strain. The smaller amplitude causes a spatially heterogeneous, temporally intermittent activity of the tracer particles, which results in extremely slow, substantially isotropic shear-induced diffusion. Above the yield strain, the shear-induced diffusive relaxation becomes anisotropic and very fast, with an observed speed-up of about 10 3 compared to the previous case. Our results with Carbopol represent the first clear observation of a truly diffusive dis-persion relation Γ(q) ∼ q 2 in fluids that are yielding during oscillatory experiments. Of note, this behavior contrasts with recent results obtained with x-ray photon correlation spectroscopy 13,59 in which a faster-than-diffusive relaxation (closer to the behavior observed in colloidal polycrystals 60 ) was found for a nanoemulsion and a colloidal gel; at the same time, our results agree with the original interpretation provided for dense emulsions 10 and colloidal glasses 11,17 , in both cases lacking an experimental determination of the dispersion relation.
Beyond providing reliable and quantitative information about the shear-induced microscopic rearrangements, our approach also enables characterizing the deformation profiles across the sample and, in turn, spotting slippage, shear bands or other unwanted effects. This is made possible by the partial coherence of the illuminating light, which not only helps in identifying the presence of the above-mentioned phenomena, but also allows imaging a selected slice of the sample. In this work, we imaged its central part, where the deformation is homogeneous, by excluding the sample boundaries where deviations from homogeneity (mainly due to wall slip) possibly occur. In case of shear-banding, the z-selectivity of the microscope can be exploited to select zones where the macroscopic deformation field is approximately homogeneous, an important advantage over traditional single and multiple light scattering approaches. For this reason, our technique appears promising for the study of the plastic dynamics in the presence of shear banding, the interplay between non-affine displacements and shear banding representing a very interesting, yet not much explored topic.. Another advantage of our approach is that it makes it possible to combine real-space pre-processing (drift correction) and Fourier-space analysis (Echo-DDM), which leads to a considerable simplification of the data interpretation.
Extensions of our approach to other imaging modes (e.g. confocal fluorescence microscopy) are easily within reach, as well as combination with a commercial rheometer provided that the sample remains optically accessible. Finally, the usual advantages (and disadvantages) of DDM over PT hold, including the capability of DDM to probe the dynamics of sub-resolution particles 2 .

Conflicts of interest
There are no conflicts to declare.