Quantitative earthquake-like statistical properties of the flow of soft materials below yield stress

The flow behavior of soft materials below the yield stress can be rich and is not fully understood. Here, we report shear-stress-induced reorganization of three-dimensional solid-like soft materials formed by closely packed nematic domains of surfactant micelles and a repulsive Wigner glass formed by anisotropic clay nano-discs having ionic interactions. The creep response of both the systems below the yield stress results in angular velocity fluctuations of the shearing plate showing large temporal burst-like events that resemble seismic foreshocks-aftershocks data measuring the ground motion during earthquake avalanches. We find that the statistical properties of the quake events inside such a burst map on to the scaling relations for magnitude and frequency distribution of earthquakes, given by Gutenberg-Richter and Omori laws, and follow a power-law distribution of the inter-occurrence waiting time. In situ polarized optical microscopy reveals that during these events the system self-organizes to a much stronger solid-like state.

E arthquakes, complex reorganization of earth crust caused mainly by the sudden release in energy from the sliding of geological faults, are probably the most severe natural phenomena to adversely affect human lives. To cut through the complexity of such a large-scale reorganization phenomenon, there have been significant recent efforts in mimicking earthquakes in controlled laboratory experiments by studying the deformation and failure in various solid materials under external loads. These include fracture of rock samples under compression [1][2][3][4][5][6] , fracture of artificial rock of sintered polystyrene beads 7 , compression of mesoporous silica ceramics 8 , avalanches in wood compression 9 and stick-slip instabilities under shear in a twodimensional assembly of polymer disks 10 . A recent study 11 demonstrates that charcoal samples damped with ethanol show avalanche events similar to earthquakes due to the internal stresses generated from ethanol evaporation. In all experiments on three-dimensional systems, the readout of the events is acoustic emission in the form of crackling noise caused by predominant irreversible deformations of the optically opaque solid samples with typical shear moduli $MPa or higher 12,13 . This raises the question of upscaling, i.e., how laws at the laboratory scale translate at the geological scale of a tens to hundreds of kilometers. Moreover, all these experiments are performed with a very high sampling rate ($100 KHz and above), whereas seismic vibrations of engineering significance occur at frequencies from <0.2 Hz to 20 Hz 14 . Furthermore, the Poisson ratio for the rock samples, porous materials, and soft wood, being significantly less than 0.5 15,16 , there is an overall volume contraction of the materials, which is counterfactual to earthquake effects on the Earth crust: most earthquakes are indeed double couples and thus conserve volume. Recent numerical simulation study of creep response of a yield-stress system of bubbles 17 indicates the possibility of observing earthquake-like statistics in soft matter systems without the need of ad hoc friction. There are a few reports of Gutenberg-Richter-like scaling laws-configurational free entropy changes in colloidal glass due to reorganizational events 18 ; probability distribution of the force drop in stick-slip motion between two polymer plates 19 and stress fluctuation in polymer network 20 . Despite these experimental and simulation studies, to our knowledge, there is no experimental study on soft, continuous, and disordered materials with solid-like yield stress of few Pa, that reports simultaneously all the three well-known scaling relations and statistics (discussed below) similar to earthquake avalanches or controlled laboratory experiments mentioned above. Such study is important, since their mechanical properties at the laboratory scale may be more adapted as analogs of the geological scales. Moreover, the mesoscopic domain structures in this broad class of materials can be easily probed, and can be reversibly controlled by an applied shear stress or an external field 21 . In general, similar tunability of domain structures is not accessible for conventional solids. In very few cases, the stress-induced evolution of domain structures have been observed in solids 13 , however, how the structural evolution change the overall shear modulus of the system is not clear. The major challenge to study reorganization events in soft yield-stress materials originates from the fact that they are much softer (shear moduli $tens of Pa, Supplementary Fig. 1) compared with those used in experiments mentioned earlier. Materials with such low values of shear moduli do not produce audible crackling noise under stress and, therefore, the statistical properties of reorganization events cannot be read out using acoustic emissions.
Here, using shear rheology along with polarized optical microscopy, we study the statistical properties of stress-induced reorganization in soft yield-stress materials. Accurate determination of the value of yield stress in soft and disordered materials is very challenging [22][23][24] . Recent experiments on a wide variety of soft materials 25 confirm the existence of a solid-like state below the yield stress consistent with the Herschel-Bulkley (HB) model. In our study, the estimated values (obtained by fitting the HB model to the experimental data) of the yield stress are of the order of a few tens of Pa (Fig. 1a). We apply perturbations in the form of shear stress and, since the system is incompressible, there is no volume strain in the sample. We directly measure and analyze the mechanical response of the sample by measuring the angular displacement (ϕ) of the shearing plate with a time resolution of 0.04 s. For constant stress well below the yield stress, the angular velocity ( _ ϕ) of the top plate (proportional to the shear rate) shows rapid temporal fluctuations around zero 26,27 . For the first time we   Supplementary Fig. 3 for the expanded version of the σ ¼ 2 Pa data), and c shows plots for four different d with same applied stress σ ¼ 2 Pa. find that, for small sample thicknesses (the gap d between shearing plates of a few tens of microns), the angular velocity fluctuations show burst-like events (Fig. 1b, c) persisting over about thousands of seconds (comparable with seismic foreshock to aftershock timescale), arising from slow building up and decay of angular velocity amplitudes. Such burst-like events closely resemble seismograph data during earthquakes resulting from failure and reorganization events of the earth crust over large length scales. In our system, these events indicate stress-induced reorganizations in the sample below yielding. Remarkably, the statistics of angular velocity fluctuations during these reorganization events follow well-known scaling relations given by Gutenberg-Richter law 28 , the direct and inverse Omori laws 29 and show power-law distributions of inter-occurrence waiting time 30 , qualitatively similar to earthquake avalanches.
In this study, we describe the creep behavior of two wellcharacterized, soft yield-stress materials: a dense nematic phase of rod-like surfactant micelles 31 with micron size randomly oriented nematic domains which are easily visible using polarized optical microscopy (POM), and a repulsive Wigner glass formed by a colloidal suspension of nanometer-sized laponite platelets 32 . The glassy state of laponite suspensions is formed due to long-range Coulomb repulsion between the charged platelets. Due to small domain structure, laponite system cannot be probed using POM. We discuss about the first system in detail and also touch upon the second one at the end.

