Source author record

Andrea J. Liu

Andrea J. Liu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2022arXiv

Structuro-elasto-plasticity (StEP) model for plasticity in disordered solids

Elastoplastic lattice models for the response of solids to deformation typically incorporate structure only implicitly via a local yield strain that is assigned to each site. However, the local yield strain can change in response to a nearby or even distant plastic event in the system. This interplay is key to understanding phenomena such as avalanches in which one plastic event can trigger another, leading to a cascade of events, but typically is neglected in elastoplastic models. To include the interplay one could calculate the local yield strain for a given particulate system and follow its evolution, but this is expensive and requires knowledge of particle interactions, which is often hard to extract from experiments. Instead, we introduce a structural quantity, "softness," obtained using machine learning to correlate with imminent plastic rearrangements. We show that softness also correlates with local yield strain. We incorporate softness to construct a "structuro-elasto-plasticity" model that reproduces particle simulation results quantitatively for several observable quantities, confirming that we capture the influence of the interplay of local structure, plasticity, and elasticity on material response.

preprint2021arXiv

Free-then-freeze: transient learning degrees of freedom for introducing function in materials

The introduction of transient learning degrees of freedom into a system can lead to novel material design and training protocols that guide a system into a desired metastable state. In this approach, some degrees of freedom, which were not initially included in the system dynamics, are first introduced and subsequently removed from the energy minimization process once the desired state is reached. Using this conceptual framework, we create stable jammed packings that exist in exceptionally deep energy minima marked by the absence of low-frequency quasilocalized modes; this added stability persists in the thermodynamic limit. The inclusion of particle radii as transient degrees of freedom leads to deeper and much more stable minima than does the inclusion of particle stiffnesses. This is because particle radii couple to the jamming transition whereas stiffnesses do not. Thus different choices for the added degrees of freedom can lead to very different training outcomes.

preprint2021arXiv

Physical learning beyond the quasistatic limit

Physical networks, such as biological neural networks, can learn desired functions without a central processor, using local learning rules in space and time to learn in a fully distributed manner. Learning approaches such as equilibrium propagation, directed aging, and coupled learning similarly exploit local rules to accomplish learning in physical networks such as mechanical, flow, or electrical networks. In contrast to certain natural neural networks, however, such approaches have so far been restricted to the quasistatic limit, where they learn on time scales slow compared to their physical relaxation. This quasistatic constraint slows down learning, limiting the use of these methods as machine learning algorithms, and potentially restricting physical networks that could be used as learning platforms. Here we explore learning in an electrical resistor network that implements coupled learning, both in the lab and on the computer, at rates that range from slow to far above the quasistatic limit. We find that up to a critical threshold in the ratio of the learning rate to the physical rate of relaxation, learning speeds up without much change of behavior or error. Beyond the critical threshold, the error exhibits oscillatory dynamics but the networks still learn successfully.

preprint2021arXiv

Quantifying the link between local structure and cellular rearrangements using information in models of biological tissues

Machine learning techniques have been used to quantify the relationship between local structural features and variations in local dynamical activity in disordered glass-forming materials. To date these methods have been applied to an array of standard (Arrhenius and super-Arrhenius) glass formers, where work on "soft spots" indicates a connection between the linear vibrational response of a configuration and the energy barriers to non-linear deformations. Here we study the Voronoi model, which takes its inspiration from dense epithelial monolayers and which displays anomalous, sub-Arrhenius scaling of its dynamical relaxation time with decreasing temperature. Despite these differences, we find that the likelihood of rearrangements can vary by several orders of magnitude within the model tissue and extract a local structural quantity, "softness" that accurately predicts the temperature-dependence of the relaxation time. We use an information-theoretic measure to quantify the extent to which softness determines impending topological rearrangements; we find that softness captures nearly all of the information about rearrangements that is obtainable from structure, and that this information is large in the solid phase of the model and decreases rapidly as state variables are varied into the fluid phase.

preprint2020arXiv

Attractive versus truncated repulsive supercooled liquids: The dynamics is encoded in the pair correlation function

We compare glassy dynamics in two liquids that differ in the form of their interaction potentials. Both systems have the same repulsive interactions but one has also an attractive part in the potential. These two systems exhibit very different dynamics despite having nearly identical pair correlation functions. We demonstrate that a properly weighted integral of the pair correlation function, which amplifies the subtle differences between the two systems, correctly captures their dynamical differences. The weights are obtained from a standard machine learning algorithm.

preprint2020arXiv

Correlation of plastic events with local structure in jammed packings across spatial dimensions

In jammed packings, it is usually thought that local structure only plays a significant role in specific regimes. The standard deviation of the relative excess coordination, $σ_Z/ Z_\mathrm{c}$, decays like $1/\sqrt{d}$, so that local structure should play no role in high spatial dimensions. Furthermore, in any fixed dimension $d \geq 2$, there are diverging length scales as the pressure vanishes approaching the unjamming transition, again suggesting that local structure should not be sufficient to describe response. Here we challenge the assumption that local structure does not matter in these cases. In simulations of jammed packings under athermal, quasistatic shear, we use machine learning to identify a local structural variable, softness, that correlates with rearrangements in dimensions $d=2$ to $d=5$. We find that softness - and even just the coordination number $Z$ - are quite predictive of rearrangements over a wide range of pressures, all the way down to unjamming, in all $d$ studied. This result provides direct evidence that local structure can play a role in higher spatial dimensions.

