Thibaut
Divoux
*n,
Elisabeth
Agoritsas
j,
Stefano
Aime
a,
Catherine
Barentin
b,
Jean-Louis
Barrat
f,
Roberto
Benzi
l,
Ludovic
Berthier
d,
Dapeng
Bi
y,
Giulio
Biroli
f,
Daniel
Bonn
c,
Philippe
Bourrianne
t,
Mehdi
Bouzid
s,
Emanuela
Del Gado
p,
Hélène
Delanoë-Ayari
b,
Kasra
Farain
c,
Suzanne
Fielding
af,
Matthias
Fuchs
u,
Jasper
van der Gucht
r,
Silke
Henkes
ag,
Maziyar
Jalaal
z,
Yogesh M.
Joshi
w,
Anaël
Lemaître
i,
Robert L.
Leheny
k,
Sébastien
Manneville
n,
Kirsten
Martens
ab,
Wilson C. K.
Poon
e,
Marko
Popović
o,
Itamar
Procaccia
hae,
Laurence
Ramos
d,
James A.
Richards
e,
Simon
Rogers
x,
Saverio
Rossi
g,
Mauro
Sbragaglia
l,
Gilles
Tarjus
g,
Federico
Toschi
mad,
Véronique
Trappe
v,
Jan
Vermant
q,
Matthieu
Wyart
j,
Francesco
Zamponi
fac and
Davoud
Zare
aaah
aMolecular, Macromolecular Chemistry, and Materials, ESPCI Paris, Paris, France
bUniv. de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
cSoft Matter Group, van der Waals-Zeeman Institute, University of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands
dLaboratoire Charles Coulomb (L2C), Université Montpellier, CNRS, Montpellier, France
eSUPA and the School of Physics and Astronomy, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, UK
fLaboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, F-75005 Paris, France
gLPTMC, CNRS-UMR 7600, Sorbonne Université, 4 Pl. Jussieu, F-75005 Paris, France
hDept. of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel
iNavier, École des Ponts, Univ Gustave Eiffel, CNRS, Marne-la-Vallée, France
jDepartment of Quantum Matter Physics (DQMP), University of Geneva, Quai Ernest-Ansermet 24, CH-1211 Geneva, Switzerland
kDepartment of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA
lDepartment of Physics & INFN, Tor Vergata University of Rome, Via della Ricerca Scientifica 1, 00133, Rome, Italy
mDepartment of Applied Physics and Science Education, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
nENSL, CNRS, Laboratoire de physique, F-69342 Lyon, France. E-mail: Thibaut.Divoux@ens-lyon.fr
oMax Planck Institute for the Physics of Complex Systems, Nöthnitzer Str.38, 01187 Dresden, Germany
pGeorgetown University, Department of Physics, Institute for Soft Matter Synthesis and Metrology, Washington, DC, USA
qDepartment of Materials, ETH Zürich, Vladimir Prelog Weg 5, 8032 Zürich, Switzerland
rPhysical Chemistry and Soft Matter, Wageningen University & Research, Stippeneng 4, 6708WE Wageningen, The Netherlands
sUniv. Grenoble Alpes, CNRS, Grenoble INP, 3SR, F-38000 Grenoble, France
tPMMH, CNRS, ESPCI Paris, Université PSL, Sorbonne Université, Université Paris Cité, Paris, France
uFachbereich Physik, Universität Konstanz, 78457 Konstanz, Germany
vDepartment of Physics, University of Fribourg, Chemin du Musée 3, Fribourg 1700, Switzerland
wDepartment of Chemical Engineering, Indian Institute of Technology, Kanpur 208016, Uttar Pradesh, India
xDepartment of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
yDepartment of Physics, Northeastern University, Boston, MA 02115, USA
zInstitute of Physics, University of Amsterdam, Science Park 904, Amsterdam, The Netherlands
aaFonterra Research and Development Centre, Dairy Farm Road, Fitzherbert, Palmerston North 4442, New Zealand
abUniv. Grenoble Alpes, CNRS, LIPhy, 38000 Grenoble, France
acDipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, 00185 Rome, Italy
adCNR-IAC, Via dei Taurini 19, 00185 Rome, Italy
aeSino-Europe Complex Science Center, School of Mathematics, North University of China, Shanxi, Taiyuan 030051, China
afDepartment of Physics, Durham University, South Road, Durham DH1 3LE, UK
agLorentz Institute, Leiden University, 2300 RA Leiden, The Netherlands
ahNestlé Institute of Food Sciences, Nestlé Research, Vers Chez les Blancs, Lausanne, Switzerland
First published on 2nd July 2024
Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.
When subjected to external forces, soft amorphous solids undergo a solid-to-liquid transition known as the yielding transition. The characteristics of this transition vary significantly depending on factors such as material properties, system preparation, the geometry confining the material, and boundary conditions. Conversely, following flow cessation, soft amorphous materials undergo a liquid-to-solid transition in diverse ways. Notably, memory effects from previous flow history can significantly impact the mechanical properties of the sample. This phenomenon is observed in various materials, including gels, soft glasses used in food products, cementitious items, and biological materials, potentially with important consequences on their practical use.
A general discussion of both yielding and memory phenomena took place during a one-week workshop at the Lorentz Center in Leiden from January 9–13, 2023. The workshop aimed at reviewing current knowledge in the field and at identifying upcoming challenges. This manuscript summarizes the enriching discussions from the workshop. It is organized as follows: Section 1 delves into the shear-induced yielding transition in soft amorphous materials, discussing its ductile or brittle nature, and the associated local scenarios. Section 2 focuses on memory effects in soft amorphous solids imprinted through shear history, exploring their impact on sample properties and relevant applications. Finally, Sections 3 and 4 provide insights into two specific examples: biological materials as soft glasses and material design through the utilization of shear history and memory effects. The manuscript concludes with a list of pressing questions that should guide the community's research agenda in the coming years.
In soft matter, where the range and hierarchy of microstructural motifs and time scales can be extreme and become easily accessible, material failure (under an imposed deformation or stress) is a complex process that may manifest itself via a rich phenomenology. Here, we focus on the behavior of soft amorphous materials under deformation: they yield and eventually flow, with yielding being often preceded by, and associated with, mechanical and flow instabilities, whose prominence and persistence depend on their mechanical or rheological history, as well as on the imposed stress or deformation rate.2 At low deformation rates, the yielding of soft amorphous solids emerges from avalanches of localized plastic events,3–5 which are activated by the externally applied stress or strain,6–8 and is often associated with extended and pronounced shear localization or banding, where part of the material remains “stuck” while the rest is already flowing.9,10
At first sight, yielding seems obviously a manifestation of the ductility of soft materials, and may be associated with a non-monotonic stress vs. strain response as fluidization takes place.11,12 Ultimately, the very same soft materials that yield may fracture,13–17 therefore raising the question of what factors and processes control, over the relevant range of length scales and time scales, the emergence of brittle versus ductile behavior in yield stress fluids and soft materials.
Fig. 1 Typical stress–strain (or time) curves measured under shear (a) in a zirconium-based metallic glass alloy (Vitreloy 1) at different temperatures (extracted from ref. 11), (b) in a model waxy crude oil system at a fixed shear rate ( = 2 s−1) for various waiting times elapsed since the rejuvenation of the sample (extracted from ref. 12), and (c) in an acid-induced casein gel (extracted from ref. 17). In (c), the inset shows the same data in semilogarithmic scale, and the Taylor–Couette cell in which the experiments are performed. The upper panel shows the images of the side view of the Couette cell at different strains recorded simultaneously to the experiment reported in the main graph; macroscopic cracks are visible beyond the stress maximum. |
Regarding experimental protocols, laboratory mechanical testing of amorphous materials in the soft matter community most often takes place in a rheometer, where a simple shear is applied in either parallel-plate or cone-and-plate geometry. However, stress–strain curves akin to the one shown in Fig. 1(b) may also be generated under compression or extension using a “universal testing machine” (UTM), or biaxial and triaxial apparatuses. A crucial issue is, therefore, whether there are differences regarding the response of a given material between different geometries. If the specimen is brittle, using compression in UTM, we expect the bulk compressive stress to rise until shear bands (or fractures) occur, usually at 45° from the compression direction as sketched in Fig. 2(a). Strain localization along a 45° plane closer to the compressing surfaces (where there are stress inhomogeneities) is suppressed by the rigid boundary itself. On the other hand, in a rheometer, the thickness of the sample may be such that stress heterogeneities induced by, e.g., irregularities in the rheometer plates, propagate throughout the sample [see Fig. 2(b)]. Furthermore, the geometry does not prevent shear deformation near the boundaries (as in a compressive test). Under these circumstances, even a sample that is structurally homogeneous may display stress inhomogeneities, and, therefore, yield in a way that may appear to be ductile in a rheometer while showing brittle yielding in a UTM [see Fig. 2(c)]. Therefore, caution must be taken when interpreting the macroscopic mechanical response obtained from rheological tests, whether in a UTM or in a rheometer.
Finally, as discussed in the next section, stress–strain curves can be predicted in computer simulations and/or phenomenological models, which, however, have mostly pertained (at least to date) to a simple shear geometry. In simulations, shear is mostly implemented using the Lees–Edwards periodic boundary condition, i.e., without boundary walls, and therefore focus on bulk properties. Extending numerical approaches to encompass compression, including boundaries, and thus model the case of a UTM, would certainly prove very useful not only for practical applications, but also for a deeper understanding of the influence of boundaries (see also Section 1.3 below).
Fig. 3 Stress–strain curves obtained from numerical computations of (a) a mean-field elasto-plastic model for increasing degree of annealing from bottom to top; the inset shows the amplitude of the stress drop vs. the degree of annealing (extracted from ref. 31), (b) a fluidity model with homogeneous flow enforced (dashed lines) and shear banding allowed (solid lines) for increasing waiting time from left to right (extracted from ref. 30), and (c) a fluidity model with decreasing initial fluidity from left to right (extracted from ref. 42). |
Still, the AQS limit can be difficult to achieve experimentally for colloidal systems. As a matter of fact, any small but finite shear rate (or the internal Brownian dynamics of the system) introduces a time scale that may affect the DBT. This time scale may not trigger the DBT,34 but it becomes important if one is interested in studying the statistical properties of the stick-slip-like dynamics that occur after the yield point. This leads to the key physical question of determining the correct theoretical framework to describe the DBT: should the DBT be described based on critical-like phenomena in the framework of statistical mechanics of disordered systems, or based on dynamical approaches as proposed in continuous models? Then, one should focus on the possible predictions, if any, which can be obtained using the different tools, and on the material properties that one should consider relevant in order to test the various approaches.
In ref. 31, a unifying picture was proposed to understand ductile vs. brittle behavior, where the crucial information is the preparation of the system. In this framework, the fictive temperature (i.e., the temperature at which the system has fallen out of equilibrium during cooling) becomes a relevant external parameter: at low fictive temperature (corresponding to very deep, very slowly cooled glassy states), the system exhibits a brittle response, whereas, at large fictive temperature (corresponding to poorly annealed systems), the solid-to-liquid transition is ductile. Similar results can be achieved for athermal systems using multiple control parameters (e.g., packing fraction, system size, annealing conditions, quenched disorder, etc.) to tune the material response.26,44 Here again, comparing possible predictions and available experimental data will be critical to draw definitive conclusions.
Finally, besides sample preparation, it remains unclear to know, at finite shear rate, what are the key differences in the DBT observed in Brownian vs. non Brownian suspensions, hard vs. soft particles, and jammed vs. glassy systems, as well as the interplay between these key parameters.45–47
In molecular dynamics simulations of a soft amorphous material subjected to AQS deformation,52 every rearrangement may only occur through the crossing of an instability threshold. As the external strain drives every local packing towards some instability threshold, the material is progressively brought into a state comprising a significant population of marginally stable regions.56 In such a system, Eshelby fields may trigger secondary events and thus give rise to plastic avalanches.3,4,57–61 These avalanches appear to be quite robust: moderate thermal fluctuations may trigger Eshelby events before they actually reach their instability threshold without changing significantly the avalanche dynamics,62 while finite strain rates introduce a cutoff on avalanche sizes.7,54 In model supercooled liquids, this avalanche regime is relevant to a very broad domain of temperatures and strain rates, and only ceases to be relevant in the Newtonian regime.62 In that context, a non-Newtonian behavior emerges from correlations between rearrangements that, like avalanches, result from the induction of secondary events via the Eshelby mechanism.
Within this picture, ductile behavior takes place when Eshelby events and plastic avalanches are broadly scattered throughout the material, while brittleness results from the localization of the plastic activity along shear bands. Avalanches appear as an intermediate between isolated events and shear bands, which can therefore be viewed as resulting from the emergence of spatial organization among avalanches,63 and not just of isolated events. The connection between Eshelby events, avalanches, and shear banding, however, remains a matter of debate and may involve the combination of elastic and structural effects.
Thus, an important challenge in the research on plasticity consists in identifying local predictors for plastic instabilities.64 In particular, it has become possible to characterize the proximity of mechanical instabilities via measurements of local yield stresses.65,66 These studies show that, in better relaxed systems, local packings present statistically larger values of local yield stresses. By identifying how the distribution of local barriers is affected by strain,38 they have also evidenced a significant role of rejuvenation: in a well-relaxed system, the first events to occur cause a local softening of the surrounding matrix, which facilitates the occurrence of more plastic events at the softened locations, which may lead to shear localization. From this perspective, hence, the DTB transition is governed by a dynamical interplay between plastic activity (including possible avalanches) and rejuvenation. It is favored by a high softness contrast between the initial state and rejuvenated packings, which explains that it is increasingly likely to take place in more relaxed (stable) system.
The identification of these effects leaves open the question of the role of elastic interactions between Eshelby quadrupoles in the formation of shear bands. This issue was partially clarified in previous studies,67,68 which have demonstrated that, with increasing shear strain γ, there is a threshold value γ* at which the minimum of the energy functional promotes an alignment of the quadrupoles along a line in two dimensions and along a plane in three dimensions. This alignment organizes the displacement field around a shear band. However, this theory had its limitations: the intensity of the quadrupoles and their core sizes were measured from simulations rather than arising directly from the theory itself.
Recent progress in understanding the role of quadrupolar plastic events in determining the mechanical response of amorphous solids69,70 presents an opportunity for achieving a self-consistent theory of shear banding. A key observation suggests that, in general, the plastic quadrupoles and their contribution to the energy functional renormalize the elastic moduli, yet they do not change the nature of linear elasticity theory. Conversely, when the quadrupolar field is non-uniform, gradients of this field act as dipoles, resulting in screening effects, the emergence of characteristic scales, and a fundamental change in mechanical response compared to the predictions of elasticity theory. These results were demonstrated in both simulations and experiments.69–74 Expanding this theory to address instabilities and shear banding necessitates further development, specifically including nonlinear quadrupolar and dipolar interactions. This ongoing work holds the potential to allow a self-consistent determination of the intensity and density of quadrupoles that are involved in creating the shear band. An example of such a theory of mechanical instability in metamaterials is outlined in ref. 75, and could be adapted for describing shear banding in amorphous solids.
Yet, a significant technical hurdle for future experiments is the detection of failure precursors. Although predictions of failure have been identified in the rheological signatures of soft materials in specific cases,76,77 rheometry typically is not sensitive enough to detect such precursors due to the small fraction of particles involved in microstructural rearrangements, which neither significantly nor reliably impact the macroscopic mechanical response.78,79 Detection should rely on multiprobe approaches that combine mechanical testing with methods that couple directly to precursors, among which X-ray and light scattering,29,55,80 electrical spectroscopy,81 acoustic emissions,82 change in the speed of sound,83 and fluorogenic mechanophores,84 to name a few, have been shown to unravel changes in the sample microstructure long before shear-induced or stress-induced failure.
However, these tools might not be sufficient. Indeed, despite the wealth of data collected by the last-generation array of seismic sensors along the San Andreas fault, no precursory signal was detected before the earthquake that took place in Parkfield.85 This episode illustrates that new strategies are required to detect precursors to failure accurately. A first strategy could be to design a setup combining different probes to monitor various length scales and time scales simultaneously. A more ambitious, second strategy entails designing innovative and more selective probes referred to as “smart probes,” capable, for instance, of detecting only non-affine motion. Such experimental probes could connect to recent theory and simulation advances that have identified microstructural features, like regions of low local yield stress or high “softness” parameter, that correlate with sites of initial failure.66,86 These technical developments will be pivotal in making significant progress in failure prediction. Finally, it appears evident that both structural heterogeneities87,88 as well the rigidity percolation transition89,90 play crucial roles in influencing the DBT in soft amorphous systems. Therefore, future experiments and simulations should examine more carefully how rigidity percolation and stress propagation, especially under well-controlled structural disorder, determine the fate of materials.
More generally, in soft amorphous materials, the DBT transition has been mainly considered when the system is forced at a constant shear rate in a simple, plane shear geometry. As noted above, it may be useful to understand the DBT in other geometries, e.g., under compression in the UTM. Moreover, it is important to understand whether or not the DBT also occurs under a constant external stress and/or in cyclic shear protocol. In order to compare the DBT in different geometries and under various testing protocols, both in experiments and in numerical computations, it is essential to introduce a quantitative measure of brittleness. Preliminary attempts focused on the energy release during brittle fracture and on stress drops in metallic glasses.91 Defining and investigating similar observables in soft glassy systems should help characterize the DBT and identify universal features, if any, common to soft and hard materials.
Finally, from a more general perspective, a key question remains about whether understanding the DBT provides any deeper insight into the shear-induced solid-to-liquid transition in amorphous systems, and into the physics of the yielding transition. Indeed, the solid-to-liquid transition is often considered as a dynamical phase transition as in other disordered systems.2,92 In the case of brittle materials, this point of view is consistent with the existence of an abrupt drop of the macroscopic stress at the yield point. In ductile materials, if the yielding transition occurs because of the nucleation of a shear band, the macroscopic stress does not show any abrupt change in time. However, the velocity gradient may exhibit a rather sharp interface in space, corresponding to a sharp change in the local shear rate. In the latter case, yielding still falls in the framework of a first-order (non-equilibrium) dynamical phase transition. Thus, one may argue that the DBT and the shear-induced solid-to-liquid transition could be discussed in a more general theoretical framework. This point deserves further investigation in the near future.
In ref. 102, criteria for the formation of shear bands were developed within a set of minimal assumptions, aimed at being independent of the particular constitutive model or material under consideration (and therefore applying not only to yield stress fluids, but complex fluids more widely). These criteria suggest that one might generically expect shear bands to form both (i) in shear startup, associated with the presence of a stress overshoot, and (ii) following creep under an imposed shear stress, when the shear rate signal curves upwards as the material yields. Both criteria are supported by observations in several experimental systems, particle based simulations and constitutive models, as reviewed in ref. 103.
However, the startup criterion (i, above) was derived within some caveats102–104 – most notably assuming only two interacting mechanical variables (the shear stress and one component of normal stress or one fluidity variable, for example) – and ref. 104 indeed confirmed it to hold in some models of polymeric fluids, but not others. It also showed that the magnitude of the banding effect depends on the solvent-to-solution viscosity ratio (in models of highly elastic polymeric fluids; the dependence is much weaker in models of yield stress fluids). Recent theoretical and numerical works105,106 confirmed that the start-up dynamics are influenced by variations in the viscoelastic constitutive model, and in the solvent-to-solution viscosity ratio;106 and that the correlation between the stress decay and the growth of linearized perturbations is not universal.105 The startup criterion has also been challenged in experiments on surfactant systems and in MCT simulations of repulsive glasses.93,107 However, wall slip plays a key role in these experiments, whereas the criterion was derived assuming zero slip.
Taken together, these works suggest that the evidence for shear banding associated with startup stress overshoot is quite widespread,103 but that further work is needed concerning its generality. The creep criterion (ii, above) holds for arbitrarily many interacting mechanical variables and is likely to be more generic.
A challenging question is to determine how many different theoretical frameworks are needed to describe the various types of shear bands and time scales associated with shear banding. For instance, in the absence of fracture, one could imagine building a consistent picture by generalizing known concepts from nucleation dynamics that occur in dynamical phase transitions (see, for instance, ref. 118). In all cases, a key quantity to predict theoretically is τf as a function of the external forcing mechanism and the material properties. The case of transient or long-lived shear bands is particularly interesting because modeling these systems is greatly simplified using a field theoretical approach, i.e., a continuous model, where, however, the choice of the order parameter is crucial. Then, in order to describe quantitatively the time evolution of the system, one needs to predict how the size of the shear band grows in time. By contrast, continuous models may not be required to describe the dynamics involved in the formation of fracture-like shear bands. Yet, in brittle materials, it is still unclear whether shear bands are the cause or the effect of the abrupt stress drop at the yielding transition. In the first case, shear band nucleation and dynamics should be described as for ductile materials,30 whereas in the second case, τf is dictated by a different mechanism, which may depend on the system state and preparation. For ductile-like materials, it was shown both in experiments and models that the nucleation of transient shear bands is the mechanism underlying the stress overshoot in start-up shear experiments.95,102,119,120 Finally, the impact of aging, possibly combined with inertia, on the yielding scenario, and on heterogeneous flows should also be considered.121–123
On more general grounds, any theoretical approach should be consistent with known robust features observed in most amorphous materials like, for instance, the aging properties of the system. This provides nontrivial constraints for possible theoretical frameworks. Moreover, such a theory should not only provide well-defined qualitative predictions to be checked but also be in quantitative agreement with existing experimental and numerical data. In particular, for ductile yielding, one should account for non-local rheological properties characterized by finite size effects and well-defined spatial scales. These features have been considered in recent developments, where the scaling properties of τf have been computed and successfully compared to experimental data.42,120,124 The presence of non-local effects naturally raises the problem of boundary conditions, which may play an important role, as discussed in more detail in Section 1.3 below.
Because experiments remain challenging, faster progress might come by comparing with numerical simulations. Recent molecular dynamics simulations have demonstrated that long-lasting heterogeneous flows and shear banding can be observed in both attractive particle systems127 and purely repulsive glasses.27,128 In the latter case, long-lasting transients are related to the presence of overconstrained microscopic domains trapped in the bulk of the sample upon preparation. Their structural reorganization under shear controls the emergence and the persistence of the shear banding. Finally, mesoscale simulations have shed some light on the delayed yielding observed during creep tests, i.e., constant stress experiments. A mean-field version of the elastoplastic scenario captured the main features of delayed yielding in soft glassy materials, including the nonmonotonic response of the shear rate, and the power-law scaling of the fluidization time, whose exponent increases for decreasing sample age.41
Fig. 4 (a) Schematic of the rheo-confocal setup for high-speed imaging, and a sketch of the shear cell, in which a second motor is positioned on the side to counter-rotate the lower glass plate. (b) Images taken during the shear-induced yielding of a depletion gel at 1 s−1. Projections in the velocity-gradient plane of shear reconstructed from fast axial scanning with a focus-tunable lens. The red arrow indicates the same particle, which drifted out of the selected orthogonal planes. Extracted from ref. 129. (c) Dimensionless master curve of the slip velocity V vs. stress at the wall σ as measured in microfluidic experiments performed on dense suspensions of microgels of various compositions. The parameters η and σE denote, respectively, the solvent viscosity and an elastic stress scale characterizing the nonlinear friction regime that is determined by fitting bulk velocity profiles. Different symbols stand for different Carbopol microgel samples. Inset: Sketch of a microgel of radius R in contact with the wall through a lubrication layer of thickness δ. Extracted from ref. 133. |
In more dilute dispersions, including gels, wall slip is also observed experimentally.99,145 To date, there is, however, no generic framework to describe wall slip in dilute systems, and it remains unclear whether its description can be coarse-grained to a mere comparison of well-defined length scales. Indeed, gels naturally form clusters due to interparticle attractive forces, which makes it hard to identify a relevant length scale between the particle size and the cluster size, which is itself a complex function of shear. Therefore, the multiscale architecture of colloidal gels, as well as the possibility for adhesion, rolling, and sliding between (aggregates of) colloids and the walls, constitute a significant challenge for models to offer a consistent description of wall slip in such dilute systems.
Finally, besides the steady state, wall slip plays a crucial role during transient flows, even in the presence of rough boundary conditions.96 For instance, in shear start-up flow, the sample may yield at the wall or in the bulk, which results in different dynamics, especially in systems that display strong aging dynamics.146 In that case, comparing two characteristic time scales, e.g., the aging time and the inverse of the plastic event rate near the wall, could be more relevant than comparing two length scales. Moreover, this picture prompts us to better understand the link between wall slip and plastic events in dilute systems, for these two concepts could be deeply related. In this context, a promising avenue lies in non-local effects, which offer an elegant way of coupling the material behavior near the wall to its bulk dynamics.120,124,147,148
Besides these fundamental considerations, addressing wall slip should help develop global strategies to suppress, reduce, and even tune wall slip for a broad variety of soft materials, from dense suspensions to colloidal gels, in steady-state flow as well as in transient regimes. One classical strategy used to reduce slip is to roughen the solid walls: for soft pastes, it has been shown empirically that a roughness comparable to the particle size efficiently reduces wall slip, whereas, for colloidal gels, the required roughness is not clearly identified, again due to the multiscale nature of such systems. A second strategy consists in modifying the interactions between the particles and the walls. For instance, it was shown that attractive wall–particle interactions may allow one to suppress the solvent lubrication layer and, consequently, wall slip for both soft pastes147 and colloidal gels.145,149 Such strategies will have a crucial impact in various engineering areas, including additive manufacturing.150–153
In control engineering, however, such systems, which relax over some characteristic time, would be considered as monostable, i.e., with no memory. Indeed, in a recent review158 of “memory formation in matter,” the word “memory” has been clearly defined as “the ability to encode, access, and erase signatures of past history in the state of a system.” A practical realization has been obtained by applying continuous oscillatory shear with a fixed strain amplitude on particulate systems, e.g., dense suspensions or bubble rafts. After a transient period where the sample microstructure is reorganized by the oscillatory strain, a steady state is reached, where oscillations do not produce any further rearrangement. A signature of the applied strain amplitude is thus “stored” by the system, and can be read out by observing at what strain the response of the system starts to be irreversible or display an abrupt change in the displacement of its constituents.159 In such a framework, a single system can store multiple memories by being submitting to more than one strain amplitude. These memories are transient and compete with each other, yet they can be stabilized if noise is added.160
Another strategy for encoding memory in soft glassy materials is based on bistability. There have been few reports to date of rheological bistability. One early example concerns lyotropic liquid crystalline phases of some aqueous surfactant solutions.161–163 There, stress first drives the solution from a low-viscosity lamellar phase to a high-viscosity phase made of an amorphous assembly of multilamellar spherulites, a.k.a., “disordered onions.” Depending on the surfactant concentration and temperature, these systems can further bifurcate to a low-viscosity phase with long-range order, either made of oriented lamellae or of a crystalline-like assembly of “ordered onions.” Such systems demonstrate “weak” bistability, in the sense that back-and-forth switching between the various states is not easy to control, as they involve complex, chaotic-like temporal dynamics and shear-banding flows.164–166 Similar examples of bistability can be found in surfactant hexagonal phases,167,168 and in the context of dense non-Brownian suspensions that exhibit discontinuous shear-thickening.169 A recent demonstration of “strong” bistability, where the system is shown to be switchable at will between two different structural states, concerns dispersions of non-Brownian repulsive particles in a background gel of much smaller, attractive colloids.170 By tuning the applied stress and strain accumulation, the system can cycle between solid and fluid states [see Fig. 5(a)–(d)]. Such a behavior allows the “writing” of memory through particular shear protocols, and the “reading” of the written state through, e.g., oscillatory rheology.
Fig. 5 Examples of memory effects in soft glassy materials. (a)–(d) Images of a fumed-silica-based binary suspension (small hydrophobic silica colloids ϕS = 0.02 and large silica microspheres ϕL = 0.3) in a glass vial, showing a liquid–solid transition under alternating high-shear vortex mixing and gentle roller mixing. Extracted from ref. 170. (e)–(h) Dependence of the stress response of a dense microgel suspension on the preshear direction. The shear rate = 0.25 s−1 is applied in the positive direction in all experiments, whereas the preshear (|0| = 20 s−1) is applied in the positive (black circles) or negative (blue open circles) direction. The sample is left at rest ( = 0) for a duration tw [= 1, 10, 1000, and 10000 s from (e) to (h)] between the preshear and the shear start-up test. Extracted from ref. 171. |
Soft amorphous systems with memory intrinsically show “history dependence.” Experimentally, this implies that specifying a well-defined protocol is vital. In this regard, one may usefully make contact with theoretical computer science,172,173 where the “algorithmic complexity” of a string is defined as the shortest algorithm that can be written to generate the string. Thus, ‘abababababab’ is generated by the pseudocode “[print ‘ab’ × 6],” while ‘aowgsiybmwsd’ is generated by “[print ‘a’, print ‘o’, …, print ‘d’],” and is more algorithmically complex. Along the same lines, one may perhaps define the concept of “rheological complexity,” which would encompass that of “rheological memory,” as the length of the shortest experimental protocol necessary for reproducing observations. Careful adherence to such protocol specifications is a prerequisite for enabling experiments to build on one another. It may also give some insight into the nature of the complexity in the systems under study.
Under continuous shear, yield stress fluids also display time-dependent features, often referred to as “memory effects,” which are closely related to thixotropy and/or viscoelasticity.178 So-called rheological hysteresis or hysteresis loops, which are obtained, e.g., by a decreasing ramp of shear rate, followed by the same increasing ramp, allow one to probe thixotropic time scales (for instance, by varying the sweep rate) that are characteristic of the material.179,180 Such time scales have been used to introduce a dimensionless number, the Mnemosyne number, defined as the product of the flow strength and the thixotropic time scale, and allows distinguishing thixotropic phenomena from other rheological responses such as viscoelasticity or plasticity.181 In that context, strong flows, i.e., large Mnemosyne numbers, are required to erase the memory from previous flows, e.g., inherited by loading the sample into the rheometer. Erasing memory has indeed long been studied in terms of defining good preshear protocols that allow for material behaviors to be studied irrespective of prior deformation.182
Last but not least, some memory can also be encoded via the direction of shear, as identified in a pioneering work on suspensions,183 and subsequently observed in cyclic shear experiments,160,184 and more recently quantified in shear start-up experiments performed along the same direction as a preshear, or the opposite [see Fig. 5(e)–(h)].171,182 In practice, a moderate shear rejuvenates a yield stress fluid only partially and may imprint some level of anisotropy in the sample microstructure, e.g., through some degree of alignment along the flow direction. Such microstructural anisotropy is then frozen upon flow cessation through the rapid liquid-to-solid transition, as we will discuss below in Section 2.2.1. Direct evidence for “directional memory”, also referred to as kinematic hardening or Bauschinger effect,185 has been reported in a broad range of systems, e.g., soft glasses,171 colloidal gels,99 or waxy crude oil.12 In that context, the sample memory is directly linked to the structural anisotropy of the yield stress material. The latter memory can then be read by orthogonal superposition rheometry.186 Finally, note that the sample memory has been also related to the recoverable strain that might be acquired prior to yielding, and which remains until deformation is applied in the opposite direction.
Finally, a key area of focus should be the experimental exploration of the local yielding properties. This entails examining the response at the microscopic or mesoscopic scales, following a loading phase at the macroscopic scale, and during the reverse experiment in which the shear is imposed in the opposite direction. Such investigations can be instrumental in identifying potential asymmetries in the local response, which may elucidate the mechanisms behind directional memory in soft glassy materials.188 Insights into strongly intermittent motion can be expected based on microrheology experiments performed in colloidal glass.189 A promising approach lies in leveraging exceptional experimental resolution to execute a systematic analysis akin to previous numerical studies.64,190
Fig. 6 Stress relaxation following flow cessation. (a) Monotonic stress relaxation vs. strain following various shear rates, reported here as a particulate Peclet numbers defined as the ratio of the strength of the shear flow with the interparticle bond force. The sample stores a non-negligible residual stress for Pe ≤ 10. Experiments performed on a gel of hard spheres of PMMA. Extracted from ref. 198. (b) Non-monotonic stress relaxation vs. time in a boehmite gel following a period of shear at constant shear rates before flow cessation. Each color corresponds to a different shear rate: from 3 s−1 (top) to 500 s−1 (bottom). Extracted from ref. 201. |
Future progress on residual stresses should come from local flow cessation experiments in which a spherical intruder is dragged through the sample of interest before being brought to a complete stop.203 Numerical simulations of active microrheology in a colloidal dispersion have already led to a micromechanical model of flow cessation, and have shown that the dispersion stress can be viewed as the duration of the microstructural memory.204,205 In general, additional experiments, particularly in gels, should be crucial in connecting the distorted microstructure inherited from flow cessation, the dynamics of local stress relaxation, and the macroscopic residual stress.
To conclude, one of the key open issues pertains to the microscopic scenario linked with the cessation of shear. While numerical results are available, there is, to date, no direct experimental evidence for structural anisotropy driving the dynamics of residual stresses. Future experiments, e.g., involving rheo-confocal or rheo-scattering experiments, are anticipated to provide an essential validation of the existing numerical scenarios and to prompt theoretical developments.
Fig. 7 (a) Phase diagram illustrating the jamming-glass transition within the framework of the self-propelled Voronoi model.225,235,236 This diagram explicates the conditions under which a confluent tissue transitions into a jammed state, maintaining a constant packing density. (b) Explores the scenario where the packing density within a tissue is variable. Under these circumstances, the jamming-glass behavior exhibits sensitivity to fluctuations in packing density. A phase diagram that encapsulates this behavior, drawing parallels to the jamming transitions observed in particulate systems, has been recently introduced.215 (a) and (b) are extracted from ref. 215 and 225, respectively. (c) Aspiration of an embryonic cellular aggregate through a constriction of size 50 μm used to generate large deformations to probe tissue elasto-capillary properties; picture extracted from ref. 237. (d) The shear response of a two-dimensional tissue using the vertex model shows the emergence of system-spanning tension chains. These structures bear a striking resemblance to the force chains observed in granular materials, suggesting a commonality in the stress response mechanisms between these disparate systems. (e) The flow curve resulting from a model tissue under shear, where the intricate balance between global external forces and local internal active drivers dictates the overall mechanical properties. In proximity to the jamming/unjamming transition, the tissue exhibits a spectrum of intriguing rheological behaviors. These include yielding, shear thinning, and both continuous (CST) and discontinuous shear thickening (DST). (d) and (e) are extracted from ref. 238. |
Furthermore, biological samples often actively respond to various stimuli, including mechanical forces.243–246 This poses a great challenge as the mechanical response of the system is also a function of the system's mechanosensing, i.e., there is feedback.247,248 This raises many open questions, among them: How do cell mechanosensing and active cell decision-making affect the emergent collective cell behavior and the mechanical properties on a tissue level? What is the detailed role of signaling pathways in such potential contributions? Accordingly, how can we computationally and experimentally study such an effect?
The vertex model is a sophisticated computational architecture designed to simulate and decode the mechanical attributes of tissue structures. It particularly focuses on the dynamics of cell packing, proliferation, and division.251–253 This model conceptualizes tissues as a two-dimensional mosaic of polygonal cells, with vertices representing cell junctions and edges denoting cell boundaries. By imbuing edges with tensional forces and cells with elastic properties, the model adeptly mimics the mechanical conduct of cells within a tissue matrix.
Vertex-based models have recently shed light on the phenomenon of glassiness in densely packed tissues.226,227 The self-propelled Voronoi model,225,235,236 in particular, has been noteworthy for its comprehensive integration of polarized cellular motion and the intricate web of cell–cell interactions within confluent tissues, where cells are tightly packed without gaps. This model, readily available as an open-source code,254 has been a cornerstone in making quantitative predictions about the glass transition phenomena in tissue biology [Fig. 7(a)]. The vertex model methodology, augmented by detail on biological feedback mechanisms, is becoming pivotal in dissecting the intricacies of morphogenesis – the developmental stages where tissues acquire their shape and structure – and in shedding light on the underlying physical principles governing these biological events.255–258
On the side of fundamentally understanding cell rheology, recent experiments on the aspiration of embryonic cell aggregates into a microfluidic constriction [see Fig. 7(c)] have revealed the existence of characteristic time scales and critical shear rates, fixing the domains in between which rearrangements can be observed.237 These experiments also showed that a viscosity could be difficult to define and measure in this context. All these different ingredients still remain to be incorporated into vertex models.
Moreover, these models are, for now, mainly used for describing cell monolayers that are considered two-dimensional systems. Nevertheless, it is not clear whether such a 2D assumption is always appropriate. Indeed, in monolayers, cells interact with a substrate and present an apico-basal polarity perpendicular to the substrate.246,259 More generally, including properly the effects of boundary conditions, such as the differential role of substrate friction and pair dissipation between cells, is crucial to describe the dynamics of deforming tissues. In experiments, the interactions with the substrate and the occurrence of slip can be adjusted through adhesive proteins, allowing cells to form biological focal adhesions.260 Furthermore, the role of boundary conditions is even more significant for developmental biology systems,261 as their entire behavior and response may depend on the rigidity of the substrate,262,263 the surrounding matrix in organoids or tumor spheroids, or on molecular properties,264 exhibiting mechanical or chemical sensitivity. This sensitivity mirrors the one of active systems, where the rigidity of the boundaries influences the dynamic response of active particles.265 This discovery raises intriguing questions about the interplay between boundary conditions and the dynamics of biological systems.
More generally, there are many open questions regarding the non-linear rheology of biological materials. For instance, do tissue models and biopolymer networks sit in the same stiffening universality class, and how does shear-stiffening influence ductility of tissues?238,268 What are the dominant/relevant length and time scales controlling the ductile–brittle transition in biological tissues? How does elastoplastic tissue flow differ from glassy elastoplastic flow, and how much information about such flow can be inferred from tissue cellular structure?269 How do cell division and apoptosis modify the mechanical noise compared to standard plasticity in passive systems, leading to non-linear rheological properties?270
The specific case of foodstuff is particularly relevant in terms of materials design, for edible materials are subject to intricate mechanical processes typically involving large deformations during production, mastication, and digestion. For example, pasta dough undergoes extrusion through a die to create long, thin strands, cheese curds are pressed to make cheese (with different textures, e.g., brittle, like the crumbly nature of feta cheese, or ductile, like the plastic flow of cream cheese), and bread dough is kneaded and folded to develop gluten structure and make a bread with a desirable texture. In the context of dairy products, the processing incorporates various principles covered in this workshop, i.e., shear-induced yielding, either brittle or ductile, as well as residual stresses and shear-induced texture frozen upon flow cessation. For instance, the preparation of stirred yogurt initially entails a sol–gel transition based on lactose fermentation that yields a brittle soft solid, which is then broken down into a dispersion of gel particles that dictates the texture and the sensory quality of the yogurt.278 Another striking example is related to the fabrication of butter. In practice, the cream is churned to separate the butterfat from the buttermilk. As a result of large deformations, the butterfat clumps together and forms butter, whose final texture is controlled by the shear history.279 The same conclusions hold for crystallization under shear in general, as evidenced in cocoa butter,280 and oil-in-water emulsions.281 These examples illustrate the relevance of studying dairy products as model systems282 to elucidate some of the open questions highlighted in the framework of this workshop.
(1) Role of shear bands at the yielding transition – experiments, numerical and theoretical results highlight the relevance of shear bands in understanding the yielding transition. For the brittle transition, one has to understand the cause-effect relation between shear band formation and the observed brittleness, and whether a possible universal picture is emerging. Somehow similar is the case for the ductile transition, where long-lasting shear bands are one possible mechanism to achieve complete fluidization. Last but not least, it remains to be understood in which conditions the latter scenario is replaced by a homogeneous fluidization that does not involve shear bands and yet shows a stress overshoot.93
(2) Quantitative measure of brittleness – although some proposals have emerged, there is a clear need to quantify brittleness in experiments and/or numerical simulations. This is a critical first step if one wants to compare results using different theoretical approaches or different material properties and/or preparations in a quantitative way.
(3) Theoretical frameworks – so far, there exist different theoretical frameworks to discuss and/or explain different material behavior at the yielding transition (e.g., DBT, creep flow, fluidization time, and rate dependence of the stress overshoot). Most of these frameworks or models show qualitative agreement with one or more features observed at the yielding transition. However, it is imperative to progress by offering quantitative predictions that can be compared against empirical data. Moreover, it is important to understand how different frameworks are linked with different material properties, identify possible limitations and generalizations, and, last but not least, suggest new experiments and/or data analysis to validate underlying theories.
(4) Microscopic view – several different topics need deeper investigations at the microscopic level, both experimentally and numerically. Here, we provide a non-exhaustive list:
• Identify microscopic signatures for the DBT, determine a microscopic scenario underpinning the existence of long-lasting shear bands, and investigate gels at the yield point.
• Investigate the role played by surface topology and microscopic interactions, if possible independently, on wall slip via numerical simulations and experiments.
• Identify particle-scale and/or mesoscopic dynamics relevant to memory effects and to the non-monotonic stress relaxation.
Moreover, a recent study283 indicates the intriguing possibility of reaching a state close to the yield point of an oxide glass by homogeneous irradiation rather than mechanical deformation. Comparing the macroscopic and microscopic properties of systems prepared in such different manners could lead to new insights concerning the yielding transition.
(5) Rheology of biological materials – there are many open questions about the interaction between single-cell behaviour and the large-scale rheology of living organisms. Most of the above points need to be investigated in detail for the specific case of biological materials. Another critical concern involves exploring potential connections between memory effects observed at the rheological level and the biochemical dynamics occurring at the single-cell level, if any such interaction exists.
(6) Materials design – use shear-induced memory effects to modify the properties of soft viscoelastic precursors employed to make hard materials, such as foodstuffs or cementitious materials.
This journal is © The Royal Society of Chemistry 2024 |