Unified study of glass and jamming rheology in soft particle systems

Nonlocal Effects in Inhomogeneous Flows of Soft Athermal Disks

We numerically investigate nonlocal effects on inhomogeneous flows of soft athermal disks close to but below their jamming transition. We employ molecular dynamics to simulate Kolmogorov flows, in which a sinusoidal flow profile with fixed wave number is externally imposed, resulting in a spatially inhomogeneous shear rate. We find that the resulting rheology is strongly wave-number-dependent, and that particle migration, while present, is not sufficient to describe the resulting stress profiles within a conventional local model. We show that, instead, stress profiles can be captured with nonlocal constitutive relations that account for gradients to fourth order. Unlike nonlocal flow in yield stress fluids, we find no evidence of a diverging length scale.

The thermal jamming transition of soft harmonic disks in two dimensions

By exploring the properties of the energy landscape of a bidisperse system of soft harmonic disks in two dimensions we determine the thermal jamming transition. To be specific, we study whether the ground state of the system where the particles do not overlap can be reached within a reasonable time. Starting with random initial configurations, the energy landscape is probed by energy minimization steps as in case of athermal jamming and in addition steps where an energy barrier can be crossed with a small but non-zero probability. For random initial conditions we find that as a function of packing fraction the thermal jamming transition, i.e. the transition from a state where all overlaps can be removed to an effectively non-ergodic state where one cannot get rid of the overlaps, occurs at a packing fraction of [Formula: see text], which is smaller than the transition packing fraction of athermal jamming at [Formula: see text]. Furthermore, we show that the thermal jamming transition is in the universality class of directed percolation and therefore is fundamentally different from the athermal jamming transition.

Marginally stable phases in mean-field structural glasses

A novel form of amorphous matter characterized by marginal stability was recently discovered in the mean-field theory of structural glasses. Using this approach, we provide complete phase diagrams delimiting the location of the marginally stable glass phase for a large variety of pair interactions and physical conditions, extensively exploring physical regimes relevant to granular matter, foams, emulsions, hard and soft colloids, and molecular glasses. We find that all types of glasses may become marginally stable, but the extent of the marginally stable phase highly depends on the preparation protocol. Our results suggest that marginal phases should be observable for colloidal and non-Brownian particles near jamming and for poorly annealed glasses. For well-annealed glasses, two distinct marginal phases are predicted. Our study unifies previous results on marginal stability in mean-field models and will be useful to guide numerical simulations and experiments aimed at detecting marginal stability in finite-dimensional amorphous materials.

The physics of jamming for granular materials: a review

Granular materials consist of macroscopic grains, interacting via contact forces, and unaffected by thermal fluctuations. They are one of a class systems that undergo jamming, i.e. a transition between fluid-like and disordered solid-like states. Roughly twenty years ago, proposals by Cates et al for the shear response of colloidal systems and by Liu and Nagel, for a universal jamming diagram in a parameter space of packing fraction, ϕ, shear stress, τ, and temperature, T raised key questions. Contemporaneously, experiments by Howell et al and numerical simulations by Radjai et al and by Luding et al helped provide a starting point to explore key insights into jamming for dry, cohesionless, granular materials. A recent experimental observation by Bi et al is that frictional granular materials have a a re-entrant region in their jamming diagram. In a range of ϕ, applying shear strain, γ, from an initially force/stress free state leads to fragile (in the sense of Cates et al), then anisotropic shear jammed states. Shear jamming at fixed ϕ is presumably conjugate to Reynolds dilatancy, involving dilation under shear against deformable boundaries. Numerical studies by Radjai and Roux showed that Reynolds dilatancy does not occur for frictionless systems. Recent numerical studies by several groups show that shear jamming occurs for finite, but not infinite, systems of frictionless grains. Shear jamming does not lead to known ordering in position space, but Sarkar et al showed that ordering occurs in a space of force tiles. Experimental studies seeking to understand random loose and random close packings (rlp and rcp) and dating back to Bernal have probed granular packings and their response to shear and intruder motion. These studies suggest that rlp's are anisotropic and shear-jammed-like, whereas rcp's are likely isotropically jammed states. Jammed states are inherently static, but the jamming diagram may provide a context for understanding rheology, i.e. dynamic shear in a variety of systems that include granular materials and suspensions.

A stability-reversibility map unifies elasticity, plasticity, yielding, and jamming in hard sphere glasses

