Long-range correlations in elastic moduli and local stresses at the unjamming transition.

We explore the behaviour of spatially heterogeneous elastic moduli as well as the correlations between local moduli in model solids with short-range repulsive potentials. We show through numerical simulations that local elastic moduli exhibit long-range correlations, similar to correlations in the local stresses. Specifically, the correlations in local shear moduli exhibit anisotropic behavior at large lengthscales characterized by pinch-point singularities in Fourier space, displaying a structural pattern akin to shear stress correlations. Focussing on two-dimensional jammed solids approaching the unjamming transition, we show that stress correlations exhibit universal properties, characterized by a quadratic p2 dependence of the correlations as the pressure p approaches zero, independent of the details of the model. In contrast, the modulus correlations exhibit a power-law dependence with different exponents depending on the specific interaction potential. Furthermore, we illustrate that while affine responses lack long-range correlations, the total modulus, which encompasses non-affine behavior, exhibits long-range correlations.

Universal stress correlations in crystalline and amorphous packings.

We present a universal characterization of stress correlations in athermal systems, across crystalline to amorphous packings. Via numerical analysis of static configurations of particles interacting through harmonic as well as Lennard-Jones potentials, for a variety of preparation protocols and ranges of microscopic disorder, we show that the properties of the stress correlations at large lengthscales are surprisingly universal across all situations, independent of structural correlations, or the correlations in orientational order. In the near-crystalline limit, we present exact results for the stress correlations for both models, which work surprisingly well at large lengthscales, even in the amorphous phase. Finally, we study the differences in stress fluctuations across the amorphization transition, where stress correlations reveal the loss of periodicity in the structure at short lengthscales with increasing disorder.

Generalized disorder averages and current fluctuations in run and tumble particles.

We present exact results for the fluctuations in the number of particles crossing the origin up to time t in a collection of noninteracting run and tumble particles in one dimension. In contrast to passive systems, such active particles are endowed with two inherent degrees of freedom, positions and velocities, which can be used to construct density and magnetization fields. We introduce generalized disorder averages associated with both these fields and perform annealed and quenched averages over various initial conditions. We show that the variance σ^{2} of the current in annealed versus quenched magnetization situations exhibits a surprising difference at short times, σ^{2}∼t vs σ^{2}∼t^{2}, respectively, with a sqrt[t] behavior emerging at large times. Our analytical results demonstrate that in the strictly quenched scenario, where both the density and magnetization fields are initially frozen, the fluctuations in the current are strongly suppressed. Importantly, these anomalous fluctuations cannot be obtained solely by freezing the density field.

Green's functions for random resistor networks.

We analyze random resistor networks through a study of lattice Green's functions in arbitrary dimensions. We develop a systematic disorder perturbation expansion to describe the weak disorder regime of such a system. We use this formulation to compute ensemble-averaged nodal voltages and bond currents in a hierarchical fashion. We verify the validity of this expansion with direct numerical simulations of a square lattice with resistances at each bond exponentially distributed. Additionally, we construct a formalism to recursively obtain the exact Green's functions for finitely many disordered bonds. We provide explicit expressions for lattices with up to four disordered bonds, which can be used to predict nodal voltage distributions for arbitrarily large disorder strengths. Finally, we introduce a novel order parameter that measures the overlap between the bond current and the optimal path (the path of least resistance), for a given resistance configuration, which helps to characterize the weak and strong disorder regimes of the system.

Enhanced vibrational stability in glass droplets.

We show through simulations of amorphous solids prepared in open-boundary conditions that they possess significantly fewer low-frequency vibrational modes compared to their periodic boundary counterparts. Specifically, using measurements of the vibrational density of states, we find that the D(ω)∼ω4 law changes to D(ω)∼ωδ with δ≈5 in two dimensions and δ≈4.5 in three dimensions. Crucially, this enhanced stability is achieved when utilizing slow annealing protocols to generate solid configurations. We perform an anharmonic analysis of the minima corresponding to the lowest frequency modes in such open-boundary systems and discuss their correlation with the density of states. A study of various system sizes further reveals that small systems display a higher degree of localization in vibrations. Lastly, we confine open-boundary solids in order to introduce macroscopic stresses in the system, which are absent in the unconfined system and find that the D(ω)∼ω4 behavior is recovered.

