Jamming transition and normal modes of polydispersed soft particle packing.

The jamming transition of soft particles characterized by narrow size distributions has been well studied by physicists. However, polydispersed systems are more relevant to engineering, and the influence of polydispersity on jamming phenomena is still unexplored. Here, we numerically investigate jamming transitions of polydispersed soft particles in two dimensions.

Moduli and modes in the Mikado model.

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, ∼ , and (ii) the networks also possess an anomalous abundance of soft (low-frequency) vibrational modes with a characteristic frequency ∼ . While it has been suggested that and are identical, the preponderance of evidence indicates that is larger than theoretical predictions for .

Stress Relaxation above and below the Jamming Transition.

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.

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.

A minimal-length approach unifies rigidity in underconstrained materials.

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 [Formula: see text], 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.

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.

Coarsening and mechanics in the bubble model for wet foams.

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}.

Normal Stresses, Contraction, and Stiffening in Sheared Elastic Networks.

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.

Softening and yielding of soft glassy materials.

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.

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.

Stress relaxation in viscous soft spheres.

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.

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], 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).

Nonlocal Elasticity near Jamming in Frictionless Soft Spheres.

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.

Contact changes of sheared systems: Scaling, correlations, and mechanisms.

We probe the onset and effect of contact changes in two-dimensional soft harmonic particle packings which are sheared quasistatically under controlled strain. First, we show that, in the majority of cases, the first contact changes correspond to the creation or breaking of contacts on a single particle, with contact breaking overwhelmingly likely for low pressures and/or small systems, and contact making and breaking equally likely for large pressures and in the thermodynamic limit. The statistics of the corresponding strains are near-Poissonian, in particular for large-enough systems.

Beyond linear elasticity: jammed solids at finite shear strain and rate.

The shear response of soft solids can be modeled with linear elasticity, provided the forcing is slow and weak. Both of these approximations must break down when the material loses rigidity, such as in foams and emulsions at their (un)jamming point - suggesting that the window of linear elastic response near jamming is exceedingly narrow. Yet precisely when and how this breakdown occurs remains unclear.

Jamming in finite systems: stability, anisotropy, fluctuations, and scaling.

Athermal packings of soft repulsive spheres exhibit a sharp jamming transition in the thermodynamic limit. Upon further compression, various structural and mechanical properties display clean power-law behavior over many decades in pressure. As with any phase transition, the rounding of such behavior in finite systems close to the transition plays an important role in understanding the nature of the transition itself.

Contact changes near jamming.

We probe the onset and effect of contact changes in soft harmonic particle packings which are sheared quasistatically. We find that the first contact changes are the creation or breaking of contacts on a single particle. We characterize the critical strain, statistics of breaking versus making a contact, and ratio of shear modulus before and after such events, and explain their finite size scaling relations.

Dynamic critical response in damped random spring networks.

The isostatic state plays a central role in organizing the response of many amorphous materials. We construct a diverging length scale in nearly isostatic spring networks that is defined both above and below isostaticity and at finite frequencies and relate the length scale to viscoelastic response. Numerical measurements verify that proximity to isostaticity controls the viscosity, shear modulus, and creep of random networks.

Soft-sphere packings at finite pressure but unstable to shear.

When are athermal soft-sphere packings jammed? Any experimentally relevant definition must, at the very least, require a jammed packing to resist shear. We demonstrate that widely used (numerical) protocols, in which particles are compressed together, can and do produce packings that are unstable to shear-and that the probability of generating such packings reaches one near jamming. We introduce a new protocol which, by allowing the system to explore different box shapes as it equilibrates, generates truly jammed packings with strictly positive shear moduli G.

Relaxations and rheology near jamming.

We determine the form of the complex shear modulus G* in soft sphere packings near jamming. Viscoelastic response at finite frequency is closely tied to a packing's intrinsic relaxational modes, which are distinct from the vibrational modes of undamped packings. We demonstrate and explain the appearance of an anomalous excess of slowly relaxing modes near jamming, reflected in a diverging relaxational density of states.

Model for the scaling of stresses and fluctuations in flows near jamming.

We probe flows of soft, viscous spheres near the jamming point, which acts as a critical point for static soft spheres. Starting from energy considerations, we find nontrivial scaling of velocity fluctuations with strain rate. Combining this scaling with insights from jamming, we arrive at an analytical model that predicts four distinct regimes of flow, each characterized by rational-valued scaling exponents.

Entropy maximization in the force network ensemble for granular solids.

A long-standing issue in the area of granular media is the tail of the force distribution, in particular, whether this is exponential, Gaussian, or even some other form. Here we resolve the issue for the case of the force network ensemble in two dimensions. We demonstrate that conservation of the total area of a reciprocal tiling, a direct consequence of local force balance, is crucial for predicting the local stress distribution.

Nonlinear elastic stress response in granular packings.

We study the nonlinear elastic response of a two-dimensional material to a localized boundary force, with the particular goal of understanding the differences observed between isotropic granular materials and those with hexagonal anisotropy. Corrections to the classical Boussinesq result for the stresses in an infinite half space of a linear, isotropic material are developed in a power series in inverse distance from the point of application of the force. The breakdown of continuum theory on scales of order of the grain size is modeled with phenomenological parameters characterizing the strengths of induced multipoles near the point of application of the external force.

Force distributions in a triangular lattice of rigid bars.

We study the uniformly weighted ensemble of force balanced configurations on a triangular network of nontensile contact forces. For periodic boundary conditions corresponding to isotropic compressive stress, we find that the probability distribution for single-contact forces decays faster than exponentially. This superexponential decay persists in lattices diluted to the rigidity percolation threshold.