Results
Burst phenomena in creep flow of CTAT nematic. Figure 1a shows the steady-state flow curves (shear stress vs shear rate) of the nematic phase formed by 39 wt% of surfactant Cetyltrimethylammonium Tosylate (CTAT) in water (see the Methods section). Below the yield point, the steady-state shear rate _ γ ¼ R d _ ϕ (R is the radius of the top plate) shows positive and negative values with a very small average magnitude jh_ γij $ 10 À4 s À1 , suggesting a solid-like state as also seen for other yield-stress fluids 24 . Despite complex temporal dynamics, the average steadystate shear rate shows very good agreement with the HB model for soft yield-stress materials: σ ¼ σ y þ αh _ γi n , with n ¼ 0:62, σ y ¼ 28 Pa for d ¼ 10 μm and n ¼ 0:81, σ y ¼ 29 Pa for d ¼ 25 μm. The value n < 1 indicates shear thinning, which is similar to the velocity-weakening nature of faults during earthquakes 33 .
We show _ ϕ as a function of time for different stress values (σ) as indicated in Fig. 1b for a sample thickness d ¼ 10 μm. For very small stress values (σ ¼ 0:01 Pa), no burst-like event is observed and _ ϕðtÞ fluctuates steadily in time. However, for σ ! 0:1 Pa and beyond, a burst occurs with time span $ 1000 s. Inside these bursts, the amplitude of _ ϕ fluctuations ($ 2 mrad s −1 ) (referred as "quakes" in this paper) is much larger than the average long-term angular velocity (h _ ϕ steady i $ 0:2 10 À4 mrad s −1 ), revealing strong stress-induced organization in the sample. The center of the burst gets shifted toward longer times as the sample thickness (d) increases (Fig. 1c). We also find that in the time window of the burst in _ ϕ, there are large jumps in the cumulative strain γðtÞ Fig. 4), signifying largescale stress-induced reorganization in the system. However, after the burst, γðtÞ varies very slowly with time, implying that in this temporal region the system enters into a state with much larger shear modulus (about two orders higher in magnitude; Supplementary Fig. 1). For higher values of stress (>30 Pa), γðtÞ increases rapidly with time. Interestingly, for small applied stress values, the long-term strain can even reverse course as a result of the strong internal reorganization in the system. Such anomalous flow behavior has been observed in the context of congested traffic flows, known as the Braess paradox 34 . One important feature of our systems is reversibility. After the burst-like event in _ ϕ, the system goes into a steady state over the remaining duration of the experiment (till 8000 s). However, if the sample is fluidized after the burst-like event by applying a stress larger than the yield stress (σ > σ y ) and then made to creep under a stress below the yielding (σ < σ y ), the system again shows a burst in the angular velocity ( Supplementary Fig. 5). This indicates that fluidization erases the history of the reorganization of the system.
Gutenberg-Richter law for angular velocity quakes. We next focus on the statistical properties of these angular velocity quakes inside these burst events. We estimate the kinetic energy associated with each quake as E ¼ _ ϕ 2 in units of mrad 2 s À2 . We find that the ensemble probability density function of kinetic energies PðEÞ shows a robust power-law decay over more than three decades in E ( Fig. 2a): PðEÞ $ 1=E ϵ , where ϵ ¼ 1 þ ðb=1:5Þ ¼ 1:6 ± 0:1; consistent with the Gutenberg-Richter law for the distribution of earthquake amplitudes 28 , and is close to the value of ϵ obtained in other experiments 2,6,10 . The value of the Gutenberg-Richter exponent b for our system (¼ 0:90 ± 0:15) is very close to the most accepted value for earthquakes (b $ 1), but slightly greater than the reported values for other fracture experiments (b $ 0:60) 8,9,11 . The exponent b obtained via the maximum likelihood method 8 for different low-energy cutoffs (E cutoff ) is very stable for our system (inset Fig. 2a). PðEÞ versus E corresponding to different applied stress values and sample thicknesses are shown in Supplementary  Fig. 6. We note that the relatively shorter range of power-law regime compared with other laboratory experiments cited earlier comes from the low-energy cutoff effects due to the background activity. The background activity is visible as the residual activity in the steady state (after 1500 s in Fig. 1b). This limitation in energy range is also seen in the probability density function with the mechanical energies of the continuously sheared granular matter 10 , with the energy released in the numerical simulations of a system of bubbles 17 , and in many other systems.
Statistics of inter-occurrence waiting time. The distributions of inter-occurrence waiting time of quakes for the three data sets corresponding to d ¼ 10 μm (having energy cutoff 2:07 10 À2 in units of mrad 2 s À2 ) follow a power-law behavior with exponent $1:5 (Fig. 2b). We show the robustness of this power-law by taking different energy cutoff values for the σ ¼ 2 Pa data (inset Fig. 2b).
Omori law. We now quantify the dynamics around the main shock of the burst. We perform a timescale analysis of the angular velocity _ ϕðtÞ using a wavelet transform, which is well suited for transient signals. The obtained wavelet coefficient at time τ and scale a reads: where Ψ stands for the conjugate of the mother wavelet, here chosen as the Morlet wavelet. The factor 1 ffiffi a p ensures the normalization of Ψ ð tÀτ a Þ. The Morlet wavelet has a dominant period T 0 , so that a daughter wavelet Ψ ð t a Þ will thus feature a dominant period T ¼ aT 0 . In the remaining of this paper, we shall use the corresponding time period T instead of the scale a, and study the wavelet coefficients C τ;T . Guided by Omori's law observed in earthquake bursts, Eq. (1) predicts Where τ s is the main-shock time. The power-law behavior of the wavelet coefficient near the singularity (t ¼ τ s ) reveals the finite-time singular power-law behavior of the angular velocity.
To check for a possible asymmetry on the left versus right side of the singularity time τ s , we use the two distinct mother wavelets Ψ L and Ψ R , where Ψ L is identical to the Morlet wavelet for t 0 and zero for t > 0, while Ψ R is chosen to be zero for t < 0 and identical to the Morlet wavelet for t ! 0. For all observed bursts, the wavelet coefficients are found to scale with the wavelet periods in a power-law fashion, almost over three decades of wavelet period with the value of the exponent α $ 1 as shown in Fig. 3a. We can also directly plot the rate of foreshocks and aftershocks as a function of time distance from the main shock to extract the Omori exponent (Fig. 3b). We find that Omori exponent of 1 fits the data very well for both fore-shock and after-shock regimes (Fig. 3b). For this system, the background amplitude during the foreshocks is greater than the one during aftershocks (so that the foreshocks E cutoff > aftershocks E cutoff ). After the background subtraction, the rate of foreshocks is lower than the rate of aftershocks for a given time distance from the main shock, similar to the rate asymmetry of real earthquakes. Only the data for σ ¼ 2 Pa with d ¼ 10 μm is plotted for clarity. For other stress values and sample thicknesses, the wavelet coefficients are plotted as a function of the wavelet period in Supplementary Fig. 7, where power-law behavior is obtained with α $ 1, similar to the statistics of the occurrence of aftershocks for the earthquakes 1,29,35 . While the inverse and direct Omori laws are recovered in their time dependence, there is an important difference with earthquakes: real foreshocks-aftershocks sequences are highly asymmetrical, the more so, the larger the main shock, with just a few events before the main shock (the number of foreshocks being independent of the size of the latter) and a variable number of aftershocks. This asymmetry reflects the productivity law, stating that the number of triggered events grows exponentially with the main-shock magnitude. This productivity law results from the extended nature of the rupture and the interactions within complex fault networks. It is thus naturally expected that below some magnitude threshold (depending on the specific productivity parameters) earthquakes feature less aftershocks than foreshocks.
Quantitative in situ rheomicroscopy of domain reorganization. Using the birefringent properties of the nematic phase of the surfactant, we capture the dynamics of nematic domains under stress by a home-built POM setup (see the Methods section). Even in the huge burst, the largest angular velocity quake corresponds to a displacement of the top plate by $1 μm, which is well below the resolution of the optical imaging system. However, we can probe the reorganization in the system by tracking the changes in the domain structure. We use transparent circular glass plates as top and bottom plates of the rheometer with a gap between them $150 μm. We capture images at a frame rate of 0.1 Hz with a CCD camera (Methods). The details of the setup are described in Supplementary Fig. 8. In Fig. 4b-e, we show the evolution of the intensity pattern as a function of time for an applied stress σ ¼ 2 Pa. Over time, these structures show a coarsening effect due to the orientational ordering of domains (Supplementary Movie 1). We find that when the applied stress is zero, no such evolution of intensity pattern is observed (Supplementary Movie 2), further confirming that such ordering in domains is due to the internal reorganization of the system under applied shear stress. Structure coarsening is quantified by calculating the spatial auto-correlation of intensity (Aðk; tÞ) in the wave vector (k) domain at different times t. For t < 2000 s, Aðk; tÞ decays with increasing k and t (Fig. 4f). After 2000 s, it does not evolve much with time. To bring out the largest length scale of the structures, we plot k À1 max as a function of time, where Aðk max ; tÞ ¼ 0:2 (Fig. 4g). It can be clearly seen that, after $ 2300 s, the structural reorganization is minimal. Another measure of this is the recently discovered computable information density (CID), which is proportional to the losslessly compressed image size 36 . CID is defined as CID LðtÞ=L, where L(t) is the size of the compressed image and L is the original size. All our original uncompressed "BMP" images are of the same size (L ¼ 602 KB).
We plot CID as a function of time, using L(t) the size of the images compressed to "PNG" format (Fig. 4h). We find that the image sizes drop until the system reaches a steady state after the burst. This indicates that the structures are getting bigger in size until the main shock appears. Since the POM captures the orientational order in the system, a larger domain structure ) before (using Ψ L as the analyzing wavelet; points on the left portion of the burst) and after (using Ψ R as the analyzing wavelet; points on the right portion of the burst) the main shock, as a function of the wavelet period T of the angular velocity for σ ¼ 2 Pa, d ¼ 10 μm (Fig. 1b). The slope suggests a singularity (the blue dashed line indicates a power-law decay with exponent α ¼ 1). The results for other bursts are shown in Supplementary Fig. 7. b The rate of quakes (N) before and after the main shock (τ s ¼ 594 s; having foreshocks E cutoff ¼ 2:31 10 À1 , aftershocks E cutoff ¼ 2:57 10 À2 ) of the same data is plotted versus time distance from the main shock (t 0 ¼ jt À τ s j). The counting time window Δt ¼ 30 s is taken to get good statistics in N. indicates a larger range of local orientational ordering in the system. This observation signifies that the domains form a jammed structure with the average nematic domain size $ 70 μm to counter the applied stress against any further rearrangements in the sample.
Flow of colloidal glass of laponite nanodiscs. To ascertain that the statistical properties of _ ϕ-quakes described so far are generic to other yield-stress solids, we have also studied a repulsive Wigner glass formed by 5 wt% laponite clay suspension. Laponite clay particles are disk-shaped having a diameter $ 20 nm and thickness 2 nm 37 . The particles have a weak positive charge along the rim but a net negative surface charge that makes the interparticle interactions long ranged 32 . Even without the application of any shear stress, these systems show considerable aging effects where viscosity can increase by orders of magnitude with waiting time 38 . An important issue while doing rheology with aging samples is to get a reproducible initial condition, because the process of sample transfer and loading in the rheometer involves shearing of the sample. We use the following protocol to get a reproducible initial condition 27 : after loading the sample on the rheometer plate, we apply a small constant shear rate of 0.1 s −1 for 1000 s. Then the sample is rejuvenated by applying a large stress of 50 Pa for 500 s. As soon as this process is over, we apply the desired stress (σ) for which we want to monitor the time dependence of the angular velocity ( _ ϕ). We summarize our experimental data of 5 wt% laponite clay suspension with σ ¼ 15 Pa and d ¼ 40 μm and 10 μm in Fig. 5. The angular velocity _ ϕ as a function of time is shown in Fig. 5a, where we again see burstlike events. In this case also, the quakes follow the Gutenberg-Richter law: PðEÞ $ 1=E ϵ ; ϵ $ 1:5 ( Fig. 5b) with Gutenberg-Richter exponent b $ 0:75. The waiting time distributions show exponent $ 1 (Fig. 5c). The wavelet coefficients are again found to scale with the wavelet periods in a power-law fashion (Eq. (2)) over three decades of wavelet period with the exponent α $ 1 as shown in Fig. 5d. Samples with different laponite concentrations show similar bursts, but with different quake amplitudes (Supplementary Fig. 12). For this system also, we can simply erase the memory by applying the pre-shear protocol (Supplementary Fig. 13). The laponite system is much more complicated due to the intrinsic aging of the sample at rest (see Supplementary Fig. 14 for experiments with different aging times). We found that the Gutenberg-Richter exponent b ¼ 0:75 ± 0:05 (summarizing all our results for Laponite) is very stable with respect to the different conditions imposed on the laponite sample here. However, after the primary burst, there are several secondary ones, unlike in the surfactant system. Such behavior probably arises from the complex structural evolution resulting from the competition between intrinsic and shearinduced aging effects. Secondary bursts appearing in d ¼ 10 μm data are analyzed, and they are also showing power-law behaviors ( Supplementary Figs. 10, 11). We would like to point out that even for this system, when the aging is sufficient to drive the system into a glassy state, the stress-induced yielding shows very similar behavior like the nematic phase of CTAT (see Fig. 1 of ref. 27 ). Laponite colloidal glasses do not show any birefringence property over the concentration range we have studied, and hence no imaging could be done.