Emergent power-law interactions in near-crystalline membranes.

We derive exact results for the fluctuations in energy produced by microscopic disorder in near-crystalline athermal systems. Our formalism captures the heterogeneity in the elastic energy of polydisperse soft disks in energy-minimized configurations. We use this to predict the distribution of interaction energy between two defects in a disordered background. We show that this interaction energy displays a disorder-averaged power-law behavior 〈δE〉∼Δ^{-4} at large distances Δ between the defects. These interactions upon disorder average also display the sixfold symmetry of the underlying reference crystal. Additionally, we show that the fluctuations in the interaction energy encode the athermal correlations introduced by the disordered background. We verify our predictions with energy-minimized configurations of polydisperse soft disks in two dimensions.

Active random walks in one and two dimensions.

We investigate active lattice walks: biased continuous time random walks which perform orientational diffusion between lattice directions in one and two spatial dimensions. We study the occupation probability of an arbitrary site on the lattice in one and two dimensions and derive exact results in the continuum limit. Next, we compute the large deviation free-energy function in both one and two dimensions, which we use to compute the moments and the cumulants of the displacements exactly at late times. Our exact results demonstrate that the cross-correlations between the motion in the x and y directions in two dimensions persist in the large deviation function. We also demonstrate that the large deviation function of an active particle with diffusion displays two regimes, with differing diffusive behaviors. We verify our analytic results with kinetic Monte Carlo simulations of an active lattice walker in one and two dimensions.

Observation of two-step aggregation kinetics of amyloid-ßproteins from fractal analysis.

Self-aggregation in proteins has long been studied and modeled due to its ubiquity and importance in many biological contexts. Several models propose a two step aggregation mechanism, consisting of linear growth of fibrils and secondary growth involving branch formation. Single molecule imaging techniques such as total internal reflection fluorescence (TIRF) microscopy can provide direct evidence of such mechanisms, however, analyzing such large data-sets is challenging. In this paper, we analyze for the first time, images of growing amyloid fibrils obtained from TIRF microscopy using the techniques of fractal geometry, which provides a natural framework to disentangle the two types of growth mechanisms at play. We find that after an initial linear growth phase, identified by a plateau in the average fractal dimension with time, the occurrence of branching events leads to a further increase in the fractal dimension, with a final saturation value≈2. This provides direct evidence of the two-step nature of the aggregation kinetics of amyloid-ßproteins, with an initial linear elongation phase followed by branching at later times.

Universal non-Debye low-frequency vibrations in sheared amorphous solids.

We study energy minimised configurations of amorphous solids with a simple shear degree of freedom. We show that the low-frequency regime of the vibrational density of states of structural glass formers is crucially sensitive to the macroscopic stress of the sampled configurations. In both two and three dimensions, shear-stabilised configurations display a D(ωmin) ∼ ω5min regime, as opposed to the ω4min regime observed under unstrained conditions. In order to isolate the source of these deviations from crystalline behaviour, we also study configurations of two dimensional, strained amorphous solids close to a plastic event. We show that the minimum eigenvalue distribution at a strain 'γ' near the plastic event occurring at 'γP' assumes a universal form that displays a collapse when scaled by , and with the number of particles as N-0.22. Notably, at low frequencies, this scaled distribution displays a robust D(ωmin) ∼ ω6min power-law regime, which survives in the large N limit. Finally, we probe the properties of these configurations through a characterisation of the second and third eigenvalues of the Hessian matrix near a plastic event.

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.

First-contact-breaking distributions in strained disordered crystals.

