Emergent disorder and mechanical memory in periodic metamaterials.

Ordered mechanical systems typically have one or only a few stable rest configurations, and hence are not considered useful for encoding memory. Multistable and history-dependent responses usually emerge from quenched disorder, for example in amorphous solids or crumpled sheets. In contrast, due to geometric frustration, periodic magnetic systems can create their own disorder and espouse an extensive manifold of quasi-degenerate configurations. Inspired by the topological structure of frustrated artificial spin ices, we introduce an approach to design ordered, periodic mechanical metamaterials that exhibit an extensive set of spatially disordered states. While our design exploits the correspondence between frustration in magnetism and incompatibility in meta-mechanics, our mechanical systems encompass continuous degrees of freedom, and thus generalize their magnetic counterparts. We show how such systems exhibit non-Abelian and history-dependent responses, as their state can depend on the order in which external manipulations were applied. We demonstrate how this richness of the dynamics enables to recognize, from a static measurement of the final state, the sequence of operations that an extended system underwent. Thus, multistability and potential to perform computation emerge from geometric frustration in ordered mechanical lattices that create their own disorder.

Bloch oscillations, Landau-Zener transition, and topological phase evolution in an array of coupled pendula.

We experimentally and theoretically study the dynamics of a one-dimensional array of pendula with a mild spatial gradient in their self-frequency and where neighboring pendula are connected with weak and alternating coupling. We map their dynamics to the topological Su-Schrieffer-Heeger model of charged quantum particles on a lattice with alternating hopping rates in an external electric field. By directly tracking the dynamics of a wave-packet in the bulk of the lattice, we observe Bloch oscillations, Landau-Zener transitions, and coupling between the isospin (i.e., the inner wave function distribution within the unit cell) and the spatial degrees of freedom (the distribution between unit cells). We then use Bloch oscillations in the bulk to directly measure the nontrivial global topological phase winding and local geometric phase of the band. We measure an overall evolution of 3.1 [Formula: see text] 0.2 radians for the geometrical phase during the Bloch period, consistent with the expected Zak phase of [Formula: see text]. Our results demonstrate the power of classical analogs of quantum models to directly observe the topological properties of the band structure and shed light on the similarities and the differences between quantum and classical topological effects.

Many-body interactions between contracting living cells.

The organization of live cells into tissues and their subsequent biological function involves inter-cell mechanical interactions, which are mediated by their elastic environment. To model this interaction, we consider cells as spherical active force dipoles surrounded by an unbounded elastic matrix. Even though we assume that this elastic medium responds linearly, each cell's regulation of its mechanical activity leads to nonlinearities in the emergent interactions between cells. We study the many-body nature of these interactions by considering several geometries that include three or more cells. We show that for different regulatory behaviors of the cells' activity, the total elastic energy stored in the medium differs from the superposition of all two-body interactions between pairs of cells within the system. Specifically, we find that the many-body interaction energy between cells that regulate their position is smaller than the sum of interactions between all pairs of cells in the system, while for cells that do not regulate their position, the many-body interaction is larger than the superposition prediction. Thus, such higher-order interactions should be considered when studying the mechanics of multiple cells in proximity.

Percolation in Networks of Liquid Diodes.

Liquid diodes are surface structures that facilitate the spontaneous flow of liquids in a specific direction. In nature, they are used to increase water collection and uptake, reproduction, and feeding. However, large networks with directional properties are exceptional and are typically limited up to a few centimeters. Here, we simulate, design, and 3D print liquid diode networks consisting of hundreds of unit cells. We provide structural and wettability guidelines for directional transport of liquids through these networks and introduce percolation theory in order to identify the threshold between a connected network, which allows fluid to reach specific points, and a disconnected network. By constructing well-defined networks with uni- and bidirectional pathways, we experimentally demonstrate the applicability of models describing isotropically directed percolation. We accurately predict the network permeability and the liquid final state. These guidelines are highly promising for the development of structures for spontaneous, yet predictable, directional liquid transport.

Mean-field interactions between living cells in linear and nonlinear elastic matrices.