preprint2020arXiv

Fluctuation Distributions of Energy Minima in Complex Landscapes

We discuss the properties of the distributions of energies of minima obtained by gradient descent in complex energy landscapes. We find strikingly similar phenomenology across several prototypical models. We particularly focus on the distribution of energies of minima in the analytically well-understood p-spin-interaction spin glass model. We numerically find non-Gaussian distributions that resemble the Tracy-Widom distributions often found in problems of random correlated variables, and non-trivial finite-size scaling. Based on this, we propose a picture of gradient descent dynamics that highlights the importance of a first-passage process in the eigenvalues of the Hessian. This picture provides a concrete link to problems in which the Tracy-Widom distribution is established. Aspects of this first-passage view of gradient-descent dynamics are generic for non-convex complex landscapes, rationalizing the commonality that we find across models.

preprint2020arXiv

Unifying framework for strong and fragile liquids via machine learning: a study of liquid silica

The fragility of a glassforming liquid characterizes how rapidly its relaxation dynamics slow down with cooling. The viscosity of strong liquids follows an Arrhenius law with a temperature-independent barrier height to rearrangements responsible for relaxation, whereas fragile liquids experience a much faster increase in their dynamics, suggesting a barrier height that increases with decreasing temperature. Strong glassformers are typically network glasses, while fragile glassformers are typically molecular or hard-sphere-like. As a result of these differences at the microscopic level, strong and fragile glassformers are usually treated separately from a theoretical point of view. Silica is the archetypal strong glassformer at low temperatures, but also exhibits a mysterious strong-to-fragile crossover at higher temperatures. Here we show that softness, a structure-based machine learned parameter that has previously been applied to fragile glassformers provides a useful description of model liquid silica in the strong and fragile regimes, and through the strong-to-fragile crossover. Just as for fragile glassformers, the relationship between softness and dynamics is invariant and Arrhenius in all regimes, but the average softness changes with temperature. The strong-to-fragile crossover in silica is not due to a sudden, qualitative change in structure, but can be explained by a simple Arrhenius form with a continuously and linearly changing local structure. Our results unify the study of liquid silica under a single simple conceptual picture.

preprint2019arXiv

Directed aging, memory and Nature's greed

Disordered materials are often out of equilibrium and evolve very slowly. This allows a memory of the imposed strains or preparation conditions to be encoded in the material. Here we consider "directed aging", where the elastic properties of a material evolve in the direction defined by its aging history. The evolution to a lower-energy configuration is controlled by steepest decent and affects stressed regions differently from unstressed ones. This process can be considered to be a "greedy algorithm" of Nature. Our experiments and simulations illustrate directed aging in examples in which the material's elasticity evolves as a direct consequence of the imposed deformation; the material itself decides how to evolve in order to produce responses that were not present inherently in the material.

preprint2019arXiv

Periodic training of creeping solids

We consider disordered solids in which the microscopic elements can deform plastically in response to stresses on them. We show that by driving the system periodically, this plasticity can be exploited to train in desired elastic properties, both in the global moduli and in local "allosteric" interactions. Periodic driving can couple an applied "source" strain to a target strain over a path in the energy landscape. This coupling allows control of the system's response even at large strains well into the nonlinear regime, where it can be difficult to achieve control simply by design.

preprint2019arXiv

Revealing structure-function relationships in functional flow networks via persistent homology

Complex networks encountered in biology are often characterized by significant structural diversity. Whether it be differences in the three-dimensional structure of allosteric proteins, or the variation among the micro-scale structures of organisms' cerebral vasculature systems, identifying relationships between structure and function often poses a difficult challenge. Here we showcase an approach to characterizing structure-function relationships in complex networks applied in the context of flow networks tuned to perform specific functions. Using persistent homology, we analyze flow networks tuned to perform complex multifunctional tasks, answering the question of how local changes in the network structure coordinate to create functionality at at the scale of the entire network. We find that the response of such networks encodes hidden topological features - sectors of uniform pressure - that are not apparent in the underlying network architectures, Regardless of differences in local connectivity, these features provide a universal topological description for all networks that perform these types of functions. We show that these features correlate strongly with the tuned response, providing a clear topological relationship between structure and function and structural insight into the limits of multifunctionality.

preprint2019arXiv

The hidden topological structure of flow network functionality