Amorphous solids, such as glasses, have complex responses to deformations, with substantial consequences in material design and applications. In this respect, two intertwined aspects are important: stability and reversibility. It is crucial to understand, on the one hand, how a glass may become unstable due to increased plasticity under shear deformations, and, on the other hand, to what extent the response is reversible, meaning how much a system is able to recover the original configuration once the perturbation is released. Here, we focus on assemblies of hard spheres as the simplest model of amorphous solids such as colloidal glasses and granular matter. We prepare glass states quenched from equilibrium supercooled liquid states, which are obtained by using the swap Monte Carlo algorithm and correspond to a wide range of structural relaxation time scales. We exhaustively map out their stability and reversibility under volume and shear strains using extensive numerical simulations. The region on the volume-shear strain phase diagram where the original glass state remains solid is bounded by the shear yielding and the shear jamming lines that meet at a yielding-jamming crossover point. This solid phase can be further divided into two subphases: the stable glass phase, where the system deforms purely elastically and is totally reversible, and the marginal glass phase, where it experiences stochastic plastic deformations at mesoscopic scales and is partially irreversible. The details of the stability-reversibility map depend strongly on the quality of annealing of the glass. This study provides a unified framework for understanding elasticity, plasticity, yielding, and jamming in amorphous solids.

Crossover between Athermal Jamming and the Thermal Glass Transition of Suspensions

The non-Newtonian flow behavior of thermal and athermal disordered systems of dispersed uniform particles at high densities have strikingly similar features. By investigating the flow curves of yield-stress fluids and colloidal glasses having different volume fractions, particle sizes, and interactions, we show that both thermal and athermal systems exhibit power-law scaling with respect to the glass and jamming point, respectively, with the same exponents. All yield-stress flow curves can be scaled onto a single universal curve using the Laplace pressure as the stress scale for athermal systems and the osmotic pressure for the thermal systems. Strikingly, the details of interparticle interactions do not matter for the rescaling, showing that they are akin to usual phase transitions of the same universality class. The rescaling allows us to predict the flow properties of these systems from the volume fraction and known material properties.

Sticky Matters: Jamming and Rigid Cluster Statistics with Attractive Particle Interactions

While the large majority of theoretical and numerical studies of the jamming transition consider athermal packings of purely repulsive spheres, real complex fluids and soft solids generically display attraction between particles. By studying the statistics of rigid clusters in simulations of soft particles with an attractive shell, we present evidence for two distinct jamming scenarios. Strongly attractive systems undergo a continuous transition in which rigid clusters grow and ultimately diverge in size at a critical packing fraction. Purely repulsive and weakly attractive systems jam via a first-order transition, with no growing cluster size. We further show that the weakly attractive scenario is a finite size effect, so that for any nonzero attraction strength, a sufficiently large system will fall in the strongly attractive universality class. We therefore expect attractive jamming to be generic in the laboratory and in nature.

Discontinuous shear thickening in Brownian suspensions

Discontinuous shear thickening in dense suspensions naturally emerges from the activation of frictional forces by shear flow in non-Brownian systems close to jamming. Yet, this physical picture is incomplete as most experiments study soft colloidal particles subject to thermal fluctuations. To characterize discontinuous shear thickening in colloidal suspensions, we use computer simulations to provide a complete description of the competition between athermal jamming, frictional forces, thermal motion, particle softness, and shear flow. We intentionally neglect hydrodynamics, electrostatics, lubrication, and inertia, but can nevertheless achieve quantitative agreement with experimental findings. In particular, shear thickening corresponds to a crossover between frictionless and frictional jamming regimes which is controlled by thermal fluctuations and particle softness and occurs at a softness dependent Péclet number. We also explore the consequences of our findings for constant pressure experiments, and critically discuss the reported emergence of "S-shaped" flow curves. Our work provides the minimal ingredients to quantitatively interpret a large body of experimental work on discontinuous shear thickening in colloidal suspensions.

Glass transition of soft colloids

We explore the glassy dynamics of soft colloids using microgels and charged particles interacting by steric and screened Coulomb interactions, respectively. In the supercooled regime, the structural relaxation time τ_{α} of both systems grows steeply with volume fraction, reminiscent of the behavior of colloidal hard spheres. Computer simulations confirm that the growth of τ_{α} on approaching the glass transition is independent of particle softness. By contrast, softness becomes relevant at very large packing fractions when the system falls out of equilibrium. In this nonequilibrium regime, τ_{α} depends surprisingly weakly on packing fraction, and time correlation functions exhibit a compressed exponential decay consistent with stress-driven relaxation. The transition to this novel regime coincides with the onset of an anomalous decrease in local order with increasing density typical of ultrasoft systems. We propose that these peculiar dynamics results from the combination of the nonequilibrium aging dynamics expected in the glassy state and the tendency of colloids interacting through soft potentials to refluidize at high packing fractions.

Different routes into the glass state for soft thermo-sensitive colloids