Living cells respond to mechanical changes in the matrix surrounding them by applying contractile forces that are in turn transmitted to distant cells. We consider simple effective geometries for the spatial arrangement of cells, we calculate the mechanical work that each cell performs in order to deform the matrix, and study how that energy changes when a contracting cell is surrounded by other cells with similar properties and behavior. Cells regulating the displacements that they generate are attracted to each other in a manner that does not depend on the cell's rigidity. Whereas cells regulating the active stress that they apply repel each other. This repulsion depends on the cell's bulk modulus in spherical geometry, while in cylindrical geometries the interaction depends also on their shear modulus. In nonlinear, strain-stiffening matrices, for displacement regulation, in the presence of other cells, cell contraction is limited due to the divergence of the shear stress. For stress regulation, the interaction energy drops at the nonlinear stiffening stress. Our theoretical work provides insight into matrix-mediated interactions between contractile cells and on the role of their mechanical regulatory behavior.

Topology Restricts Quasidegeneracy in Sheared Square Colloidal Ice.

Recovery of ground-state degeneracy in two-dimensional square ice is a significant challenge in the field of geometric frustration with far-reaching fundamental implications, such as realization of vertex models and understanding the effect of dimensionality reduction. We combine experiments, theory, and numerical simulations to demonstrate that sheared square colloidal ice partially recovers the ground-state degeneracy for a wide range of field strengths and lattice shear angles. Our method could inspire engineering a novel class of frustrated microstructures and nanostructures based on sheared magnetic lattices in a wide range of soft- and condensed-matter systems.

Non-Newtonian Topological Mechanical Metamaterials Using Feedback Control.

We introduce a method to design topological mechanical metamaterials that are not constrained by Newtonian dynamics. The unit cells in a mechanical lattice are subjected to active feedback forces that are processed through autonomous controllers preprogrammed to generate the desired local response in real time. As an example, we focus on the quantum Haldane model, which is a two-band system with nonreciprocal coupling terms, the implementation of which in mechanical systems requires violating Newton's third law. We demonstrate that the required topological phase characterized by chiral edge modes can be achieved in an analogous mechanical system only with closed-loop control. We then show that our approach enables us to realize, a modified version of the Haldane model in a mechanical metamaterial. Here, the complex-valued couplings are polarized in a way that modes on opposite edges of a lattice propagate in the same direction, and are balanced by counterpropagating bulk modes. The proposed method is general and flexible, and could be used to realize arbitrary lattice parameters, such as nonlocal or nonlinear couplings, time-dependent potentials, non-Hermitian dynamics, and more, on a single platform.

Constraint relaxation leads to jamming.

Adding transitions to an equilibrium system increases the activity. Naively, one would expect this to hold also in out-of-equilibrium systems. We demonstrate, using relatively simple models, how adding transitions to an out of equilibrium system may in fact reduce the activity and even cause it to vanish. This surprising effect is caused by adding heretofore forbidden transitions into less and less active states. We investigate six related kinetically constrained lattice gas models, some of which behave as naively expected while others exhibit this nonintuitive behavior. These models exhibit an absorbing state phase transition, which is also affected by the added transitions. We introduce a semi-mean-field approximation describing the models, which agrees qualitatively with our numerical simulation.

Kinetically constrained model for gravity-driven granular flow and clogging.

We add extreme driving to the Kob-Andersen kinetically constrained lattice-gas model in order to mimic the effect of gravity on dense granular systems. For low particle densities, the current that develops in the system agrees at arbitrary field intensity with a mean-field theory. At intermediate densities, spatial correlations give rise to nonmonotonic dependence of the current on field intensity. At higher densities, the current ultimately vanishes at a finite, field-dependent density. We supplement the study of this bulk behavior with an investigation of the current through a narrow hole. There, lateral flow decreases the local density in front of the hole. Remarkably, the current through the hole quantitatively agrees with a theoretical prediction based on the bulk current at the measured local density.

Elastic interactions between anisotropically contracting circular cells.