The ability to reroute and control flow is vital to the function of venation networks across a wide range of organisms. By modifying individual edges in these networks, either by adjusting edge conductances or creating and destroying edges, organisms can robustly control the propagation of inputs to perform specific tasks. However, a fundamental disconnect exists between the structure and function of these networks: networks with different local architectures can perform the same functions. Here we answer the question of how structural changes at the microscopic level are able to collectively create functionality at the scale of an entire network. Using persistent homology, we analyze networks tuned to perform complex multifunctional tasks. We find that the responses of such networks encode a hidden topological structure composed of sectors of uniform pressure. Although these sectors are not apparent in the underlying network architectures, we find that they nonetheless correlate strongly with the tuned function. We conclude that the connectivity of these sectors, rather than that of the individual nodes, provides a quantitative relationship between structure and function in flow networks. Finally, we use this topological description to place a bound on the limits of task complexity.

preprint2016arXiv

Pinning Susceptibility: The effect of dilute, quenched disorder on jamming

We study the effect of dilute pinning on the jamming transition. Pinning reduces the average contact number needed to jam unpinned particles and shifts the jamming threshold to lower densities, leading to a pinning susceptibility, $χ_p$. Our main results are that this susceptibility obeys scaling form and diverges in the thermodynamic limit as $χ_p \propto |ϕ- ϕ_c^\infty|^{-γ_p}$ where $ϕ_c^\infty$ is the jamming threshold in the absence of pins. Finite-size scaling arguments yield these values with associated statistical (systematic) errors $γ_p = 1.018 \pm 0.026 (0.291) $ in $d=2$ and $γ_p =1.534 \pm 0.120 (0.822)$ in $d=3$. Logarithmic corrections raise the exponent in $d=2$ to close to the $d=3$ value, although the systematic errors are very large.

preprint2016arXiv

Spatial structure of states of self stress in jammed systems

States of self stress, organizations of internal forces in many-body systems that are in equilibrium with an absence of external forces, can be thought of as the constitutive building blocks of the elastic response of a material. In overconstrained disordered packings they have a natural mathematical correspondence with the zero-energy vibrational modes in underconstrained systems. While substantial attention in the literature has been paid to diverging length scales associated with zero- and finite-energy vibrational modes in jammed systems, less is known about the spatial structure of the states of self stress. In this work we define a natural way in which a unique state of self stress can be associated with each bond in a disordered spring network derived from a jammed packing, and then investigate the spatial structure of these bond-localized states of self stress. This allows for an understanding of how the elastic properties of a system would change upon changing the strength or even existence of any bond in the system.

preprint2015arXiv

A random matrix definition of the boson peak

The density of vibrational states for glasses and jammed solids exhibits universal features, including an excess of modes above the Debye prediction known as the boson peak located at a frequency $ω^*$ . We show that the eigenvector statistics for boson peak modes are universal, and develop a new definition of the boson peak based on this universality that displays the previously observed characteristic scaling $ω^*\sim p^{-1/2}$ . We identify a large new class of random matrices that obey a generalized global tranlational invariance constraint and demonstrate that members of this class also have a boson peak with precisely the same universal eigenvector statistics. We denote this class as boson peak random matrices, and conjecture it comprises a new universality class. We characterize the eigenvector statistics as a function of coordination number, and find that one member of this new class reproduces the scaling of $ω^{*}$ with coordination number that is observed near the jamming transition.

preprint2015arXiv

Deformation-Driven Diffusion and Plastic Flow in Two-Dimensional Amorphous Granular Pillars

We report a combined experimental and simulation study of deformation-induced diffusion in compacted two-dimensional amorphous granular pillars, in which thermal fluctuations play negligible role. The pillars, consisting of bidisperse cylindrical acetal plastic particles standing upright on a substrate, are deformed uniaxially and quasistatically by a rigid bar moving at a constant speed. The plastic flow and particle rearrangements in the pillars are characterized by computing the best-fit affine transformation strain and non-affine displacement associated with each particle between two stages of deformation. The non-affine displacement exhibits exponential crossover from ballistic to diffusive behavior with respect to the cumulative deviatoric strain, indicating that in athermal granular packings, the cumulative deviatoric strain plays the role of time in thermal systems and drives effective particle diffusion. We further study the size-dependent deformation of the granular pillars by simulation, and find that different-sized pillars follow self-similar shape evolution during deformation. In addition, the yield stress of the pillars increases linearly with pillar size. Formation of transient shear lines in the pillars during deformation becomes more evident as pillar size increases. The width of these elementary shear bands is about twice the diameter of a particle, and does not vary with pillar size.

preprint2015arXiv

Divergence of Voronoi cell anisotropy vector: A threshold-free characterization of local structure in amorphous materials

Characterizing structural inhomogeneity is an essential step in understanding the mechanical response of amorphous materials. We introduce a threshold-free measure based on the field of vectors pointing from the center of each particle to the centroid of the Voronoi cell in which the particle resides. These vectors tend to point in toward regions of high free volume and away from regions of low free volume, reminiscent of sinks and sources in a vector field. We compute the local divergence of these vectors, for which positive values correspond to overpacked regions and negative values identify underpacked regions within the material. Distributions of this divergence are nearly Gaussian with zero mean, allowing for structural characterization using only the moments of the distribution. We explore how the standard deviation and skewness vary with packing fraction for simulations of bidisperse systems and find a kink in these moments that coincides with the jamming transition.