We report an experimental and theoretical investigation of glass formation in soft thermo-sensitive colloids following two different routes: a gradual increase of the particle number density at constant temperature and an increase of the radius in a fixed volume at constant particle number density. Confocal microscopy experiments and the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory consistently show that the two routes lead to a dynamically comparable state at sufficiently long aging times. However, experiments reveal the presence of moderate but persistent structural differences. Successive cycles of radius decrease and increase lead instead to a reproducible glass state, indicating a suitable route to obtain rejuvenation without using shear fields.

No unjamming transition in a Voronoi model of biological tissue

Vertex models are a popular approach to simulating the mechanical and dynamical properties of dense biological tissues, describing the tissue as a network of polygons. Recently a class of two-dimensional vertex models was shown to exhibit a disordered rigidity transition controlled by the preferred cellular geometry, which was subsequently echoed by experimental findings. An attractive variant of these models uses a Voronoi tessellation to describe the cells, which reduces the number of degrees of freedom as compared the original vertex model. The Voronoi model was also endowed with a non-equilibrium model of cellular motility, leading to rich, glassy behavior. This glassy behavior was suggested to be inextricably linked to an underlying jamming transition. We test this conjecture, exploring the low-effective-temperature limit of the 2D Voronoi model by studying cell trajectories from detailed dynamical simulations in combination with rigidity measurements of energy-minimized disordered cell configurations. We find that the zero-temperature limit of this model has no unjamming transition. We show that this absence of an unjamming transition is intimately linked to the marginality of the model, i.e. the fact that the constraints imposed on cell areas and perimeters precisely balance the number of degrees of freedom in the model. Our work suggests that constraint counting arguments are useful to understand rigidity in a broad class of models of dense biological tissues.

Criticality of the zero-temperature jamming transition probed by self-propelled particles

We perform simulations of athermal systems of self-propelled particles (SPPs) interacting via harmonic repulsion in the vicinity of the zero-temperature jamming transition at point J. Every particle is propelled by a constant force f pointing to a randomly assigned and fixed direction. When f is smaller than the yield force f_y, the system is statically jammed. We find that f_y increases with packing fraction and exhibits finite size scaling, implying the criticality of point J. When f > f_y, SPPs flow forever and their velocities satisfy the k-Gamma distribution. Velocity distributions at various packing fractions and f collapse when the particle velocity is scaled by the average velocity v[combining macron], suggesting that v[combining macron] is a reasonable quantity to characterize the response to f. We thus define a response function R(ϕ,f) = v[combining macron](ϕ,f)/f. The function exhibits critical scaling nicely, implying again the criticality of point J. Our analysis and results indicate that systems of SPPs behave analogically to sheared systems, although their driving mechanisms are apparently distinct.

Ergodicity breaking transition in a glassy soft sphere system at small but non-zero temperatures

While the glass transition at non-zero temperature seems to be hard to access for experimental, theoretical, or simulation studies, jamming at zero temperature has been studied in great detail. Motivated by the exploration of the energy landscape that has been successfully used to investigate athermal jamming, we introduce a new method that includes the possibility of the thermally excited crossing of energy barriers. We then determine whether the ground state configurations of a soft sphere system are accessible or not and as a consequence whether the system is ergodic or effectively non-ergodic. Interestingly, we find an transition where the system becomes effectively non-ergodic if the density is increased. The transition density in the limit of small but non-zero temperatures is independent of temperature and below the transition density of athermal jamming. This confirms recent computer simulation studies where athermal jamming occurs deep inside the glass phase. In addition, we show that the ergodicity breaking transition is in the universality class of directed percolation. Therefore, our approach not only makes the transition from an ergodic to an effectively non-ergodic systems easily accessible and helps to reveal its universality class but also shows that it is fundamentally different from athermal jamming.

Viscous forces and bulk viscoelasticity near jamming

When weakly jammed packings of soft, viscous, non-Brownian spheres are probed mechanically, they respond with a complex admixture of elastic and viscous effects. While many of these effects are understood for specific, approximate models of the particles' interactions, there are a number of proposed force laws in the literature, especially for viscous interactions. We numerically measure the complex shear modulus G* of jammed packings for various viscous force laws that damp relative velocities between pairs of contacting particles or between a particle and the continuous fluid phase. We find a surprising sensitive dependence of G* on the viscous force law: the system may or may not display dynamic critical scaling, and the exponents describing how G* scales with frequency can change. We show that this sensitivity is closely linked to manner in which viscous damping couples to floppy-like, non-affine motion, which is prominent near jamming.

The liquid-glass-jamming transition in disordered ionic nanoemulsions