We study interactions between biological cells that apply anisotropic active mechanical forces on an elastic substrate. We model the cells as thin disks that along their perimeters apply radial, but angle-dependent forces on the substrate. We obtain analytical expressions for the elastic energy stored in the substrate as a function of the distance between the cells, the Fourier modes of applied forces, and their phase angles. We show how the relative phases of the forces applied by the cells can switch the interaction between attractive and repulsive, and relate our results to those for linear force dipoles. For long enough distances, the interaction energy decays in magnitude as a power law of the cell-cell distance with an integer exponent that generally increases with the Fourier modes of the applied forces.

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.

Mechanical Interaction between Cells Facilitates Molecular Transport.

In vivo, eukaryotic cells are embedded in a matrix environment, where they grow and develop. Generally, this extracellular matrix (ECM) is an anisotropic fibrous structure, through which macromolecules and biochemical signaling molecules at the nanometer scale diffuse. The ECM is continuously remodeled by cells, via mechanical interactions, which lead to a potential link between biomechanical and biochemical cell-cell interactions. Here, it is studied how cell-induced forces applied on the ECM impact the biochemical transport of molecules between distant cells. It is experimentally observed that cells remodel the ECM by increasing fiber alignment and density of the matrix between them over time. Using random walk simulations on a 3D lattice, elongated fixed obstacles are implemented that mimic the fibrous ECM structure. Both diffusion of a tracer molecule and the mean first-passage time a molecule secreted from one cell takes to reach another cell are measured. The model predicts that cell-induced remodeling can lead to a dramatic speedup in the transport of molecules between cells. Fiber alignment and densification cause reduction of the transport dimensionality from a 3D to a much more rapid 1D process. Thus, a novel mechanism of mechano-biochemical feedback in the regulation of long-range cell-cell communication is suggested.

Attraction Controls the Entropy of Fluctuations in Isosceles Triangular Networks.

We study two-dimensional triangular-network models, which have degenerate ground states composed of straight or randomly-zigzagging stripes and thus sub-extensive residual entropy. We show that attraction is responsible for the inversion of the stable phase by changing the entropy of fluctuations around the ground-state configurations. By using a real-space shell-expansion method, we compute the exact expression of the entropy for harmonic interactions, while for repulsive harmonic interactions we obtain the entropy arising from a limited subset of the system by numerical integration. We compare these results with a three-dimensional triangular-network model, which shows the same attraction-mediated selection mechanism of the stable phase, and conclude that this effect is general with respect to the dimensionality of the system.

Attraction Controls the Inversion of Order by Disorder in Buckled Colloidal Monolayers.

We show how including attraction in interparticle interactions reverses the effect of fluctuations in ordering of a prototypical artificial frustrated system. Buckled colloidal monolayers exhibit the same ground state as the Ising antiferromagnet on a deformable triangular lattice, but it is unclear if ordering in the two systems is driven by the same geometric mechanism. By a real-space expansion we find that, for buckled colloids, bent stripes constitute the stable phase, whereas in the Ising antiferromagnet straight stripes are favored. For generic pair potentials we show that attraction governs this selection mechanism, in a manner that is linked to local packing considerations. This supports the geometric origin of entropy in jammed sphere packings and provides a tool for designing self-assembled colloidal structures.

Hydrodynamics in kinetically constrained lattice-gas models.

Kinetically constrained models are lattice-gas models that are used for describing glassy systems. By construction, their equilibrium state is trivial and there are no equal-time correlations between the occupancy of different sites. We drive such models out of equilibrium by connecting them to two reservoirs of different densities, and we measure the response of the system to this perturbation. We find that under the proper coarse-graining, the behavior of these models may be expressed by a nonlinear diffusion equation, with a model- and density-dependent diffusion coefficient. We find a simple approximation for the diffusion coefficient, and we show that the relatively mild discrepancy between the approximation and our numerical results arises due to non-negligible correlations that appear as the system is driven out of equilibrium, even when the density gradient is infinitesimally small. Similar correlations appear when such kinetically constrained models are driven out of equilibrium by applying a uniform external force. We suggest that these correlations are the reason for the same discrepancy between the approximate diffusion coefficient and the numerical results for a broader group of models-nongradient lattice-gas models-for which kinetically constrained models are arguably the simplest example thereof.

