Scholarship list
1–10 of 97 results
Journal article
Published 07/28/2026
Soft matter
Force transmission at large length scales is crucial for such biological functions as cell motility and morphogenesis. The networks that transmit these forces are malleable, patterned by active forces generated at the microscale by biological motors. In this paper we explore a simple model of a non-linear fiber network which has only two modes of deformation, but exhibits diverse mechanical phases with distinct large-scale response, tuned by the strength of a microscopic force dipole. We demonstrate, numerical simulations, that the network is remodeled by organized patterns of buckling. In addition, we use a mechanical screening theory of applied forces to extract an effective parameter that characterizes the spatial response of the stress field. This parameter shows strong variations that accompany changes in the buckling pattern. The emergent behavior at large length scales indicates that the medium screens out the details of the force dipole.
Preprint
Contact-network organization and motion statistics in shear-thickening suspensions
Posted to a preprint site 06/03/2026
We use lubricated-flow discrete-element-method (LF-DEM) simulations to examine how contact-network organization shapes particle motion in dense shear-thickening suspensions. The primary system studied is a two-dimensional bidisperse monolayer where rigid clusters are identified by the(3,3)pebble game; three-dimensional simulations are shown to have qualitatively similar rotational velocity statistics. Across the stress–solid-fraction state diagram, frictional contact number,k≥ 3percolation, and rigid-cluster fluctuations all strengthen in the same region where translational velocity correlations grow, consistent with rigid clusters translating coherently while the surrounding non-rigid particles accommodate a disproportionate share of the local velocity gradient. Rotational motion provides a complementary view: non-affine angular-velocity distributions broaden, near-contact rotations become increasingly anti-correlated, and rigid and non-rigid particles carry distinct statistics. Connectivity, rigidity, and velocity correlations are related but distinct signatures of the constrained collective motion that accompanies shear-thickening and the approach to shear jamming.
Journal article
Acoustic tuning of shear thickening suspensions: A universal scaling analysis
Published 03/2026
Journal of rheology (New York : 1978), 70, 2, 419 - 425
Tuning shear thickening behavior is a longstanding problem in the field of dense suspensions. Acoustic perturbations offer a convenient way to control shear thickening in real time, opening the door to a new class of smart materials. However, complete control over shear thickening requires a quantitative description for how suspension viscosity varies under acoustic perturbation. Here, we achieve this goal by experimentally probing suspensions with acoustic perturbations and incorporating their effect on the suspension viscosity into a universal scaling framework where the viscosity is described by a scaling function, which captures a crossover from the frictionless jamming critical point to a frictional shear jamming critical point. Our analysis reveals that the effect of acoustic perturbations may be explained by the introduction of an effective interparticle repulsion whose magnitude is roughly equal to the acoustic energy density. Furthermore, we show how this scaling framework may be leveraged to produce explicit predictions for the viscosity of a dense suspension under acoustic perturbation. Our results demonstrate the utility of the scaling framework for experimentally manipulating shear thickening systems.
Preprint
The Jammed Phase of Infinitely Persistent Active Matter
Published 02/24/2026
We study an extreme active matter system, which is essentially a dense assembly of athermal, soft and infinitely persistent active particles. Using extensive numerical simulations we obtain jammed configurations of this system in two dimensions and probe the stability of such structures under increasing active forcing magnitude. We show that the critical active forcing magnitude for the jammed phase to yield scales with virial pressure asf_(c)∼ pᵅ , withα=1.17 , describing the yielding line. Using a Laplacian framework, we redistribute the active forces into a modified contact force network. By analysing the statistics of these redistributed forces, we obtain a very robust scaling law consistent with the passive limit, not just near the unjamming line, but in the entire jammed active phase. The probability distribution of the magnitude of the contact force deviates from the power-law form found in passive systems for values smaller than the active force. Moreover, within the jammed phase, the system displays elastic, plastic, and yielding events with increasing active forcing. This active plasticity appears abruptly and can not be captured by the continuous softening of the Hessian spectrum. However, we demonstrate that the Hessian still retains the ability to predict relaxation times. These results clarify how activity modifies force distributions and leads to deformation, plasticity and yielding in dense, jammed, infinitely persistent active matter.
Journal article
Contact network structures and rigidity development in two-dimensional bidisperse suspensions
Published 01/23/2026
Journal of rheology, 70, 2, 343 - 359
Dense non-Brownian suspensions with conservative repulsive forces between the particles are known to exhibit shear thickening, where viscosity increases with applied stress due to a change in the dominant stress mechanism. At low stress, repulsion maintains liquid films that
lubricate particle interactions, while higher stress overcomes the repulsion to generate frictional contacts and leads to greater flow resistance.
Here, shear-thickened suspensions are studied in stress-controlled simulations incorporating hydrodynamic, electrostatic double-layer repulsion, and frictional forces; two-dimensional monolayers are studied for monodisperse and bidisperse suspensions with size ratios δ ¼ as=al
from 1.0 to 4.0, where as and al are the small and large particle radii. Small-particle fractions ζ ¼ fs=f ¼ 0:25, 0:50, and 0:75 are considered. Total area fractions of 0:71 f 0:82 are studied, with the larger values at greater size ratios. Flow curves for mono- and bidisperse
systems under varying stress are analyzed, along with detailed structural comparisons for different interparticle friction. We examine the
approach to shear jamming, through the emergence of rigid local clusters generated by the reduction of degrees of freedom by frictional contacts. The variance of the fraction of particles in rigid clusters increases sharply near the jamming solid fraction, consistent with a secondorder phase transition description of the phenomenon. The contact fabric tensor is determined to provide a measure of the structural anisotropy.
Preprint
Universal scaling of shear thickening suspensions under acoustic perturbation
Published 10/13/2025
arXiv (Cornell University)
Tuning shear thickening behavior is a longstanding problem in the field of dense suspensions. Acoustic perturbations offer a convenient way to control shear thickening in real time, opening the door to a new class of smart materials. However, complete control over shear thickening requires a quantitative description for how suspension viscosity varies under acoustic perturbation. Here, we achieve this goal by experimentally probing suspensions with acoustic perturbations and incorporating their effect on the suspension viscosity into a universal scaling framework where the viscosity is described by a scaling function, which captures a crossover from the frictionless jamming critical point to a frictional shear jamming critical point. Our analysis reveals that the effect of acoustic perturbations may be explained by the introduction of an effective interparticle repulsion whose magnitude is roughly equal to the acoustic energy density. Furthermore, we demonstrate how this scaling framework may be leveraged to produce explicit predictions for the viscosity of a dense suspension under acoustic perturbation. Our results demonstrate the utility of the scaling framework for experimentally manipulating shear thickening systems.
Preprint
Stress Response of Jammed Solids: Prestress and Screening
Published 09/17/2025
arXiv (Cornell University)
Unlike classical elasticity, where stresses arise from deformations relative to a stress-free reference configuration, rigidity in amorphous systems is maintained by disordered force networks that generate internal prestress. Previously, we introduced a ''stress-only'' formulation, where mechanical equilibrium resembles Gauss's law in a rank-2 tensor electrostatics with vector charges, and demonstrated that the mechanical response of jammed solids is described by the dielectric response of this gauge-theoretic formulation. Here, we extend this framework by incorporating scale-dependent screening that captures both dielectric and Debye-type behaviour. This introduces a characteristic length scale in stress correlations as well as in the response to external forces. Through numerical simulations of soft-sphere packings, we show that this length scale is set by the particle size, thus providing a natural ultraviolet cutoff while preserving long-wavelength emergent elasticity. We show that this lengthscale remains finite for all pressures, with no evidence for an emergent Debye-like screening near the frictionless unjamming transition. We demonstrate that although individual realisations show strong fluctuations, disorder averaging at fixed macroscopic conditions yields a robust dielectric-like response that persists up to unjamming. Finally, we also provide a physical interpretation of the gauge field within the electrostatic mapping: relative grain displacements in response to localised external perturbations correspond to difference in the gauge field, linking the field-theoretic description to particle-level mechanics.
Journal article
Pinch-point Singularities in Stress-Stress Correlations Reveal Rigidity in Colloidal Gels
Published 06/28/2025
Soft matter, 24
We demonstrate that the spatial correlations of microscopic stresses in 2D model colloidal gels obtained in computer simulations can be quantitatively
described by the predictions of a theory for emergent elasticity of pre-stressed solids (vector charge theory). By combining a rigidity analysis with the characterization provided by the stress correlations, we show that the theoretical predictions are able to distinguish rigid from floppy gels, and quantify that distinction in terms of the size of a pinch-point singularity emerging at large length scales, which, in the theory, directly derives from the constraints imposed by mechanical equilibrium on the internal forces. We also use the theoretical predictions to investigate the coupling between stress-transmission and rigidity, and we explore the possibility of a Debye-like screening mechanism that would modify the theory predictions below a characteristic length scale.
Journal article
Pinch-point singularities in stressstress correlations reveal rigidity in colloidal gels
Published 06/19/2025
Soft matter, 21, 24, 4812 - 4821
We demonstrate that the spatial correlations of microscopic stresses in 2D model colloidal gels obtained in computer simulations can be quantitatively described by the predictions of a theory for emergent elasticity of pre-stressed solids (vector charge theory). By combining a rigidity analysis with the characterization provided by the stress correlations, we show that the theoretical predictions are able to distinguish rigid from floppy gels, and quantify that distinction in terms of the size of a pinch-point singularity emerging at large length scales, which, in the theory, directly derives from the constraints imposed by mechanical equilibrium on the internal forces. We also use the theoretical predictions to investigate the coupling between stresstransmission and rigidity, and we explore the possibility of a Debye-like screening mechanism that would modify the theory predictions below a characteristic length scale. The loss of rigidity in particulate gels can be quantified by the disappearance of a characteristic pinch-point singularity in the gel stress-stress correlation function.
Journal article
Rigidity transitions in anisotropic networks: a crossover scaling analysis
Published 04/03/2025
Soft matter
We study how the rigidity transition in a triangular lattice changes as a function of anisotropy by preferentially filling bonds on the lattice in one direction. We discover that the onset of rigidity in anisotropic spring networks on a regular triangular lattice arises in at least two steps, reminiscent of the two-step melting transition in two dimensional crystals. In particular, our simulations demonstrate that the percolation of stress-supporting bonds happens at different critical volume fractions along different directions. By examining each independent component of the elasticity tensor, we determine universal exponents and develop universal scaling functions to analyze isotropic rigidity percolation as a multicritical point. Our crossover scaling approach is applicable to anisotropic biological materials ( cellular cytoskeletons, extracellular networks of tissues like tendons), and extensions to this analysis are important for the strain stiffening of these materials.