In quenched disordered out-of-equilibrium many-body colloidal systems, there are important distinctions between the glass transition, which is related to the onset of nonergodicity and loss of low-frequency relaxations caused by crowding, and the jamming transition, which is related to the dramatic increase in elasticity of the system caused by the deformation of constituent objects. For softer repulsive interaction potentials, these two transitions become increasingly smeared together, so measuring a clear distinction between where the glass ends and where jamming begins becomes very difficult or even impossible. Here, we investigate droplet dynamics in concentrated silicone oil-in-water nanoemulsions using light scattering. For zero or low NaCl electrolyte concentrations, interfacial repulsions are soft and longer in range, this transition sets in at lower concentrations, and the glass and the jamming regimes are smeared. However, at higher electrolyte concentrations the interactions are stiffer, and the characteristics of the glass-jamming transition resemble more closely the situation of disordered elastic spheres having sharp interfaces, so the glass and jamming regimes can be distinguished more clearly.

Dense fluidized granular media in microgravity

Handling and transport of granular media are inevitably governed by the settling of particles. Settling into a dense state is one of the defining characteristics of granular media, among dissipation and absence of thermal agitation. Hence, settling complicates the adaptation of microscopic theories from atomic, molecular, or colloidal media to granular media. It is desirable to provide experiments in which selectively one of the granular characteristics is tuned to test suitable adaptation of a theory. Here we show that gas fluidization of granular media in microgravity is a suitable approach to achieve steady states closer to thermally agitated systems free of settling. We use diffusing-wave spectroscopy to compare the spatial homogeneity and the microscopic dynamics of gas-fluidized granular media on the ground and in drop tower flights with increasing packing densities up to full arrest. The gas fluidization on the ground leads to inhomogeneous states as known from fluidized beds, and partial arrest occurs at packing fractions lower than the full arrested packing. The granular medium in microgravity in contrast attains a homogeneous state with complete mobilization even close to full arrest. Fluidized granular media thus can be studied in microgravity with dynamics and packing fractions not achievable on the ground.

On the apparent yield stress in non-Brownian magnetorheological fluids

We use simulations to probe the flow properties of dense two-dimensional magnetorheological fluids. Prior results from both experiments and simulations report that the shear stress σ scales with strain rate [small gamma, Greek, dot above] as σ ∼ [small gamma, Greek, dot above]^1-Δ, with values of the exponent ranging between 2/3 < Δ ≤ 1. However it remains unclear what properties of the system select the value of Δ, and in particular under what conditions the system displays a yield stress (Δ = 1). To address these questions, we perform simulations of a minimalistic model system in which particles interact via long ranged magnetic dipole forces, finite ranged elastic repulsion, and viscous damping. We find a surprising dependence of the apparent exponent Δ on the form of the viscous force law. For experimentally relevant values of the volume fraction ϕ and the dimensionless Mason number Mn (which quantifies the competition between viscous and magnetic stresses), models using a Stokes-like drag force show Δ ≈ 0.75 and no apparent yield stress. When dissipation occurs at the contact, however, a clear yield stress plateau is evident in the steady state flow curves. In either case, increasing ϕ towards the jamming transition suffices to induce a yield stress. We relate these qualitatively distinct flow curves to clustering mechanisms at the particle scale. For Stokes-like drag, the system builds up anisotropic, chain-like clusters as Mn tends to zero (vanishing strain rate and/or high field strength). For contact damping, by contrast, there is a second clustering mechanism due to inelastic collisions.

Jamming and overpacking fuzzy microgels: Deformation, interpenetration, and compression

Tuning the solubility of fuzzy polymer microgels by external triggers, such as temperature or pH, provides a unique mechanism for controlling the porosity and size of colloidal particles on the nanoscale. As a consequence, these smart microgel particles are being considered for applications ranging from viscosity modifiers and sensing to drug delivery and as models for the glass and the jamming transition. Despite their widespread use, little is known about how these soft particles adapt their shape and size under strong mechanical compression. We use a combination of precise labeling protocols and two-color superresolution microscopy to unravel the behavior of tracer microgels inside densely packed soft solids. We find that interpenetration and shape deformation are dominant until, in the highly overpacked state, this mechanism saturates and the only remaining way to further densify the system is by isotropic compression.

Echoes of the Glass Transition in Athermal Soft Spheres

Recent theoretical advances have led to the creation of a unified phase diagram for the thermal glass and athermal jamming transitions. This diagram makes clear that, while related, the mode-coupling-or dynamic-glass transition is distinct from the jamming transition, occurring at a finite temperature and significantly lower density than the jamming transition. Nonetheless, we demonstrate a prejamming transition in athermal frictionless spheres which occurs at the same density as the mode-coupling transition and is marked by percolating clusters of locally rigid particles. At this density in both the thermal and athermal systems, individual motions of an extensive number of particles become constrained, such that only collective motion is possible. This transition, which is well below jamming, exactly matches the definition of collective behavior at the dynamical transition of glasses. Thus, we reveal that the genesis of rigidity in both thermal and athermal systems is governed by the same underlying topological transition in their shared configuration space.