Combinatorial design of textured mechanical metamaterials.

The structural complexity of metamaterials is limitless, but, in practice, most designs comprise periodic architectures that lead to materials with spatially homogeneous features. More advanced applications in soft robotics, prosthetics and wearable technology involve spatially textured mechanical functionality, which requires aperiodic architectures. However, a naive implementation of such structural complexity invariably leads to geometrical frustration (whereby local constraints cannot be satisfied everywhere), which prevents coherent operation and impedes functionality. Here we introduce a combinatorial strategy for the design of aperiodic, yet frustration-free, mechanical metamaterials that exhibit spatially textured functionalities. We implement this strategy using cubic building blocks-voxels-that deform anisotropically, a local stacking rule that allows cooperative shape changes by guaranteeing that deformed building blocks fit together as in a three-dimensional jigsaw puzzle, and three-dimensional printing. These aperiodic metamaterials exhibit long-range holographic order, whereby the two-dimensional pixelated surface texture dictates the three-dimensional interior voxel arrangement. They also act as programmable shape-shifters, morphing into spatially complex, but predictable and designable, shapes when uniaxially compressed. Finally, their mechanical response to compression by a textured surface reveals their ability to perform sensing and pattern analysis. Combinatorial design thus opens up a new avenue towards mechanical metamaterials with unusual order and machine-like functionalities.

Relation between structure of blocked clusters and relaxation dynamics in kinetically constrained models.

We investigate the relation between the cooperative length and relaxation time, represented, respectively, by the culling time and the persistence time, in the Fredrickson-Andersen, Kob-Andersen, and spiral kinetically constrained models. By mapping the dynamics to diffusion of defects, we find a relation between the persistence time, τ_{p}, which is the time until a particle moves for the first time, and the culling time, τ_{c}, which is the minimal number of particles that need to move before a specific particle can move, τ_{p}=τ_{c}^{γ}, where γ is model- and dimension-dependent. We also show that the persistence function in the Kob-Andersen and Fredrickson-Andersen models decays subexponentially in time, P(t)=exp[-(t/τ)^{ß}], but unlike previous works, we find that the exponent ß appears to decay to 0 as the particle density approaches 1.

Response of adherent cells to mechanical perturbations of the surrounding matrix.

We present a generic and unified theory to explain how cells respond to perturbations of their mechanical environment such as the presence of neighboring cells, slowly applied stretch, or gradients of matrix rigidity. Motivated by experiments, we calculate the local balance of forces that give rise to a tendency for the cell to locally move or reorient, with a focus on the contribution of feedback and homeostasis to cell contractility (manifested by a fixed displacement, strain or stress) that acts on the adhesions at the cell boundary. These forces can be either reinforced or diminished by elastic stresses due to mechanical perturbations of the matrix. Our model predicts these changes and how their balance with local protrusive forces that act on the cell's leading edge either increase or decrease the tendency of the cell to locally move (toward neighboring cells or rigidity gradients) or reorient (in the direction of slowly applied stretch or rigidity gradients).

Finite-density effects in the Fredrickson-Andersen and Kob-Andersen kinetically-constrained models.

We calculate the corrections to the thermodynamic limit of the critical density for jamming in the Kob-Andersen and Fredrickson-Andersen kinetically-constrained models, and find them to be finite-density corrections, and not finite-size corrections. We do this by introducing a new numerical algorithm, which requires negligible computer memory since contrary to alternative approaches, it generates at each point only the necessary data. The algorithm starts from a single unfrozen site and at each step randomly generates the neighbors of the unfrozen region and checks whether they are frozen or not. Our results correspond to systems of size greater than 10(7) × 10(7), much larger than any simulated before, and are consistent with the rigorous bounds on the asymptotic corrections. We also find that the average number of sites that seed a critical droplet is greater than 1.