Discussion
In our experiments, the observation of both positive and negative angular velocity fluctuations as a result of a unidirectional drive (constant applied stress) can be rationalized by the physical picture proposed by Nandi et al. 39 , based on a stochastic attachment and detachment dynamics of surfactant micelles/ laponite particles with the rotating top plate of the rheometer. As the sample deforms, micellar domains/particle clusters create transient bridges between the two plates that form and break intermittently, leading to a very jerky motion (strain) of the drive similar to stick-slip dynamics 10,40 , with the angular velocity showing large fluctuations (quakes) of both positive and negative signs. Within this physical picture, the average strainstrengthening behavior after the burst can be associated with jamming of such bridges that leads to an increase of the average long-time viscosity. The bridges constitute force paths, similar to the force chains in proposed analogy between granular media, tectonic deformations, and spinglasses 41,42 and to the stress field in model of fault network in the Earth crust 43,44 . The coarsening of nematic domains in the micellar system with time under an applied stress likely indicates the growth of such jammed regions in the system whose reorganization leads to a stronger solidlike state.
While not being perfect analogs, our experimental systems can be put in correspondence with the seismogenic crust via dimensional scaling of lengths, times, and mechanical properties, such as strain and stress. For this, we introduce the scaling factors L sf ¼ L=L Ã and T sf ¼ T=T Ã relating the length L Ã and timescales T Ã of the model to those, L and T, in the Earth crust. The spatial range over which tectonic shear stresses occur in the Earth crust due to tectonic motions is of the order of 100 km. The width of the gap over which the shear stress is applied in the experiment is $10 μm. Thus L sf $ 10 10 . The typical duration of earthquake clusters (including foreshocks-mainshock-aftershocks) is of the order of months to years. The typical time span of a burst in the laboratory experiment is about 1000 s, thus we have T sf $ 10 5 . From these two scaling ratios L sf $ 10 10 and T sf $ 10 5 , we proceed to derive the dimensional scaling laws of all other relevant physical quantities, allowing to map our laboratory experiment onto the Earth crust.
An earthquake in the crust is associated with the driving shear stress overpassing the frictional force along faults. The frictional force is equal to the normal stress, controlled by the lithospheric pressure, times the friction coefficient, in the range 0.1-0.6. Hence, the scale of the stresses at which fault sliding is triggered is governed by ρgh, where ρ is the average density of the Earth crust, g is the acceleration of gravity, and h is the width of the seismogenic crust. As ρ is about three times that of our suspensions, and g is the same, this predicts that the shear modulus in our experiments should be about 1=L sf that of the Earth crust (roughly 10 10 to 10 12 Pa), i.e., 1-100 Pa in the experiments, which is of the same order as in our laboratory observations. The typical rupture sliding velocity during an earthquake is $ 1 m s −1 , which should be L sf =T sf $ 10 5 times the velocity of local peaks of the burst in our experiments. Indeed, the typical velocity of the shear motion at a local peak is $ R _ ϕ $1 cm 1 mrad s À1 ¼ 10 À5 m s −1 . Lastly, since strain is dimensionless, the strain rate of $ 0:1 s À1 in our sample at a peak during a burst should correspond to T sf ¼ 10 5 times the strain rate during an earthquake, which is typically 1 m of slip over a fault length of 100 km occurring in 10 s of slip, which yields the strain rate $ 10 À6 , which is in perfect agreement with the prediction from the scaling laws.
We have shown a good semi-quantitative description of our experimental data with Gutenberg-Richter law, Omori scaling law, and the power-law distribution of inter-occurrence waiting time over a wide range of stress values and system sizes. The robustness of the scaling exponents observed in our experiments and their similarity with the values obtained for other systems having orders of magnitude different moduli suggest an underlying universal mechanism of stress-induced reorganization in different materials via the general process of avalanches triggered by threshold dynamics 45 . Our experiments on yield-stress materials may help to correlate such scaling laws (and especially deviations from them) directly to the microscopic deformations in the system.
The reversal of the low-frequency angular velocity is reminiscent of an analog in congested traffic flows, known as the Braess paradox 34 : when a new link is added to a random network (a new cluster or bridge forms), the capacity can be augmented or decreased with equal probability 46 . The relevance of this mechanism for mechanical networks has been exemplified by Cohen and Horowitz 47 . In addition to the quantitative analogies with the inverse and direct Omori laws, one can observe that the angular velocity exhibit seismogram-like structure having many quakes in it. This parallels the succession of repetitive seismic bursts, accelerating with time, and increased low-frequency seismic noise, which was observed for one of the best-recorded large earthquakes to date, the 1999 7.6 Mw Izmit (Turkey) earthquake 48 .
It is tempting to interpret the log-periodic variation of wavelet coefficients seen near the singularity (Fig. 3). Such behavior has been observed in case of earthquakes 35 , material failures 49 , finance, and population dynamics 50 , indicating a discrete-scale invariance property of the system as opposed to continuous-scale invariance given by ordinary power laws. However, such study is outside the scope of this paper. We hope that our work will motivate further experimental and theoretical studies in such subtle but fascinating stress-induced reorganizations in a wide range of materials.

