Modeling the Transition between Localized and Extended Deposition in Flow Networks through Packings of Glass Beads.

We use a theoretical model to explore how fluid dynamics, in particular, the pressure gradient and wall shear stress in a channel, affect the deposition of particles flowing in a microfluidic network. Experiments on transport of colloidal particles in pressure-driven systems of packed beads have shown that at lower pressure drop, particles deposit locally at the inlet, while at higher pressure drop, they deposit uniformly along the direction of flow. We develop a mathematical model and use agent-based simulations to capture these essential qualitative features observed in experiments. We explore the deposition profile over a two-dimensional phase diagram defined in terms of the pressure and shear stress threshold, and show that two distinct phases exist. We explain this apparent phase transition by drawing an analogy to simple one-dimensional mass-aggregation models in which the phase transition is calculated analytically.

Stress-stress correlations reveal force chains in gels.

We investigate the spatial correlations of microscopic stresses in soft particulate gels using 2D and 3D numerical simulations. We use a recently developed theoretical framework predicting the analytical form of stress-stress correlations in amorphous assemblies of athermal grains that acquire rigidity under an external load. These correlations exhibit a pinch-point singularity in Fourier space. This leads to long-range correlations and strong anisotropy in real space, which are at the origin of force-chains in granular solids. Our analysis of the model particulate gels at low particle volume fractions demonstrates that stress-stress correlations in these soft materials have characteristics very similar to those in granular solids and can be used to identify force chains. We show that the stress-stress correlations can distinguish floppy from rigid gel networks and that the intensity patterns reflect changes in shear moduli and network topology, due to the emergence of rigid structures during solidification.

Dilatancy, shear jamming, and a generalized jamming phase diagram of frictionless sphere packings.

Granular packings display the remarkable phenomenon of dilatancy, wherein their volume increases upon shear deformation. Conventional wisdom and previous results suggest that dilatancy, also being the related phenomenon of shear-induced jamming, requires frictional interactions. Here, we show that the occurrence of isotropic jamming densities φj above the minimal density (or the J-point density) φJ leads both to the emergence of shear-induced jamming and dilatancy in frictionless packings. Under constant pressure shear, the system evolves into a steady-state at sufficiently large strains, whose density only depends on the pressure and is insensitive to the initial jamming density φj. In the limit of vanishing pressure, the steady-state exhibits critical behavior at φJ. While packings with different φj values display equivalent scaling properties under compression, they exhibit striking differences in rheological behaviour under shear. The yield stress under constant volume shear increases discontinuously with density when φj > φJ, contrary to the continuous behaviour in generic packings that jam at φJ. Our results thus lead to a more coherent, generalised picture of jamming in frictionless packings, which also have important implications on how dilatancy is understood in the context of frictional granular matter.

Athermal Fluctuations in Disordered Crystals.

We analyze the fluctuations in particle positions and interparticle forces in disordered crystals composed of jammed soft particles in the limit of weak disorder. We demonstrate that such athermal systems are fundamentally different from their thermal counterparts, characterized by constrained fluctuations of forces perpendicular to the lattice directions. We develop a disorder perturbation expansion in polydispersity about the crystalline state, which we use to derive exact results to linear order. We show that constrained fluctuations result as a consequence of local force balance conditions, and are characterized by non-Gaussian distributions, which we derive exactly. We analytically predict several properties of such systems, including the scaling of the average coordination with polydispersity and packing fraction, which we verify with numerical simulations using soft disks with one-sided harmonic interactions.

Emergent Elasticity in Amorphous Solids.

The mechanical response of naturally abundant amorphous solids such as gels, jammed grains, and biological tissues are not described by the conventional paradigm of broken symmetry that defines crystalline elasticity. In contrast, the response of such athermal solids are governed by local conditions of mechanical equilibrium, i.e., force and torque balance of its constituents. Here we show that these constraints have the mathematical structure of a generalized electromagnetism, where the electrostatic limit successfully captures the anisotropic elasticity of amorphous solids. The emergence of elasticity from local mechanical constraints offers a new paradigm for systems with no broken symmetry, analogous to emergent gauge theories of quantum spin liquids. Specifically, our U(1) rank-2 symmetric tensor gauge theory of elasticity translates to the electromagnetism of fractonic phases of matter with the stress mapped to electric displacement and forces to vector charges. We corroborate our theoretical results with numerical simulations of soft frictionless disks in both two and three dimensions, and experiments on frictional disks in two dimensions. We also present experimental evidence indicating that force chains in granular media are subdimensional excitations of amorphous elasticity similar to fractons.

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.

Stress fluctuations in transient active networks.