preprint2015arXiv

Scaling theory for the jamming transition

We propose a scaling ansatz for the elastic energy of a system near the critical jamming transition in terms of three relevant fields: the compressive strain $Δϕ$ relative to the critical jammed state, the shear strain $ε$, and the inverse system size $1/N$. We also use $ΔZ$, the number of contacts relative to the minimum required at jamming, as an underlying control parameter. Our scaling theory predicts new exponents, exponent equalities and scaling collapses for energy, pressure and shear stress that we verify with numerical simulations of jammed packings of soft spheres. It also yields new insight into why the shear and bulk moduli exhibit different scalings; the difference arises because the shear stress vanishes as $1/\sqrt{N}$ while the pressure approaches a constant in the thermodynamic limit. The success of the scaling ansatz implies that the jamming transition exhibits an emergent scale invariance, and that it should be possible to develop a renormalization-group theory for jamming.

preprint2015arXiv

Strain fluctuations and elastic moduli in disordered solids

Recently there has been a surge in interest in using video-microscopy techniques to infer the local mechanical properties of disordered solids. One common approach is to minimize the difference between particle vibrational displacements in a local coarse-graining volume and the displacements that would result from a best-fit affine deformation. Effective moduli are then be inferred under the assumption that the components of this best-fit affine deformation tensor have a Boltzmann distribution. In this paper, we combine theoretical arguments with experimental and simulation data to demonstrate that the above does not reveal information about the true elastic moduli of jammed packings and colloidal glasses.

preprint2015arXiv

Tuning by pruning: exploiting disorder for global response and the principle of bond-level independence

We exploit the intrinsic difference between disordered and crystalline solids to create systems with unusual and exquisitely tuned mechanical properties. To demonstrate the power of this approach, we design materials that are either virtually incompressible or completely auxetic. Disordered networks can be efficiently driven to these extreme limits by removing a very small fraction of bonds via a selected-bond removal procedure that is both simple and experimentally relevant. The procedure relies on the nearly complete absence of any correlation between the contributions of an individual bond to different elastic moduli. A new principle unique to disordered solids underlies this lack of correlation: independence of bond-level response.

preprint2014arXiv

Comment on 'Repulsive Contact Interactions Make Jammed Particulate Systems Inherently Nonharmonic'

In their Letter, Schreck, Bertrand, O'Hern and Shattuck [Phys. Rev. Lett. 107, 078301 (2011)] study nonlinearities in jammed particulate systems that arise when contacts are altered. They conclude that there is "no harmonic regime in the large system limit for all compressions" and "at jamming onset for any system size." Their argument rests on the claim that for finite-range repulsive potentials, of the form used in studies of jamming, the breaking or forming of a single contact is sufficient to destroy the linear regime. We dispute these conclusions and argue that linear response is both justified and essential for understanding the nature of the jammed solid.

preprint2014arXiv

Contact nonlinearities and linear response in jammed particulate packings

Packings of frictionless athermal particles that interact only when they overlap experience a jamming transition as a function of packing density. Such packings provide the foundation for the theory of jamming. This theory rests on the observation that, despite the multitude of disordered configurations, the mechanical response to linear order depends only on the distance to the transition. We investigate the validity and utility of such measurements that invoke the harmonic approximation and show that, despite particles coming in and out of contact, there is a well-defined linear regime in the thermodynamic limit.

preprint2014arXiv

Disordered surface vibrations in jammed sphere packings

We study the vibrational properties near a free surface of disordered spring networks derived from jammed sphere packings. In bulk systems, without surfaces, it is well understood that such systems have a plateau in the density of vibrational modes extending down to a frequency scale $ω^*$. This frequency is controlled by $ΔZ = \langle Z \rangle - 2d$, the difference between the average coordination of the spheres and twice the spatial dimension, $d$, of the system, which vanishes at the jamming transition. In the presence of a free surface we find that there is a density of disordered vibrational modes associated with the surface that extends far below $ω^*$. The total number of these low-frequency surface modes is controlled by $ΔZ$, and the profile of their decay into the bulk has two characteristic length scales, which diverge as $ΔZ^{-1/2}$ and $ΔZ^{-1}$ as the jamming transition is approached.

preprint2014arXiv

Identifying structural flow defects in disordered solids using machine learning methods

We use machine learning methods on local structure to identify flow defects - or regions susceptible to rearrangement - in jammed and glassy systems. We apply this method successfully to two disparate systems: a two dimensional experimental realization of a granular pillar under compression, and a Lennard-Jones glass in both two and three dimensions above and below its glass transition temperature. We also identify characteristics of flow defects that differentiate them from the rest of the sample. Our results show it is possible to discern subtle structural features responsible for heterogeneous dynamics observed across a broad range of disordered materials.