Advances and challenges in the rheology of concentrated emulsions and nanoemulsions

We review advances that have been made in the rheology of concentrated emulsions and nanoemulsions, which can serve as model soft materials that have highly tunable viscoelastic properties at droplet volume fractions near and above the glass transition and jamming point. As revealed by experiments, simulations, and theoretical models, interfacial and positional structures of droplets can depend on the applied flow history and osmotic pressure that an emulsion has experienced, thereby influencing its key rheological properties such as viscoelastic moduli, yield stress and strain, and flow behavior. We emphasize studies of monodisperse droplets, since these have led to breakthroughs in the fundamental understanding of dispersed soft matter. This review also covers the rheological properties of attractive emulsions, which can exhibit a dominant elasticity even at droplet volume fractions far below maximal random jamming of hard spheres.

Discontinuous change of shear modulus for frictional jammed granular materials

The shear modulus of jammed frictional granular materials with harmonic repulsive interaction under an oscillatory shear is numerically investigated. It is confirmed that the storage modulus, the real part of the shear modulus, for frictional grains with sufficiently small strain amplitude γ_{0} discontinuously emerges at the jamming transition point. The storage modulus for small γ_{0} differs from that of frictionless grains even in the zero friction limit, whereas they are almost identical with each other for sufficiently large γ_{0}, where the transition becomes continuous. The stress-strain curve exhibits a hysteresis loop even for a small strain, which connects a linear region for sufficiently small strain to another linear region for larger strain. We propose a scaling law to interpolate between the states of small and large γ_{0}.

Oscillatory rheology of dense, athermal suspensions of nearly hard spheres below the jamming point

The viscosity of a dense suspension has contributions from hydrodynamics and particle interactions, both of which depend upon the flow-induced arrangement of particles into fragile structures. Here, we study the response of nearly hard sphere suspensions to oscillatory shear using simulations and experiments in the athermal, non-inertial limit. Three distinct regimes are observed as a function of the strain amplitude γ₀. For γ₀ < 10^-1, initially non-contacting particles remain separated and the suspension behaves similarly to a Newtonian fluid, with the shear stress proportional to the strain rate, and the normal stresses close to zero. For γ₀ > 10¹, the microstructure becomes well-established at the beginning of each shear cycle and the rheology is quasi-Newtonian: the shear stress varies with the rate, but flow-induced structures lead to non-zero normal stresses. At intermediate γ₀, particle-particle contacts break and reform across entire oscillatory cycles, and we probe a non-linear regime that reveals the fragility of the material. Guided by these features, we further show that oscillatory shear may serve as a diagnostic tool to isolate specific stress contributions in dense suspensions, and more generally in those materials whose rheology has contributions with both hydrodynamic and non-hydrodynamic origin.

Scaling Theory for the Frictionless Unjamming Transition

We develop a scaling theory of the unjamming transition of soft frictionless disks in two dimensions by defining local areas, which can be uniquely assigned to each contact. These serve to define local order parameters, whose distribution exhibits divergences as the unjamming transition is approached. We derive scaling forms for these divergences from a mean-field approach that treats the local areas as noninteracting entities, and demonstrate that these results agree remarkably well with numerical simulations. We find that the asymptotic behavior of the scaling functions arises from the geometrical structure of the packing while the overall scaling with the compression energy depends on the force law. We use the scaling forms of the distributions to determine the scaling of the total grain area A_{G} and the total number of contacts N_{C}.

Shear Yielding and Shear Jamming of Dense Hard Sphere Glasses

We investigate the response of dense hard sphere glasses to a shear strain in a wide range of pressures ranging from the glass transition to the infinite-pressure jamming point. The phase diagram in the density-shear strain plane is calculated analytically using the mean-field infinite-dimensional solution. We find that just above the glass transition, the glass generically yields at a finite shear strain. The yielding transition in the mean-field picture is a spinodal point in presence of disorder. At higher densities, instead, we find that the glass generically jams at a finite shear strain: the jamming transition prevents yielding. The shear yielding and shear jamming lines merge in a critical point, close to which the system yields at extremely large shear stress. Around this point, highly nontrivial yielding dynamics, characterized by system-spanning disordered fractures, is expected.

Introduction to the Colloidal Glass Transition