Inspired by experiments on dynamic extensile gels of biofilaments and motors, we propose a model of a network of linear springs with kinetics consisting of growth at a prescribed rate, death after a lifetime drawn from a distribution, and birth at a randomly chosen node. The model captures features such as the build-up of self-stress, that are not easily incorporated into hydrodynamic theories. We study the model numerically and show that our observations can largely be understood through a stochastic effective-medium model. The resulting dynamically extending force-dipole network displays many features of yielded plastic solids, and offers a way to incorporate strongly non-affine effects into theories of active solids. A rather distinctive form for the stress distribution, and a Herschel-Bulkley dependence of stress on activity, are our major predictions.

Motion of active tracer in a lattice gas with cross-shaped particles.

We analyze the dynamics of an active tracer particle embedded in a thermal lattice gas. All particles are subject to exclusion up to third nearest neighbors on the square lattice, which leads to slow dynamics at high densities. For the case with no rotational diffusion of the tracer, we derive an analytical expression for the resulting drift velocity v of the tracer in terms of non-equilibrium density correlations involving the tracer particle and its neighbors, which we verify using numerical simulations. We show that the properties of the passive system alone do not adequately describe even this simple system of a single non-rotating active tracer. For large activity and low density, we develop an approximation for v. For the case where the tracer undergoes rotational diffusion independent of its neighbors, we relate its diffusion coefficient to the thermal diffusion coefficient and v. Finally, we study dynamics where the rotation of the tracer is limited by the presence of neighboring particles. We find that the effect of this rotational locking may be quantitatively described in terms of a reduction in the rotation rate.

Microscopic Origin of Frictional Rheology in Dense Suspensions: Correlations in Force Space.

We develop a statistical framework for the rheology of dense, non-Brownian suspensions, based on correlations in a space representing forces, which is dual to position space. Working with the ensemble of steady state configurations obtained from simulations of suspensions in two dimensions, we find that the anisotropy of the pair correlation function in force space changes with confining shear stress (σ_{xy}) and packing fraction (ϕ). Using these microscopic correlations, we build a statistical theory for the macroscopic friction coefficient: the anisotropy of the stress tensor, µ=σ_{xy}/P. We find that µ decreases (i) as ϕ is increased and (ii) as σ_{xy} is increased. Using a new constitutive relation between µ and viscosity for dense suspensions that generalizes the rate-independent one, we show that our theory predicts a discontinuous shear thickening flow diagram that is in good agreement with numerical simulations, and the qualitative features of µ that lead to the generic flow diagram of a discontinuous shear thickening fluid observed in experiments.

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

Jamming by shear.

A broad class of disordered materials including foams, glassy molecular systems, colloids and granular materials can form jammed states. A jammed system can resist small stresses without deforming irreversibly, whereas unjammed systems flow under any applied stresses. The broad applicability of the Liu-Nagel jamming concept has attracted intensive theoretical and modelling interest but has prompted less experimental effort. In the Liu-Nagel framework, jammed states of athermal systems exist only above a certain critical density. Although numerical simulations for particles that do not experience friction broadly support this idea, the nature of the jamming transition for frictional grains is less clear. Here we show that jamming of frictional, disk-shaped grains can be induced by the application of shear stress at densities lower than the critical value, at which isotropic (shear-free) jamming occurs. These jammed states have a much richer phenomenology than the isotropic jammed states: for small applied shear stresses, the states are fragile, with a strong force network that percolates only in one direction. A minimum shear stress is needed to create robust, shear-jammed states with a strong force network percolating in all directions. The transitions from unjammed to fragile states and from fragile to shear-jammed states are controlled by the fraction of force-bearing grains. The fractions at which these transitions occur are statistically independent of the density. Jammed states with densities lower than the critical value have an anisotropic fabric (contact network). The minimum anisotropy of shear-jammed states vanishes as the density approaches the critical value from below, in a manner reminiscent of an order-disorder transition.

Synchronization patterns in geometrically frustrated rings of relaxation oscillators.

Diffusively coupled chemical oscillators can exhibit a wide variety of complex spatial patterns. In this paper, we show that a ring of relaxation oscillators diffusively coupled through the inhibitory species leads to remarkable spatiotemporal patterns in the regime where there is a large separation of time scales between the activator and the inhibitor dynamics. The origin of these complex patterns can be traced back to a preponderance of antiphase synchronized states in the space of attractors. We provide an analytical explanation for the existence and stability of the antiphase synchronized states by examining the limit of extreme time scale separation. Numerical results on rings with small numbers of oscillators show that an explosion of patterns occurs for a ring with five oscillators.

Origin of rigidity in dry granular solids.

Solids are distinguished from fluids by their ability to resist shear. In traditional solids, the resistance to shear is associated with the emergence of broken translational symmetry as exhibited by a nonuniform density pattern. In this work, we focus on the emergence of shear rigidity in a class of solids where this paradigm is challenged. Dry granular materials have no energetically or entropically preferred density modulations. We show that, in contrast to traditional solids, the emergence of shear rigidity in these granular solids is a collective process, which is controlled solely by boundary forces, the constraints of force and torque balance, and the positivity of the contact forces. We develop a theoretical framework based on these constraints, which connects rigidity to broken translational symmetry in the space of forces, not positions of grains. We apply our theory to experimentally generated shear-jammed states and show that these states are indeed characterized by a persistent, non-uniform density modulation in force space, which emerges at the shear-jamming transition.