preprint2014arXiv

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. The situation for jamming is surprisingly rich: the assumption that jammed packings are isotropic is only strictly true in the large-size limit, and finite-size has a profound effect on the very meaning of jamming. Here, we provide a comprehensive numerical study of finite-size effects in sphere packings above the jamming transition, focusing on stability as well as the scaling of the contact number and the elastic response.

preprint2014arXiv

Predicting plasticity with soft vibrational modes: from dislocations to glasses

We show that quasi localized low-frequency modes in the vibrational spectrum can be used to construct soft spots, or regions vulnerable to rearrangement, which serve as a universal tool for the identification of flow defects in solids. We show that soft spots not only encode spatial information, via their location, but also directional information, via directors for particles within each soft spot. Single crystals with isolated dislocations exhibit low-frequency phonon modes that localize at the core, and their polarization pattern predicts the motion of atoms during elementary dislocation glide in exquisite detail. Even in polycrystals and disordered solids, we find that the directors associated with particles in soft spots are highly correlated with the direction of particle displacements in rearrangements.

preprint2014arXiv

Understanding plastic deformation in thermal glasses from single-soft-spot dynamics

By considering the low-frequency vibrational modes of amorphous solids, Manning and Liu [Phys. Rev. Lett. 107, 108302 (2011)] showed that a population of "soft spots" can be identified that are intimately related to plasticity at zero temperature under quasistatic shear. In this work we track individual soft spots with time in a two-dimensional sheared thermal Lennard Jones glass at temperatures ranging from deep in the glassy regime to above the glass transition temperature. We show that the lifetimes of individual soft spots are correlated with the timescale for structural relaxation. We additionally calculate the number of rearrangements required to destroy soft spots, and show that most soft spots can survive many rearrangements. Finally, we show that soft spots are robust predictors of rearrangements at temperatures well into the super-cooled regime. Altogether, these results pave the way for mesoscopic theories of plasticity of amorphous solids based on dynamical behavior of individual soft spots.

preprint2013arXiv

Mechanical signaling via nonlinear wavefront propagation in a mechanically-excitable medium

Models that invoke nonlinear wavefront propagation in a chemically excitable medium are rife in the biological literature. Indeed, the idea that wavefront propagation can serve as a signaling mechanism has often been invoked to explain synchronization of developmental processes. In this paper we suggest a new kind of signaling based not on diffusion of a chemical species but on the propagation of mechanical stress. We construct a theoretical approach to describe mechanical signaling as a nonlinear wavefront propagation problem and study its dependence on key variables such as the effective elasticity and damping of the medium.

preprint2013arXiv

Phonon Dispersion and Elastic Moduli of Two-Dimensional Disordered Colloidal Packings of Soft Particles with Frictional Interactions

Particle tracking and displacement covariance matrix techniques are employed to investigate the phonon dispersion relations of two-dimensional colloidal glasses composed of soft, thermoresponsive microgel particles whose temperature-sensitive size permits \textit{in situ} variation of particle packing fraction. Bulk, $B$, and shear, $G$, moduli of the colloidal glasses are extracted from the dispersion relations as a function of packing fraction, and variation of the ratio $G/B$ with packing fraction is found to agree quantitatively with predictions for jammed packings of frictional soft particles. In addition, $G$ and $B$ individually agree with numerical predictions for frictional particles. This remarkable level of agreement enabled us to extract an energy scale for the inter-particle interaction from the individual elastic constants and to derive an approximate estimate for the inter-particle friction coefficient.

preprint2013arXiv

Phonons in two-dimensional soft colloidal crystals

The vibrational modes of pristine and polycrystalline monolayer colloidal crystals composed of thermosensitive microgel particles are measured using video microscopy and covariance matrix analysis. At low frequencies, the Debye relation for two dimensional harmonic crystals is observed in both crystal types; at higher frequencies, evidence for van Hove singularities in the phonon density of states is significantly smeared out by experimental noise and measurement statistics. The effects of these errors are analyzed using numerical simulations. We introduce methods to correct for these limitations, which can be applied to disordered systems as well as crystalline ones, and we show that application of the error correction procedure to the experimental data leads to more pronounced van Hove singularities in the pristine crystal. Finally, quasi-localized low-frequency modes in polycrystalline two-dimensional colloidal crystals are identified and demonstrated to correlate with structural defects such as dislocations, suggesting that quasi-localized low-frequency phonon modes may be used to identify local regions vulnerable to rearrangements in crystalline as well as amorphous solids.

preprint2013arXiv

Stability of jammed packings I: the rigidity length scale

In 2005, Wyart et al. (Europhys. Lett., 72 (2005) 486) showed that the low frequency vibrational properties of jammed amorphous sphere packings can be understood in terms of a length scale, called l*, that diverges as the system becomes marginally unstable. Despite the tremendous success of this theory, it has been difficult to connect the counting argument that defines l* to other length scales that diverge near the jamming transition. We present an alternate derivation of l* based on the onset of rigidity. This phenomenological approach reveals the physical mechanism underlying the length scale and is relevant to a range of systems for which the original argument breaks down. It also allows us to present the first direct numerical measurement of l*.