Colloids are suspensions of small solid particles in a liquid and exhibit glassy behavior when the particle concentration is high. In these samples, the particles are roughly analogous to individual molecules in a traditional glass. This model system has been used to study the glass transition since the 1980s. In this Viewpoint I summarize some of the intriguing behaviors of the glass transition in colloids and discuss open questions.

Shear-induced solidification of athermal systems with weak attraction

We find that unjammed packings of frictionless particles with rather weak attraction can always be driven into solidlike states by shear. The structure of shear-driven solids evolves continuously with packing fraction from gel-like to jamminglike, but is almost independent of the shear stress. In contrast, both the density of vibrational states (DOVS) and force network evolve progressively with the shear stress. There exists a packing fraction independent shear stress σ_{c}, at which the shear-driven solids are isostatic and have a flattened DOVS. Solidlike states induced by a shear stress greater than σ_{c} possess properties of marginally jammed solids and are thus strictly defined shear jammed states. Below σ_{c}, shear-driven solids with rather different structures are all under isostaticity and share common features in the DOVS and force network. Our study leads to a jamming phase diagram for weakly attractive particles, which reveals the significance of the shear stress in determining properties of shear-driven solids and the connections and distinctions between jamminglike and gel-like states.

Effect of particle collisions in dense suspension flows

We study nonlocal effects associated with particle collisions in dense suspension flows, in the context of the Affine Solvent Model, known to capture various aspects of the jamming transition. We show that an individual collision changes significantly the velocity field on a characteristic volume Ω_{c}∼1/δz that diverges as jamming is approached, where δz is the deficit in coordination number required to jam the system. Such an event also affects the contact forces between particles on that same volume Ω_{c}, but this change is modest in relative terms, of order f_{coll}∼f[over ¯]^{0.8}, where f[over ¯] is the typical contact force scale. We then show that the requirement that coordination is stationary (such that a collision has a finite probability to open one contact elsewhere in the system) yields the scaling of the viscosity (or equivalently the viscous number) with coordination deficit δz. The same scaling result was derived [E. DeGiuli, G. Düring, E. Lerner, and M. Wyart, Phys. Rev. E 91, 062206 (2015)PLEEE81539-375510.1103/PhysRevE.91.062206] via different arguments making an additional assumption. The present approach gives a mechanistic justification as to why the correct finite size scaling volume behaves as 1/δz and can be used to recover a marginality condition known to characterize the distributions of contact forces and gaps in jammed packings.

Slow and long-ranged dynamical heterogeneities in dissipative fluids

A two-dimensional bidisperse granular fluid is shown to exhibit pronounced long-ranged dynamical heterogeneities as dynamical arrest is approached. Here we focus on the most direct approach to study these heterogeneities: we identify clusters of slow particles and determine their size, Nc, and their radius of gyration, RG. We show that , providing direct evidence that the most immobile particles arrange in fractal objects with a fractal dimension, df, that is observed to increase with packing fraction ϕ. The cluster size distribution obeys scaling, approaching an algebraic decay in the limit of structural arrest, i.e., ϕ→ϕc. Alternatively, dynamical heterogeneities are analyzed via the four-point structure factor S4(q,t) and the dynamical susceptibility χ4(t). S4(q,t) is shown to obey scaling in the full range of packing fractions, 0.6 ≤ϕ≤ 0.805, and to become increasingly long-ranged as ϕ→ϕc. Finite size scaling of χ4(t) provides a consistency check for the previously analyzed divergences of χ4(t) ∝ (ϕ-ϕc)(-γχ) and the correlation length ξ∝ (ϕ-ϕc)(-γξ). We check the robustness of our results with respect to our definition of mobility. The divergences and the scaling for ϕ→ϕc suggest a non-equilibrium glass transition which seems qualitatively independent of the coefficient of restitution.

Equilibrium Sampling of Hard Spheres up to the Jamming Density and Beyond

We implement and optimize a particle-swap Monte Carlo algorithm that allows us to thermalize a polydisperse system of hard spheres up to unprecedentedly large volume fractions, where previous algorithms and experiments fail to equilibrate. We show that no glass singularity intervenes before the jamming density, which we independently determine through two distinct nonequilibrium protocols. We demonstrate that equilibrium fluid and nonequilibrium jammed states can have the same density, showing that the jamming transition cannot be the end point of the fluid branch.

Clustering and flow around a sphere moving into a grain cloud