Segregation of polymers in confined spaces.

We investigate the motion of two overlapping polymers confined in a 2D box. A statistical model is constructed using blob-free-energy arguments. We find spontaneous segregation under the condition L > R([parallel]), and mixing under L < R([parallel]), where L is the length of the box and R([parallel]) is the polymer extension in an infinite slit. The segregation time (τ) is determined by solving a mean first-passage time problem and by performing Monte Carlo simulations. Both show a minimum in τ as a function of L. Although our results are restricted to 2D, the basic mechanism of competition between entropy and confinement leading to the minimum is suggestive of an evolutionary driving force for size selection.

Tensor electromagnetism and emergent elasticity in jammed solids.

The theory of mechanical response and stress transmission in disordered, jammed solids poses several open questions of how nonperiodic networks-apparently indistinguishable from a snapshot of a fluid-sustain shear. We present a stress-only theory of emergent elasticity for a nonthermal amorphous assembly of grains in a jammed solid, where each grain is subjected to mechanical constraints of force and torque balance. These grain-level constraints lead to the Gauss's law of an emergent U(1) tensor electromagnetism, which then accounts for the mechanical response of such solids. This formulation of amorphous elasticity has several immediate consequences. The mechanical response maps exactly to the static, dielectric response of this tensorial electromagnetism with the polarizability of the medium mapping to emergent elastic moduli. External forces act as vector electric charges, whereas the tensorial magnetic fields are sourced by momentum density. The dynamics in the electric and magnetic sectors naturally translate into the dynamics of the rigid jammed network and ballistic particle motion, respectively. The theoretical predictions for both stress-stress correlations and responses are borne out by the results of numerical simulations of frictionless granular packings in the static limit of the theory in both 2D and 3D.

Single polymer confinement in a tube: correlation between structure and dynamics.

In this paper, we construct an effective model for the dynamics of an excluded-volume chain under confinement by extending the formalism of Rouse modes. We make specific predictions about the behavior of the modes for a single polymer confined to a tube. The results are tested against Monte Carlo simulations using the bond-fluctuation algorithm which uses a lattice representation of the polymer chain with excluded-volume effects.

Statistical mechanics framework for static granular matter.

The physical properties of granular materials have been extensively studied in recent years. So far, however, there exists no theoretical framework which can explain the observations in a unified manner beyond the phenomenological jamming diagram. This work focuses on the case of static granular matter, where we have constructed a statistical ensemble which mirrors equilibrium statistical mechanics. This ensemble, which is based on the conservation properties of the stress tensor, is distinct from the original Edwards ensemble and applies to packings of deformable grains. We combine it with a field theoretical analysis of the packings, where the field is the Airy stress function derived from the force and torque balance conditions. In this framework, Point J characterized by a diverging stiffness of the pressure fluctuations. Separately, we present a phenomenological mean-field theory of the jamming transition, which incorporates the mean contact number as a variable. We link both approaches in the context of the marginal rigidity picture proposed by Wyart and others.

Dynamics in the metabasin space of a Lennard-Jones glass former: connectivity and transition rates.

Using simulations, we construct the effective dynamics in metabasin space for a Lennard-Jones glass former. Metabasins are identified via a scheme that measures transition rates between inherent structures and generates clusters of inherent structures by drawing in branches that have the largest transition rates. This construction is fundamentally different from the stochastic approach based on molecular-dynamics trajectories. The effective dynamics in this metabasin space is shown to be Markovian but to differ significantly from the simplest trap models. We specifically show that retaining information about the connectivity in this metabasin space is crucial for reproducing the slow dynamics observed in this system.

Growing length scale in gravity-driven dense granular flow.

We report simulations of a two-dimensional, dense, bidisperse system of inelastic hard disks falling down a vertical tube under the influence of gravity. We examine the approach to jamming as the average flow of particles down the tube is slowed by making the outlet narrower. Defining coarse-grained velocity and stress fields, we study two-point temporal and spatial correlation functions of these fields in a region of the tube where the time-averaged velocity is spatially uniform. We find that fluctuations in both velocity and stress become increasingly correlated as the system approaches jamming. We extract a growing length scale and time scale from these correlations.

Shapes of semiflexible polymers in confined spaces.

Biological macromolecules, living in the confines of a cell, often adopt conformations that are unlikely to occur in free space. In this paper, we investigate the effects of confinement on the shape of a semiflexible chain. Results of Monte Carlo simulations show the existence of a shape transition when the persistence length of the polymer becomes comparable to the dimensions of the box. An order parameter is introduced to quantify this behavior. A simple model is constructed to study the effect of the shape transition on the effective persistence length of the polymer.