preprint2013arXiv

Stability of jammed packings II: the transverse length scale

As a function of packing fraction at zero temperature and applied stress, an amorphous packing of spheres exhibits a jamming transition where the system is sensitive to boundary conditions even in the thermodynamic limit. Upon further compression, the system should become insensitive to boundary conditions provided it is sufficiently large. Here we explore the linear response to a large class of boundary perturbations in 2 and 3 dimensions. We consider each finite packing with periodic-boundary conditions as the basis of an infinite square or cubic lattice and study properties of vibrational modes at arbitrary wave vector. We find that the stability of such modes be understood in terms of a competition between plane waves and the anomalous vibrational modes associated with the jamming transition; infinitesimal boundary perturbations become irrelevant for systems that are larger than a length scale that characterizes the transverse excitations. This previously identified length diverges at the jamming transition.

preprint2013arXiv

The syncytial Drosophila embryo as a mechanically excitable medium

Mitosis in the early syncytial Drosophila embryo is highly correlated in space and time, as manifested in mitotic wavefronts that propagate across the embryo. In this paper we investigate the idea that the embryo can be considered a mechanically-excitable medium, and that mitotic wavefronts can be understood as nonlinear wavefronts that propagate through this medium. We study the wavefronts via both image analysis of confocal microscopy videos and theoretical models. We find that the mitotic waves travel across the embryo at a well-defined speed that decreases with replication cycle. We find two markers of the wavefront in each cycle, corresponding to the onsets of metaphase and anaphase. Each of these onsets is followed by displacements of the nuclei that obey the same wavefront pattern. To understand the mitotic wavefronts theoretically we analyze wavefront propagation in excitable media. We study two classes of models, one with biochemical signaling and one with mechanical signaling. We find that the dependence of wavefront speed on cycle number is most naturally explained by mechanical signaling, and that the entire process suggests a scenario in which biochemical and mechanical signaling are coupled.

preprint2012arXiv

Finite-Size Scaling at the Jamming Transition

We present an analysis of finite-size effects in jammed packings of N soft, frictionless spheres at zero temperature. There is a 1/N correction to the discrete jump in the contact number at the transition so that jammed packings exist only above isostaticity. As a result, the canonical power-law scalings of the contact number and elastic moduli break down at low pressure. These quantities exhibit scaling collapse with a non-trivial scaling function, demonstrating that the jamming transition can be considered a phase transition. Scaling is achieved as a function of N in both 2 and 3 dimensions, indicating an upper critical dimension of 2.

preprint2011arXiv

Direct determination of the size of basins of attraction of jammed solids

We propose a free-energy based Monte-Carlo method to measure the volume of potential-energy basins in configuration space. Using this approach we can estimate the number of distinct potential-energy minima, even when this number is much too large to be sampled directly. We validate our approach by comparing our results with the direct enumeration of distinct jammed states in small packings of frictionless spheres. We find that the entropy of distinct packings is extensive and that the entropy of distinct hard-sphere packings must have a maximum as a function of packing fraction.

preprint2011arXiv

Mapping the glassy dynamics of soft spheres onto hard-sphere behavior

We show that the dynamics of soft-sphere systems with purely repulsive interactions can be described by introducing an effective hard-sphere diameter determined using the Andersen-Weeks-Chandler approximation. We find that this approximation, known to describe static properties of liquids, also gives a good description of a dynamical quantity, the relaxation time, even in the vicinity of the glass transition.

preprint2011arXiv

Measurement of correlations between low-frequency vibrational modes and particle rearrangements in quasi-two-dimensional colloidal glasses

We investigate correlations between low-frequency vibrational modes and rearrangements in two-dimensional colloidal glasses composed of thermosensitive microgel particles which readily permit variation of sample packing fraction. At each packing fraction, the particle displacement covariance matrix is measured and used to extract the vibrational spectrum of the "shadow" colloidal glass (i.e., the particle network with the same geometry and interactions as the sample colloid but absent damping). Rearrangements are induced by successive, small reductions in packing fraction. The experimental results suggest that low-frequency quasi-localized phonon modes in colloidal glasses, i.e., modes that present low energy barriers for system rearrangements, are spatially correlated with rearrangements in this thermal system.

preprint2011arXiv

Self-Diffusiophoresis in the Advection Dominated Regime

In both biological and artificial systems, concentration gradients can serve as a convenient mechanism for manipulating particles and generating motility. Particles that interact with a solute will move along its gradient; if they themselves generate the gradient, this mechanism provides a means of self-propulsion. We consider a version of this type of motility appropriate to certain biological systems where polymeric filaments provide the concentration gradient. As the filament diffusion is small, this corresponds to a regime of large Péclet number where the motion is dominated by the effects of fluid advection. The nature of such concentration-gradient-driven motion in the advective regime differs in certain fundamental respects from the same process at low Péclet number. In particular, we show that out of four broad scenarios of steady state motion at low Péclet number, only two remain viable in the strongly advecting limit.