Methods
Procedure of rheology experiments. All the rheology experiments are carried out using a MCR 102 stress-controlled rheometer (Anton Paar, Austria) fitted with PT temperature controller using 25-mm steel parallel-plate (P-P) geometry. A humidity chamber is used to minimize the evaporation of water from the sample during the rheology experiments. The bottom plate of the instrument is fixed, and the top plate is connected to the torque transducer as shown in Supplementary  Fig. 2. The minimum angular resolution of this instrument is 0.01 μrad. The details of the experimental setup are same as in ref. 26 . Due care is taken to cutoff all external sources of vibrations using vibration isolation arrangements. For Rheo-POM (rheology-POM) measurement, steel plates are replaced by glass plates with the polarizer axis along the vorticity direction ( Supplementary Fig. 8). The images were grabbed at a 5 mm interior from the edge of the top plate, with an eight-bit color CCD camera (Lumenera, 0.75C, 640 480 pixels) fitted with a microscope (objective: M Plan Apo 5X, N.A. 0.14) and analyzed in ImageJ.
Preparation of CTAT-water system. Cetyltrimethylammonium Tosylate (CTAT)-water samples are prepared by dissolving known amount of CTAT (Sigma Aldrich) in ultrapure deionized double-distilled water (Millipore) at pH 10 in 10ml vials, and are kept sealed for equilibration for two weeks at 60°C.
Preparation of laponite-water system. The laponite powder as procured from Southern Clay Products Inc. is mixed with deionized water for 10 min with a magnetic stirrer until a clear solution is obtained. The solution is then passed through a 0.45 μm pore size filter to obtain the final sample. During the filtration process, the sample is under very high stress/strain rate, and hence we define the zero of the aging time for the sample after the filtration process is over.

Data availability
All the data that support the findings of this study are available from the corresponding author upon reasonable request.