Network materials, foams, and emulsions are ubiquitous in our daily life. We have a good intuition about how they respond as we handle them, but our theoretical understanding is poor. One of their most interesting features is that they are unusually fragile and appear to switch between solid and liquid state seamlessly. In fact, foams and emulsions undergo a non-equilibrium phase transition as their packing fraction increases - this is the jamming transition. Networks show a similar transition as their connectivity increases, where the material switches from sloppy to rigid. The fact that these materials undergo a phase transition, opens up the theoretical toolset of statistical mechanics. An important part of current research is therefore dedicated to finding diverging length and time scales and investigating the critical behavior of the systems in detail. Because the systems in question are highly disordered, analytical modeling is challenging. At the same time there are significant experimental obstacles to approaching the critical point closely. For this reason, the development of simulation software plays an important role - all data presented in this thesis is generated through simulations. As the subtitle of this dissertation suggests, our findings concern length, strain, and time scales which can be found in the linear response to external forces.@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

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

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

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

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

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

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