preprint2011arXiv

The jamming scenario - an introduction and outlook

The jamming scenario of disordered media, formulated about 10 years ago, has in recent years been advanced by analyzing model systems of granular media. This has led to various new concepts that are increasingly being explored in in a variety of systems. This chapter contains an introductory review of these recent developments and provides an outlook on their applicability to different physical systems and on future directions. The first part of the paper is devoted to an overview of the findings for model systems of frictionless spheres, focussing on the excess of low-frequency modes as the jamming point is approached. Particular attention is paid to a discussion of the cross-over frequency and length scales that govern this approach. We then discuss the effects of particle asphericity and static friction, the applicability to bubble models for wet foams in which the friction is dynamic, the dynamical arrest in colloids, and the implications for molecular glasses.

preprint2010arXiv

Electrostatic cluster formation in lipid monolayers

We study phase separation in mixed monolayers of neutral and highly negatively charged lipids, induced by the addition of divalent positively charged counterions. We find good agreement between experiments on mixtures of pip2 and sopc and simulations of a simplified model in which only the essential electrostatic interactions are retained. Thus, our results support an interpretation of pip2 clustering as governed primarily by electrostatic interactions, in which divalent ions such as calcium mediate an effective attraction between like-charged lipids. Surprisingly, the mediated attractions are strong enough to give nearly complete phase separation, so that clusters can even form when the overall concentration of pip2 is low, as is the case in the cell membrane.

preprint2010arXiv

Far-From-Equilibrium Physics: An Overview

Isolated systems tend to evolve towards equilibrium, a special state that has been the focus of many-body research for a century. Yet much of the richness of the world around us arises from conditions far from equilibrium. Phenomena such as turbulence, earthquakes, fracture, and life itself occur only far from equilibrium. Subjecting materials to conditions far from equilibrium leads to otherwise unattainable properties. For example, rapid cooling is a key process in manufacturing the strongest metallic alloys and toughest plastics. Processes that occur far from equilibrium also create some of the most intricate structures known, from snowflakes to the highly organized structures of life. While much is understood about systems at or near equilibrium, we are just beginning to uncover the basic principles governing systems far from equilibrium.

preprint2010arXiv

Jamming Mechanisms and Density Dependence in a Kinetically-Constrained Model

We add relaxation mechanisms that mimic the effect of temperature and non-equilibrium driving to the recently-proposed spiral model which jams at a critical density rho_c < 1. This enables us to explore unjamming by temperature or driving at rho_c < rho < 1. We numerically calculate the relaxation time of the persistence function and its spatial heterogeneity. We disentangle the three different relaxation mechanisms responsible for unjamming when varying density, temperature, and driving strength, respectively. We show that the spatial scale of dynamic heterogeneity depends on density much more strongly than on temperature and driving.

preprint2010arXiv

Low-frequency vibrations of soft colloidal glasses

We conduct experiments on two-dimensional packings of colloidal thermosensitive hydrogel particles whose packing fraction can be tuned above the jamming transition by varying the temperature. By measuring displacement correlations between particles, we extract the vibrational properties of a corresponding "shadow" system with the same configuration and interactions, but for which the dynamics of the particles are undamped. The vibrational spectrum and the nature of the modes are very similar to those predicted for zero-temperature idealized sphere models and found in atomic and molecular glasses; there is a boson peak at low frequency that shifts to higher frequency as the system is compressed above the jamming transition.

preprint2010arXiv

Rotational and Translational Phonon Modes in Glasses Composed of Ellipsoidal Particles

The effects of particle shape on the vibrational properties of colloidal glasses are studied experimentally. 'Ellipsoidal glasses' are created by stretching polystyrene spheres to different aspect ratios and then suspending the resulting ellipsoidal particles in water at high packing fraction. By measuring displacement correlations between particles, we extract vibrational properties of the corresponding "shadow" ellipsoidal glass with the same geometric configuration and interactions as the 'source' suspension but without damping. Low frequency modes in glasses composed of ellipsoidal particles with major/minor axis aspect ratios $\sim$1.1 are observed to have predominantly rotational character. By contrast, low frequency modes in glasses of ellipsoidal particles with larger aspect ratios ($\sim$3.0) exhibit a mix of rotational and translational character. All glass samples were characterized by a distribution of particles with different aspect ratios. Interestingly, even within the same sample it was found that small-aspect-ratio particles participate relatively more in rotational modes, while large-aspect-ratio particles tend to participate relatively more in translational modes.

preprint2010arXiv

Universal Jamming Phase Diagram in the Hard-Sphere Limit