A bidimensional simulation of a sphere moving at constant velocity into a cloud of smaller spherical grains far from any boundaries and without gravity is presented with a non-smooth contact dynamics method. A dense granular "cluster" zone builds progressively around the moving sphere until a stationary regime appears with a constant upstream cluster size. The key point is that the upstream cluster size increases with the initial solid fraction [Formula: see text] but the cluster packing fraction takes an about constant value independent of [Formula: see text]. Although the upstream cluster size around the moving sphere diverges when [Formula: see text] approaches a critical value, the drag force exerted by the grains on the sphere does not. The detailed analysis of the local strain rate and local stress fields made in the non-parallel granular flow inside the cluster allows us to extract the local invariants of the two tensors: dilation rate, shear rate, pressure and shear stress. Despite different spatial variations of these invariants, the local friction coefficient μ appears to depend only on the local inertial number I as well as the local solid fraction, which means that a local rheology does exist in the present non-parallel flow. The key point is that the spatial variations of I inside the cluster do not depend on the sphere velocity and explore only a small range around the value one.

The glass and jamming transitions of soft polyelectrolyte microgel suspensions

We explore the influence of particle softness on the state diagram of well characterized polyelectrolyte microgel suspensions using dynamic light scattering and rheology. Upon increasing the polymer concentration, we cross successively the well defined glass and jamming transitions which delimit four different states: dilute colloidal suspension, entropic glass, jammed glass, and dense glass. Each state has a specific dynamical fingerprint dictated by two key ingredients related to particle softness: elastic contact interactions, and osmotic or steric deswelling. Soft interactions control yielding and flow of the jammed glasses. The shrinkage of the microgels makes the glass transition look smoother than in hard sphere suspensions. We quantify the relationship between the polymer concentration and the volume fraction, and show that the glass transition behaviour of soft microgels can be mapped to that of hard sphere glasses once the volume fraction is used as the control parameter.

Motility-driven glass and jamming transitions in biological tissues

Cell motion inside dense tissues governs many biological processes, including embryonic development and cancer metastasis, and recent experiments suggest that these tissues exhibit collective glassy behavior. To make quantitative predictions about glass transitions in tissues, we study a self-propelled Voronoi (SPV) model that simultaneously captures polarized cell motility and multi-body cell-cell interactions in a confluent tissue, where there are no gaps between cells. We demonstrate that the model exhibits a jamming transition from a solid-like state to a fluid-like state that is controlled by three parameters: the single-cell motile speed, the persistence time of single-cell tracks, and a target shape index that characterizes the competition between cell-cell adhesion and cortical tension. In contrast to traditional particulate glasses, we are able to identify an experimentally accessible structural order parameter that specifies the entire jamming surface as a function of model parameters. We demonstrate that a continuum Soft Glassy Rheology model precisely captures this transition in the limit of small persistence times, and explain how it fails in the limit of large persistence times. These results provide a framework for understanding the collective solid-to-liquid transitions that have been observed in embryonic development and cancer progression, which may be associated with Epithelial-to-Mesenchymal transition in these tissues.

Grains unchained: local fluidization of a granular packing by focused ultrasound

We report experimental results on the dynamics of a granular packing submitted to high-intensity focused ultrasound. Acoustic radiation pressure is shown to remotely induce local rearrangements within a pile as well as global motion around the focal spot in an initially jammed system. We demonstrate that this fluidization process is intermittent for a range of acoustic pressures and hysteretic when the pressure is cycled. Such a first-order-like unjamming transition is reproduced in numerical simulations in which the acoustic pressure field is modeled by a localized external force. Further analysis of the simulated packings suggests that in the intermittent regime unjamming is not associated with any noticeable prior structural signature. A simple two-state model based on effective temperatures is proposed to account for these findings.

Finite-temperature mechanical instability in disordered lattices

Mechanical instability takes different forms in various ordered and disordered systems and little is known about how thermal fluctuations affect different classes of mechanical instabilities. We develop an analytic theory involving renormalization of rigidity and coherent potential approximation that can be used to understand finite-temperature mechanical stabilities in various disordered systems. We use this theory to study two disordered lattices: a randomly diluted triangular lattice and a randomly braced square lattice. These two lattices belong to two different universality classes as they approach mechanical instability at T=0. We show that thermal fluctuations stabilize both lattices. In particular, the triangular lattice displays a critical regime in which the shear modulus scales as G∼T(1/2), whereas the square lattice shows G∼T(2/3). We discuss generic scaling laws for finite-T mechanical instabilities and relate them to experimental systems.

Shear thickening regimes of dense non-Brownian suspensions

We propose a unifying rheological framework for dense suspensions of non-Brownian spheres, predicting the onsets of particle friction and particle inertia as distinct shear thickening mechanisms, while capturing quasistatic and soft particle rheology at high volume fractions and shear rates respectively. Discrete element method simulations that take suitable account of hydrodynamic and particle-contact interactions corroborate the model predictions, demonstrating both mechanisms of shear thickening, and showing that they can occur concurrently with carefully selected particle surface properties under certain flow conditions. Microstructural transitions associated with frictional shear thickening are presented. We find very distinctive divergences of both microstructural and dynamic variables with respect to volume fraction in the thickened and non-thickened states.