We derive exact probability distributions for the strain (ε) at which the first stress drop event occurs in uniformly strained disordered crystals, with quenched disorder introduced through polydispersity in particle sizes. We characterize these first stress drop events numerically as well as theoretically and identify them with the first-contact-breaking event in the system. Our theoretical results are corroborated with numerical simulations of quasistatic volumetric strain applied to disordered near-crystalline configurations of athermal soft particles. We develop a general technique to determine the distribution of strains at which the first stress drop events occur, through an exact mapping between the cumulative distribution of first-contact-breaking events and the volume of a convex polytope whose dimension is determined by the number of defects N_{d} in the system. An exact numerical computation of this polytope volume for systems with small numbers of defects displays a remarkable match with the distribution of strains generated through direct numerical simulations. Finally, we derive the distribution of strains at which the first stress drop occurs, assuming that individual contact-breaking events are uncorrelated, which accurately reproduces distributions obtained from direct numerical simulations.

Disorder perturbation expansion for athermal crystals.

We introduce a perturbation expansion for athermal systems that allows an exact determination of displacement fields away from the crystalline state as a response to disorder. We show that the displacement fields in energy-minimized configurations of particles interacting through central potentials with microscopic disorder can be obtained as a series expansion in the strength of the disorder. We introduce a hierarchy of force-balance equations that allows an order-by-order determination of the displacement fields, with the solutions at lower orders providing sources for the higher-order solutions. This allows the simultaneous force-balance equations to be solved, within a hierarchical perturbation expansion to arbitrary accuracy. We present exact results for an isotropic defect introduced into the crystalline ground state at linear order and second order in our expansion. We show that the displacement fields produced by the defect display interesting self-similar properties at every order. We derive a |δr|∼1/r and |δf|∼1/r^{2} decay for the displacement fields and excess interparticle forces at large distances r away from the defect. Finally, we derive nonlinear corrections introduced by the interactions between defects at second order in our expansion. We verify our exact results with displacement fields obtained from energy-minimized configurations of soft disks.

Long-range correlations in pinned athermal networks.

We derive exact results for displacement fields that develop as a response to external pinning forces in two-dimensional athermal networks. For a triangular lattice arrangement of particles interacting through soft potentials, we develop a Green's function formalism which we use to derive exact results for displacement fields produced by localized external forces. We show that in the continuum limit the displacement fields decay as 1/r at large distances r away from a force dipole. Finally, we extend our formulation to study correlations in the displacement fields produced by the external pinning forces. We show that uncorrelated pinned forces at each vertex give rise to long-range correlations in displacements in athermal systems, with a nontrivial system size dependence. We verify our predictions with numerical simulations of athermal networks in two dimensions.

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.

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.

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

Gaps between avalanches in one-dimensional random-field Ising models.

We analyze the statistics of gaps (ΔH) between successive avalanches in one-dimensional random-field Ising models (RFIMs) in an external field H at zero temperature. In the first part of the paper we study the nearest-neighbor ferromagnetic RFIM. We map the sequence of avalanches in this system to a nonhomogeneous Poisson process with an H-dependent rate ρ(H). We use this to analytically compute the distribution of gaps P(ΔH) between avalanches as the field is increased monotonically from -∞ to +∞. We show that P(ΔH) tends to a constant C(R) as ΔH→0^{+}, which displays a nontrivial behavior with the strength of disorder R. We verify our predictions with numerical simulations. In the second part of the paper, motivated by avalanche gap distributions in driven disordered amorphous solids, we study a long-range antiferromagnetic RFIM. This model displays a gapped behavior P(ΔH)=0 up to a system size dependent offset value ΔH_{off}, and P(ΔH)∼(ΔH-ΔH_{off})^{θ} as ΔH→H_{off}^{+}. We perform numerical simulations on this model and determine θ≈0.95(5). We also discuss mechanisms which would lead to a nonzero exponent θ for general spin models with quenched random fields.

Columnar order and Ashkin-Teller criticality in mixtures of hard squares and dimers.

We show that critical exponents of the transition to columnar order in a mixture of 2×1 dimers and 2×2 hard squares on the square lattice depends on the composition of the mixture in exactly the manner predicted by the theory of Ashkin-Teller criticality, including in the hard-square limit. This result settles the question regarding the nature of the transition in the hard-square lattice gas. It also provides the first example of a polydisperse system whose critical properties depend on composition. Our ideas also lead to some interesting predictions for a class of frustrated quantum magnets that exhibit columnar ordering of the bond energies at low temperature.