We present a new formulation of the jamming phase diagram for a class of glass-forming fluids consisting of spheres interacting via finite-ranged repulsions at temperature $T$, packing fraction $ϕ$ or pressure $p$, and applied shear stress $Σ$. We argue that the natural choice of axes for the phase diagram are the dimensionless quantities $T/pσ^3$, $pσ^3/ε$, and $Σ/p$, where $T$ is the temperature, $p$ is the pressure, $Σ$ is the stress, $σ$ is the sphere diameter, $ε$ is the interaction energy scale, and $m$ is the sphere mass. We demonstrate that the phase diagram is universal at low $pσ^3/ε$; at low pressure, observables such as the relaxation time are insensitive to details of the interaction potential and collapse onto the values for hard spheres, provided the observables are non-dimensionalized by the pressure. We determine the shape of the jamming surface in the jamming phase diagram, organize previous results in relation to the jamming phase diagram, and discuss the significance of various limits.

preprint2010arXiv

Vibrational modes identify soft spots in a sheared disordered packing

We analyze low-frequency vibrational modes in a two-dimensional, zero-temperature, quasistatically sheared model glass to identify a population of structural "soft spots" where particle rearrangements are initiated. The population of spots evolves slowly compared to the interval between particle rearrangements and the soft spots are structurally different from the rest of the system. Our results suggest that disordered solids flow via localized rearrangements that tend to occur at soft spots, which are analogous to dislocations in crystalline solids.

preprint2009arXiv

Anharmonicity and quasi-localization of the excess low-frequency vibrations in jammed solids

We compare the harmonic and anharmonic properties of the vibrational modes in 3-dimensional jammed packings of frictionless spheres interacting via repulsive, finite range potentials. A crossover frequency is apparent in the density of states, the diffusivity and the participation ratio of the modes. At this frequency, which shifts to zero at the jamming threshold, the vibrational modes have a very small participation ratio implying that the modes are quasi-localized. The most anharmonic modes occur at low frequency which is opposite to what is normally found in crystals. The lowest frequency modes have the strongest response to the pressure and the lowest energy barriers to mechanical failure.

preprint2009arXiv

Heat transport in model jammed solids

We calculate numerically the normal modes of vibrations in 3D jammed packings of soft spheres as a function of the packing fraction and obtain the energy diffusivity, a spectral measure of transport that controls sound propagation and thermal conductivity. The crossover frequency between weak and strong phonon scattering is controlled by the coordination and shifts to zero as the system is decompressed towards the critical packing fraction at which rigidity is lost. Below the crossover, the diffusivity displays a power-law divergence with inverse frequency, which suggests that the vibrational modes are primarily transverse waves, weakly scattered by disorder. Above it, a large number of modes appear whose diffusivity plateaus at a nearly constant value independent of the inter-particle potential, before dropping to zero above the Anderson localization frequency. The thermal conductivity of a marginally jammed solid just above the rigidity threshold is calculated and related to the one measured experimentally at room temperature for most glasses.

preprint2008arXiv

Normal Modes in Model Jammed Systems in Three Dimensions

Vibrational spectra and normal modes of mechanically stable particle packings in three dimensions are analyzed over a range of compressions, from near the jamming transition, where the packings lose their rigidity, to far above it. At high frequency, the normal modes are localized at all compressions. At low frequency, the nature of the modes depends somewhat on compression. At large compressions, far from the transition, the lowest-frequency normal modes have some plane-wave character, though less than one would expect for a crystalline or isotropic solid. At low compressions near the jamming transition, the lowest-frequency modes are neither plane-wave-like nor localized. We characterize these differences, highlighting the unusual dispersion behavior that emerges for marginally jammed solids.

preprint2007arXiv

Activated dynamics and effective temperature in a steady state sheared glass

We conduct nonequilibrium molecular dynamics simulations to measure the shear stress, the average inherent structure energy, and the effective temperature $T_{eff}$ of a sheared model glass as a function of bath temperature $T$ and shear strain rate. For $T$ above the glass transition temperature $T_0$, the rheology approaches a Newtonian limit and $T_{eff}$ approaches $T$ as the strain rate approaches zero, while for $T<T_0$, the shear stress approaches a yield stress and $T_{eff}$ approaches a limiting value near $T_0$. In the shear-dominated regime at high $T$, high strain rate or at low $T$, we find that the shear stress and the average inherent structure energy each collapse onto a single curve as a function of $T_{eff}$. This indicates that $T_{eff}$ is controlling behavior in this regime.

preprint2005arXiv

Structural Polymorphism of the Cytoskeleton: A Model of Linker-Assisted Filament Aggregation

The phase behavior of charged rods in the presence of inter-rod linkers is studied theoretically as a model for the equilibrium behavior underlying the organization of actin filaments by linker proteins in the cytoskeleton. The presence of linkers in the solution modifies the effective inter-rod interaction and can lead to inter-filament attraction. Depending on the system's composition and physical properties such as linker binding energies, filaments will either orient perpendicular or parallel to each other, leading to network-like or bundled structures. We show that such a system can have one of three generic phase diagrams, one dominated by bundles, another by networks, and the third containing both bundle and network-like phases. The first two diagrams can be found over a wide range of interaction energies, while the third occurs only for a narrow range. These results provide theoretical understanding of the classification of linker proteins as bundling proteins or crosslinking proteins. In addition, they suggest possible mechanisms by which the cell may control cytoskeletal morphology.