Two-scale evolution during shear reversal in dense suspensions

We use shear-reversal simulations to explore the rheology of dense, non-Brownian, noninertial, suspensions, resolving lubrication forces between neighboring particles and modeling particle surface contacts. The transient stress response to an abrupt reversal of the direction of shear shows rate-independent, nonmonotonic behavior, capturing the salient features of the corresponding classical experiments. Based on analyses of the hydrodynamic and particle contact stresses and related contact networks, we demonstrate distinct responses at small and large strains, associated with contact breakage and structural reorientation, respectively, emphasizing the importance of particle contacts. Consequently, the hydrodynamic and contact stresses evolve over disparate strain scales and with opposite trends, resulting in nonmonotonic behavior when combined. We further elucidate the roles of particle roughness and repulsion in determining the microstructure and hence the stress response at each scale.

Out-of-Contact Elastohydrodynamic Deformation due to Lubrication Forces

We characterize the spatiotemporal deformation of an elastic film during the radial drainage of fluid from a narrowing gap. Elastic deformation of the film takes the form of a dimple and prevents full contact to be reached. With a thinner elastic film the stress becomes increasingly supported by the underlying rigid substrate and the dimple formation is suppressed, which allows the surfaces to reach full contact. We highlight the lag due to viscoelasticity on the surface profiles, and that for a given fluid film thickness deformation leads to stronger hydrodynamic forces than for rigid surfaces.

Soft matter food physics--the physics of food and cooking

This review discusses the (soft matter) physics of food. Although food is generally not considered as a typical model system for fundamental (soft matter) physics, a number of basic principles can be found in the interplay between the basic components of foods, water, oil/fat, proteins and carbohydrates. The review starts with the introduction and behavior of food-relevant molecules and discusses food-relevant properties and applications from their fundamental (multiscale) behavior. Typical food aspects from 'hard matter systems', such as chocolates or crystalline fats, to 'soft matter' in emulsions, dough, pasta and meat are covered and can be explained on a molecular basis. An important conclusion is the point that the macroscopic properties and the perception are defined by the molecular interplay on all length and time scales.

Constant Stress and Pressure Rheology of Colloidal Suspensions

We study the constant stress and pressure rheology of dense hard-sphere colloidal suspensions using Brownian dynamics simulation. Expressing the flow behavior in terms of the friction coefficient-the ratio of shear to normal stress-reveals a shear arrest point from the collapse of the rheological data in the non-Brownian limit. The flow curves agree quantitatively (when scaled) with the experiments of Boyer et al. [Phys. Rev. Lett. 107, 188301 (2011)]. Near suspension arrest, both the shear and the incremental normal viscosities display a universal power law divergence. This work shows the important role of jamming on the arrest of colloidal suspensions and illustrates the care needed when conducting and analyzing experiments and simulations near the flow-arrest transition.

Structure and energy of liquid foams

We present an overview of recent advances in the understanding of foam structure and energy and their dependence on liquid volume fraction. We consider liquid foams in equilibrium for which the relevant energy is surface energy. Measurements of osmotic pressure can be used to determine this as a function of liquid fraction in good agreement with results from computer simulations. This approach is particularly useful in the description of foams with high liquid content, so-called wet foams. For such foams X-ray tomography proves to be an important technique in analysing order and disorder. Much of the discussion in this article is also relevant to bi-liquid foams, i.e. emulsions, and to solid foams, provided that the solidification preserves the structure of the initially liquid foam template.

The structural origin of the hard-sphere glass transition in granular packing

Glass transition is accompanied by a rapid growth of the structural relaxation time and a concomitant decrease of configurational entropy. It remains unclear whether the transition has a thermodynamic origin, and whether the dynamic arrest is associated with the growth of a certain static order. Using granular packing as a model hard-sphere glass, we show the glass transition as a thermodynamic phase transition with a ‘hidden' polytetrahedral order. This polytetrahedral order is spatially correlated with the slow dynamics. It is geometrically frustrated and has a peculiar fractal dimension. Additionally, as the packing fraction increases, its growth follows an entropy-driven nucleation process, similar to that of the random first-order transition theory. Our study essentially identifies a long-sought-after structural glass order in hard-sphere glasses.

Glass transition shows dramatic dynamic slowdown, but its origin remains unclear. Here, Xia et al. observe in granular systems the rapid growth of a geometrically frustrated polytetrahedral order with packing fraction, which is spatially correlated with the slow dynamics.

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

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

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter.