We show that foams and emulsions can display a fundamentally different normal response to a simple shear deformation. While foams dilate or push outwards on the shearing surfaces, known as a positive Poynting effect, in emulsions the Poynting effect can have either sign and can be tuned by changing the emulsion properties. We relate the sign of the Poynting effect to the presence of a compressible contact network supported by adhesive contacts. When the concentration of surfactant in the continuous phase is low, the emulsions are nonadhesive and push outward on their shearing surfaces, as do the foams. When the surfactant concentration is increased, the emulsions become adhesive due to depletion interactions, and the Poynting effect changes sign. We argue that the adhesive contact network develops a shear modulus that stiffens in response to dilation, which leads to the negative Poynting effect.

@en

We determine how low frequency vibrational modes control the elastic shear modulus of Mikado networks, a minimal mechanical model for semi-flexible fiber networks. From prior work it is known that when the fiber bending modulus is sufficiently small, (i) the shear modulus of 2D Mikado networks scales as a power law in the fiber line density, G ∼ ρα+1, and (ii) the networks also possess an anomalous abundance of soft (low-frequency) vibrational modes with a characteristic frequency ωκ ∼ ρβ/2. While it has been suggested that α and β are identical, the preponderance of evidence indicates that α is larger than theoretical predictions for β. We resolve this inconsistency by measuring the vibrational density of states in Mikado networks for the first time. Supported by these results, we then demonstrate analytically that α = β + 1. In so doing, we uncover new insights into the coupling between soft modes and shear, as well as the origin of the crossover from bending- to stretching-dominated response. This journal is

@en

Numerous soft materials jam into an amorphous solid at a high packing fraction. This nonequilibrium phase transition is best understood in a model system where particles repel when they overlap. Recently, however, it was shown that introducing any finite amount of attraction between particles changes the universality class of the transition. The properties of this "sticky jamming"class remain almost entirely unexplored. We use molecular dynamics simulations and scaling analysis to determine the shear modulus, bulk modulus, and coordination of marginal solids close to the sticky jamming point. Each observable differs not just quantitatively but also qualitatively from the purely repulsive case.

@en

We numerically investigate stress relaxation in soft athermal disks to reveal critical slowing down when the system approaches the jamming point. The exponents describing the divergence of the relaxation time differ dramatically depending on whether the transition is approached from the jammed or unjammed phase. This contrasts sharply with conventional dynamic critical scaling scenarios, where a single exponent characterizes both sides. We explain this surprising difference in terms of the vibrational density of states, which is a key ingredient of linear viscoelastic theory. The vibrational density of states exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.

@en

We investigate irreversibility in soft frictionless disk packings on approach to the unjamming transition. Using simulations of shear reversal tests, we study the relationship between plastic work and irreversible rearrangements of the contact network. Infinitesimal strains are reversible, while any finite strain generates plastic work and contact changes in a sufficiently large packing. The number of irreversible contact changes grows with strain, and the stress–strain curve displays a crossover from linear to increasingly nonlinear response when the fraction of irreversible contact changes approaches unity.

@en

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.

@en

We present an approach to understand geometric-incompatibility–induced rigidity in underconstrained materials, including subisostatic 2D spring networks and 2D and 3D vertex models for dense biological tissues. We show that in all these models a geometric criterion, represented by a minimal length ℓ¯min, determines the onset of prestresses and rigidity. This allows us to predict not only the correct scalings for the elastic material properties, but also the precise magnitudes for bulk modulus and shear modulus discontinuities at the rigidity transition as well as the magnitude of the Poynting effect. We also predict from first principles that the ratio of the excess shear modulus to the shear stress should be inversely proportional to the critical strain with a prefactor of 3. We propose that this factor of 3 is a general hallmark of geometrically induced rigidity in underconstrained materials and could be used to distinguish this effect from nonlinear mechanics of single components in experiments. Finally, our results may lay important foundations for ways to estimate ℓ¯min from measurements of local geometric structure and thus help develop methods to characterize large-scale mechanical properties from imaging data.
@en

Yield-stress materials form an interesting class of materials that behave like solids at small stresses, but start to flow once a critical stress is exceeded. It has already been reported both in experimental and simulation work that flow curves of different yield-stress materials can be scaled with the distance to jamming or with the confining pressure. However, different scaling exponents are found between experiments and simulations. In this paper we identify sources of this discrepancy. We numerically relate the volume fraction with the confining pressure and discuss the similarities and differences between rotational and oscillatory measurements. Whereas simulations are performed in the elastic response regime close to the jamming transition and with very small amplitudes to calculate the scaling exponents, these conditions are hardly possible to achieve experimentally. Measurements are often performed far away from the critical volume fraction and at large amplitudes. We show that these differences are the underlying reason for the different exponents for rescaling flow curves.@en

Aqueous foams are an important model system that displays coarsening dynamics. Coarsening in dispersions and foams is well understood in the dilute and dry limits, where the gas fraction tends to zero and one, respectively. However, foams are known to undergo a jamming transition from a fluidlike to a solidlike state at an intermediate gas fraction φc. Much less is known about coarsening dynamics in wet foams near jamming, and the link to mechanical response, if any, remains poorly understood. Here we probe coarsening and mechanical response using numerical simulations of a variant of the Durian bubble model for wet foams. As in other coarsening systems we find a steady state scaling regime with an associated particle size distribution. We relate the time rate of evolution of the coarsening process to the wetness of the foam and identify a characteristic coarsening time that diverges approaching jamming. We further probe mechanical response of the system to strain while undergoing coarsening. There are two competing timescales, namely the coarsening time and the mechanical relaxation time. We relate these to the evolution of the elastic response and the mechanical structure.

@en

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.

@en

When elastic solids are sheared, a nonlinear effect named after Poynting gives rise to normal stresses or changes in volume. We provide a novel relation between the Poynting effect and the microscopic Grüneisen parameter, which quantifies how stretching shifts vibrational modes. By applying this relation to random spring networks, a minimal model for, e.g., biopolymer gels and solid foams, we find that networks contract or develop tension because they vibrate faster when stretched. The amplitude of the Poynting effect is sensitive to the network's linear elastic moduli, which can be tuned via its preparation protocol and connectivity. Finally, we show that the Poynting effect can be used to predict the finite strain scale where the material stiffens under shear.@en

Many soft matter systems are confined in some but not all dimensions; examples include microfluidic channels and inclined plane flows. Hence it is important to characterize finite size effects not only as a function of volume, but also for varying aspect ratio. For soft sphere packings close to the jamming transition, finite size effects are well understood, but only in square and cubic systems. In these cases there is clear evidence for a critical volume that diverges at jamming, but it is equally clear that this picture must break down for extreme aspect ratios. We perform simulations of soft spheres near jamming in two and three dimensions for aspect ratios as large as 1024. In addition to the previously identified critical volume, we find evidence for a non-trivial length scale that diverges at the jamming point.@en

Solids deform and fluids flow, but soft glassy materials, such as emulsions, foams, suspensions, and pastes, exhibit an intricate mix of solid- and liquid-like behavior. While much progress has been made to understand their elastic (small strain) and flow (infinite strain) properties, such understanding is lacking for the softening and yielding phenomena that connect these asymptotic regimes. Here we present a comprehensive framework for softening and yielding of soft glassy materials, based on extensive numerical simulations of oscillatory rheological tests, and show that two distinct scenarios unfold depending on the material's packing density. For dense systems, there is a single, pressure-independent strain where the elastic modulus drops and the particle motion becomes diffusive. In contrast, for weakly jammed systems, a two-step process arises: at an intermediate softening strain, the elastic and loss moduli both drop down and then reach a new plateau value, whereas the particle motion becomes diffusive at the distinctly larger yield strain. We show that softening is associated with an extensive number of microscopic contact changes leading to a non-analytic rheological signature. Moreover, the scaling of the softening strain with pressure suggest the existence of a novel pressure scale above which softening and yielding coincide, and we verify the existence of this crossover scale numerically. Our findings thus evidence the existence of two distinct classes of soft glassy materials – jamming dominated and dense – and show how these can be distinguished by their rheological fingerprint.@en

Materials like foams and emulsions display complex rhelogical behavior close to their jamming transition. When driven too hard the initial linear stress-strain response breaks down and the material softens. Using simulations of soft repulsive spheres, we characterize the softening crossover by establishing the relevant strain scale below which linear response is valid. We further perform shear reversal tests to investigate the interplay between proximity to jamming and the onset of irreversibility.@en

We use simulations of frictionless soft sphere packings to identify novel constitutive relations for linear elasticity near the jamming transition. By forcing packings at varying wavelengths, we directly access their transverse and longitudinal compliances. These are found to be wavelength dependent, in violation of conventional (local) linear elasticity. Crossovers in the compliances select characteristic length scales, which signify the appearance of nonlocal effects. Two of these length scales diverge as the pressure vanishes, indicating that critical effects near jamming control the breakdown of local elasticity. We expect these nonlocal constitutive relations to be applicable to a wide range of weakly jammed solids, including emulsions, foams, and granulates.

@en

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.@en

We report the results of molecular dynamics simulations of stress relaxation tests in athermal viscous soft sphere packings close to their unjamming transition. By systematically and simultaneously varying both the amplitude of the applied strain step and the pressure of the initial condition, we access both linear and nonlinear response regimes and control the distance to jamming. Stress relaxation in viscoelastic solids is characterized by a relaxation time τ∗ that separates short time scales, where viscous loss is substantial, from long time scales, where elastic storage dominates and the response is essentially quasistatic. We identify two distinct plateaus in the strain dependence of the relaxation time, one each in the linear and nonlinear regimes. The height of both plateaus scales as an inverse power law with the distance to jamming. By probing the time evolution of particle velocities during relaxation, we further identify a correlation between mechanical relaxation in the bulk and the degree of non-affinity in the particle velocities on the micro scale.@en

While jamming is best understood in the context of purely repulsive soft spheres, emulsions and other experimental realizations of the soft sphere model commonly display weakly cohesive forces. We perform simulations of soft spheres with a finite-ranged attractive tail in the pair potential. The resulting attractive soft sphere packings can be stable at volume fractions below the purely repulsive jamming point. These new jammed states have counter-intuitive properties -- for example, while attraction introduces tensile forces, their presence leads to an increase in the compressive stress. We use critical scaling analysis to characterize the geometry and mechanics of attractive soft sphere packings as a function of both the volume fraction and the range of the attractive interaction.@en

In many soft matter systems with a network structure, the mean connectivity plays a key role in determining the mechanical response. In the typical scenario for central force networks, there is a critical connectivity (the isostatic point) above which the system can withstand shear and bulk deformations. However, there are several different ways to induce the material to rigidify below the critical connectivity. Here we explore three such mechanisms: bending rigidity in the Mikado model for semi-flexible fiber networks, pre-tensioning in random spring networks, and finite-ranged attraction in soft sphere packings. Here we show common features in the rigidity transition of all three seemingly disparate systems. In particular we identify a band of low frequency normal modes for low perturbation strengths whose height and width are related to the distance to isostaticity. These in turn control the elastic moduli, as we explain with simple scaling arguments.@en

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 as σ ∼ ^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.

@en