Coherent potential approximation of random nearly isostatic kagome lattice

Zero modes activation to reconcile floppiness, rigidity, and multistability into an all-in-one class of reprogrammable metamaterials

Existing mechanical metamaterials are typically designed to either withstand loads as a stiff structure, shape morph as a floppy mechanism, or trap energy as a multistable matter, distinct behaviours that correspond to three primary classes of macroscopic solids. Their stiffness and stability are sealed permanently into their architecture, mostly remaining immutable post-fabrication due to the invariance of zero modes. Here, we introduce an all-in-one reprogrammable class of Kagome metamaterials that enable the in-situ reprogramming of zero modes to access the apparently conflicting properties of all classes. Through the selective activation of metahinges via self-contact, their architecture can be switched to acquire on-demand rigidity, floppiness, or global multistability, bridging the seemingly uncrossable gap between structures, mechanisms, and multistable matters. We showcase the versatile generalizations of the metahinge and remarkable reprogrammability of zero modes for a range of properties including stiffness, mechanical signal guiding, buckling modes, phonon spectra, and auxeticity, opening a plethora of opportunities for all-in-one materials and devices.

Transformable Superisostatic Crystals Self-Assembled from Segment Colloidal Rods

Transformable mechanical structures can switch between distinct mechanical states. Whether this kind of structure can be self-assembled from simple building blocks at microscale is a question to be answered. In this work, we propose a self-assembly strategy for these structures based on a nematic monolayer of segmented colloidal rods with lateral cutting. By using Monte Carlo simulation, we find that rods with different cutting degrees can self-assemble into different crystals characterized by bond coordination z that varies from 3 to 6. Among these, we identify a transformable superisostatic structure with pgg symmetry and redundant bonds (z = 5). We show that this structure can support either soft bulk modes or soft edge modes depending on its Poisson's ratio, which can be tuned from positive to negative through a uniform soft deformation. We also prove that the bulk soft modes are associated with states of self-stress along the direction of zero strain during uniform soft deformation. The self-assembled transformable structures may act as mechanical metamaterials with potential applications in micromechanical engineering.

Tunable assembly of host-guest colloidal crystals

Entropy compartmentalization provides new self-assembly routes to colloidal host-guest (HG) structures. Leveraging host particle shape to drive the assembly of HG structures has only recently been proposed and demonstrated. However, the extent to which the guest particles can dictate the structure of the porous network of host particles has not been explored. In this work, by modifying only the guest shape, we show athermal, binary mixtures of star-shaped host particles and convex polygon-shaped guest particles assemble as many as five distinct crystal structures, including rotator and discrete rotator guest crystals, two homoporous host crystals, and one heteroporous host crystal. Edge-to-edge alignment of neighboring stars results in the formation of three distinct pore motifs, whose preferential formation is controlled by the size and shape of the guest particles. Finally, we confirm, via free volume calculations, that assembly is driven by entropy compartmentalization, where the hosts and guests contribute differently to the free energy of the system; free volume calculations also explain differences in assembly based on guest shape. These results provide guest design rules for assembling colloidal HG structures, especially on surfaces and interfaces.

Dynamics of Self-Dual Kagome Metamaterials and the Emergence of Fragile Topology

Recent years have seen the discovery of systems featuring fragile topological states. These states of matter lack certain protection attributes typically associated with topology and are therefore characterized by weaker signatures that make them elusive to observe. Moreover, they are typically confined to special symmetry classes and, in general, rarely studied in the context of phononic media. In this Letter, we theoretically predict the emergence of fragile topological bands in the spectrum of a twisted kagome elastic lattice with threefold rotational symmetry, in the so-called self-dual configuration. A necessary requirement is that the lattice is a structural metamaterial, in which the role of the hinges is played by elastic finite-thickness ligaments. The interplay between the edge modes appearing in the band gaps bounding the fragile topological states is also responsible for the emergence of corner modes at selected corners of a finite hexagonal domain, which qualifies the lattice as a second-order topological insulator. We demonstrate our findings through a series of experiments via 3D scanning laser doppler vibrometry conducted on a physical prototype. The selected configuration stands out for its remarkable geometric simplicity and ease of physical implementation in the panorama of dynamical systems exhibiting fragile topology.

Universal scaling for disordered viscoelastic matter near the onset of rigidity

The onset of rigidity in interacting liquids, as they undergo a transition to a disordered solid, is associated with a rearrangement of the low-frequency vibrational spectrum. In this Letter, we derive scaling forms for the singular dynamical response of disordered viscoelastic networks near both jamming and rigidity percolation. Using effective-medium theory, we extract critical exponents, invariant scaling combinations, and analytical formulas for universal scaling functions near these transitions. Our scaling forms describe the behavior in space and time near the various onsets of rigidity, for rigid and floppy phases and the crossover region, including diverging length scales and timescales at the transitions.

Duality, Hidden Symmetry, and Dynamic Isomerism in 2D Hinge Structures

Recently, a new type of duality was reported in some deformable mechanical networks that exhibit Kramers-like degeneracy in phononic spectrum at the self-dual point. In this work, we clarify the origin of this duality and propose a design principle of 2D self-dual structures with arbitrary complexity. We find that this duality originates from the partial central inversion (PCI) symmetry of the hinge, which belongs to a more general end-fixed scaling transformation. This symmetry gives the structure an extra degree of freedom without modifying its dynamics. This results in dynamic isomers, i.e., dissimilar 2D mechanical structures, either periodic or aperiodic, having identical dynamic modes, based on which we demonstrate a new type of wave guide without reflection or loss. Moreover, the PCI symmetry allows us to design various 2D periodic isostatic networks with hinge duality. At last, by further studying a 2D nonmechanical magnonic system, we show that the duality and the associated hidden symmetry should exist in a broad range of Hamiltonian systems.

Self-Assembly of Isostatic Self-Dual Colloidal Crystals

Self-dual structures whose dual counterparts are themselves possess unique hidden symmetry, beyond the description of classical spatial symmetry groups. Here we propose a strategy based on a nematic monolayer of attractive half-cylindrical colloids to self-assemble these exotic structures. This system can be seen as a 2D system of semidisks. By using Monte Carlo simulations, we discover two isostatic self-dual crystals, i.e., an unreported crystal with pmg space-group symmetry and the twisted kagome crystal. For the pmg crystal approaching the critical point, we find the double degeneracy of the full phononic spectrum at the self-dual point and the merging of two tilted Weyl nodes into one critically tilted Dirac node. The latter is "accidentally" located on the high-symmetry line. The formation of this unconventional Dirac node is due to the emergence of the critical flatbands at the self-dual point, which are linear combinations of "finite-frequency" floppy modes. These modes can be understood as mechanically coupled self-dual rhombus chains vibrating in some unique uncoupled ways. Our work paves the way for designing and fabricating self-dual materials with exotic mechanical or phononic properties.

Multifunctional twisted kagome lattices: Tuning by pruning mechanical metamaterials

This article investigates phonons and elastic response in randomly diluted lattices constructed by combining (via the addition of next-nearest bonds) a twisted kagome lattice, with bulk modulus B=0 and shear modulus G>0, with either a generalized untwisted kagome lattice with B>0 and G>0 or with a honeycomb lattice with B>0 and G=0. These lattices exhibit jamming-like critical endpoints at which B, G, or both B and G jump discontinuously from zero while the remaining moduli (if any) begin to grow continuously from zero. Pairs of these jamming points are joined by lines of continuous rigidity percolation transitions at which both B and G begin to grow continuously from zero. The Poisson ratio and G/B can be continuously tuned throughout their physical range via random dilution in a manner analogous to "tuning by pruning" in random jammed lattices. These lattices can be produced with modern techniques, such as three-dimensional printing, for constructing metamaterials.

Signatures of Topological Phonons in Superisostatic Lattices

Soft topological surface phonons in idealized ball-and-spring lattices with coordination number z=2d in d dimensions become finite-frequency surface phonons in physically realizable superisostatic lattices with z>2d. We study these finite-frequency modes in model lattices with added next-nearest-neighbor springs or bending forces at nodes with an eye to signatures of the topological surface modes that are retained in the physical lattices. Our results apply to metamaterial lattices, prepared with modern printing techniques, that closely approach isostaticity.

Jamming as a Multicritical Point

The discontinuous jump in the bulk modulus B at the jamming transition is a consequence of the formation of a critical contact network of spheres that resists compression. We introduce lattice models with underlying undercoordinated compression-resistant spring lattices to which next-nearest-neighbor springs can be added. In these models, the jamming transition emerges as a kind of multicritical point terminating a line of rigidity-percolation transitions. Replacing the undercoordinated lattices with the critical network at jamming yields a faithful description of jamming and its relation to rigidity percolation.

Topological Edge Floppy Modes in Disordered Fiber Networks

Disordered fiber networks are ubiquitous in a broad range of natural (e.g., cytoskeleton) and manmade (e.g., aerogels) materials. In this Letter, we discuss the emergence of topological floppy edge modes in two-dimensional fiber networks as a result of deformation or active driving. It is known that a network of straight fibers exhibits bulk floppy modes which only bend the fibers without stretching them. We find that, interestingly, with a perturbation in geometry, these bulk modes evolve into edge modes. We introduce a topological index for these edge modes and discuss their implications in biology.

Elasticity of randomly diluted honeycomb and diamond lattices with bending forces

We use numerical simulations and an effective-medium theory to study the rigidity percolation transition of the honeycomb and diamond lattices when weak bond-bending forces are included. We use a rotationally invariant bond-bending potential, which, in contrast to the Keating potential, does not involve any stretching. As a result, the bulk modulus does not depend on the bending stiffness κ. We obtain scaling functions for the behavior of some elastic moduli in the limits of small ΔP = 1-P, and small δP = P-Pc, where P is an occupation probability of each bond, and Pc is the critical probability at which rigidity percolation occurs. We find good quantitative agreement between effective-medium theory and simulations for both lattices for P close to one.

Nonlinear elasticity of disordered fiber networks

Disordered biopolymer gels have striking mechanical properties including strong nonlinearities. In the case of athermal gels (such as collagen-I) the nonlinearity has long been associated with a crossover from a bending dominated to a stretching dominated regime of elasticity. The physics of this crossover is related to the existence of a central-force isostatic point and to the fact that for most gels the bending modulus is small. This crossover induces scaling behavior for the elastic moduli. In particular, for linear elasticity such a scaling law has been demonstrated [Broedersz et al. Nat. Phys., 2011 7, 983]. In this work we generalize the scaling to the nonlinear regime with a two-parameter scaling law involving three critical exponents. We test the scaling law numerically for two disordered lattice models, and find a good scaling collapse for the shear modulus in both the linear and nonlinear regimes. We compute all the critical exponents for the two lattice models and discuss the applicability of our results to real systems.

Finite-temperature mechanical instability in disordered lattices

Mechanical instability takes different forms in various ordered and disordered systems and little is known about how thermal fluctuations affect different classes of mechanical instabilities. We develop an analytic theory involving renormalization of rigidity and coherent potential approximation that can be used to understand finite-temperature mechanical stabilities in various disordered systems. We use this theory to study two disordered lattices: a randomly diluted triangular lattice and a randomly braced square lattice. These two lattices belong to two different universality classes as they approach mechanical instability at T=0. We show that thermal fluctuations stabilize both lattices. In particular, the triangular lattice displays a critical regime in which the shear modulus scales as G∼T(1/2), whereas the square lattice shows G∼T(2/3). We discuss generic scaling laws for finite-T mechanical instabilities and relate them to experimental systems.

Finite-temperature buckling of an extensible rod

Thermal fluctuations can play an important role in the buckling of elastic objects at small scales, such as polymers or nanotubes. In this paper, we study the finite-temperature buckling transition of an extensible rod by analyzing fluctuation corrections to the elasticity of the rod. We find that, in both two and three dimensions, thermal fluctuations delay the buckling transition, and near the transition, there is a critical regime in which fluctuations are prominent and make a contribution to the effective force that is of order √T. We verify our theoretical prediction of the phase diagram with Monte Carlo simulations.

Rigidity percolation by next-nearest-neighbor bonds on generic and regular isostatic lattices

We study rigidity percolation transitions in two-dimensional central-force isostatic lattices, including the square and the kagome lattices, as next-nearest-neighbor bonds ("braces") are randomly added to the system. In particular, we focus on the differences between regular lattices, which are perfectly periodic, and generic lattices with the same topology of bonds but whose sites are at random positions in space. We find that the regular square and kagome lattices exhibit a rigidity percolation transition when the number of braces is ∼LlnL, where L is the linear size of the lattice. This transition exhibits features of both first-order and second-order transitions: The whole lattice becomes rigid at the transition, and a diverging length scale also exists. In contrast, we find that the rigidity percolation transition in the generic lattices occur when the number of braces is very close to the number obtained from Maxwell's law for floppy modes, which is ∼L. The transition in generic lattices is a very sharp first-order-like transition, at which the addition of one brace connects all small rigid regions in the bulk of the lattice, leaving only floppy modes on the edge. We characterize these transitions using numerical simulations and develop analytic theories capturing each transition. Our results relate to other interesting problems, including jamming and bootstrap percolation.

Self-assembly of three-dimensional open structures using patchy colloidal particles

Open structures can display a number of unusual properties, including a negative Poisson's ratio, negative thermal expansion, and holographic elasticity, and have many interesting applications in engineering. However, it is a grand challenge to self-assemble open structures at the colloidal scale, where short-range interactions and low coordination number can leave them mechanically unstable. In this paper we discuss the self-assembly of three-dimensional open structures using triblock Janus particles, which have two large attractive patches that can form multiple bonds, separated by a band with purely hard-sphere repulsion. Such surface patterning leads to open structures that are stabilized by orientational entropy (in an order-by-disorder effect) and selected over close-packed structures by vibrational entropy. For different patch sizes the particles can form into either tetrahedral or octahedral structural motifs which then compose open lattices, including the pyrochlore, the hexagonal tetrastack and the perovskite lattices. Using an analytic theory, we examine the phase diagrams of these possible open and close-packed structures for triblock Janus particles and characterize the mechanical properties of these structures. Our theory leads to rational designs of particles for the self-assembly of three-dimensional colloidal structures that are possible using current experimental techniques.

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 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 interparticle interaction from the individual elastic constants and to derive an approximate estimate for the interparticle friction coefficient.

Shear shocks in fragile networks

A minimal model for studying the mechanical properties of amorphous solids is a disordered network of point masses connected by unbreakable springs. At a critical value of its mean connectivity, such a network becomes fragile: it undergoes a rigidity transition signaled by a vanishing shear modulus and transverse sound speed. We investigate analytically and numerically the linear and nonlinear visco-elastic response of these fragile solids by probing how shear fronts propagate through them. Our approach, which we tentatively label shear front rheology, provides an alternative route to standard oscillatory rheology. In the linear regime, we observe at late times a diffusive broadening of the fronts controlled by an effective shear viscosity that diverges at the critical point. No matter how small the microscopic coefficient of dissipation, strongly disordered networks behave as if they were overdamped because energy is irreversibly leaked into diverging nonaffine fluctuations. Close to the transition, the regime of linear response becomes vanishingly small: the tiniest shear strains generate strongly nonlinear shear shock waves qualitatively different from their compressional counterparts in granular media. The inherent nonlinearities trigger an energy cascade from low to high frequency components that keep the network away from attaining the quasi-static limit. This mechanism, reminiscent of acoustic turbulence, causes a superdiffusive broadening of the shock width.

Entropic effects in the self-assembly of open lattices from patchy particles

Open lattices are characterized by low-volume-fraction arrangements of building blocks, low coordination number, and open spaces between building blocks. The self-assembly of these lattices faces the challenge of mechanical instability due to their open structures. We theoretically investigate the stabilizing effects of entropy in the self-assembly of open lattices from patchy particles. A preliminary account of these findings and their comparison to experiment was presented recently [Mao, Chen, and Granick, Nat. Mater. 12, 217 (2013)]. We found that rotational entropy of patchy particles can provide mechanical stability to open lattices, whereas vibrational entropy of patchy particles can lower the free energy of open lattices and, thus, enables the selection of open lattices verses close-packed lattices which have the same potential energy. These effects open the door to significant simplifications of possible future designs of patchy particles for open-lattice self-assembly.

Elasticity of a filamentous kagome lattice

The diluted kagome lattice, in which bonds are randomly removed with probability 1-p, consists of straight lines that intersect at points with a maximum coordination number of 4. If lines are treated as semiflexible polymers and crossing points are treated as cross-links, this lattice provides a simple model for two-dimensional filamentous networks. Lattice-based effective-medium theories and numerical simulations for filaments modeled as elastic rods, with stretching modulus μ and bending modulus κ, are used to study the elasticity of this lattice as functions of p and κ. At p=1, elastic response is purely affine, and the macroscopic elastic modulus G is independent of κ. When κ=0, the lattice undergoes a first-order rigidity-percolation transition at p=1. When κ>0, G decreases continuously as p decreases below one, reaching zero at a continuous rigidity-percolation transition at p=p(b)≈0.605 that is the same for all nonzero values of κ. The effective-medium theories predict scaling forms for G, which exhibit crossover from bending-dominated response at small κ/μ to stretching-dominated response at large κ/μ near both p=1 and p(b), that match simulations with no adjustable parameters near p=1. The affine response as p→1 is identified with the approach to a state with sample-crossing straight filaments treated as elastic rods.

Entropy favours open colloidal lattices

Burgeoning experimental and simulation activity seeks to understand the existence of self-assembled colloidal structures that are not close-packed. Here we describe an analytical theory based on lattice dynamics and supported by experiments that reveals the fundamental role entropy can play in stabilizing open lattices. The entropy we consider is associated with the rotational and vibrational modes unique to colloids interacting through extended attractive patches. The theory makes predictions of the implied temperature, pressure and patch-size dependence of the phase diagram of open and close-packed structures. More generally, it provides guidance for the conditions at which targeted patchy colloidal assemblies in two and three dimensions are stable, thus overcoming the difficulty in exploring by experiment or simulation the full range of conceivable parameters.

Surface phonons, elastic response, and conformal invariance in twisted kagome lattices

Model lattices consisting of balls connected by central-force springs provide much of our understanding of mechanical response and phonon structure of real materials. Their stability depends critically on their coordination number z. d-dimensional lattices with z = 2d are at the threshold of mechanical stability and are isostatic. Lattices with z < 2d exhibit zero-frequency "floppy" modes that provide avenues for lattice collapse. The physics of systems as diverse as architectural structures, network glasses, randomly packed spheres, and biopolymer networks is strongly influenced by a nearby isostatic lattice. We explore elasticity and phonons of a special class of two-dimensional isostatic lattices constructed by distorting the kagome lattice. We show that the phonon structure of these lattices, characterized by vanishing bulk moduli and thus negative Poisson ratios (equivalently, auxetic elasticity), depends sensitively on boundary conditions and on the nature of the kagome distortions. We construct lattices that under free boundary conditions exhibit surface floppy modes only or a combination of both surface and bulk floppy modes; and we show that bulk floppy modes present under free boundary conditions are also present under periodic boundary conditions but that surface modes are not. In the long-wavelength limit, the elastic theory of all these lattices is a conformally invariant field theory with holographic properties (characteristics of the bulk are encoded on the sample boundary), and the surface waves are Rayleigh waves. We discuss our results in relation to recent work on jammed systems. Our results highlight the importance of network architecture in determining floppy-mode structure.

Nonlinear effective-medium theory of disordered spring networks

Disordered soft materials, such as fibrous networks in biological contexts, exhibit a nonlinear elastic response. We study such nonlinear behavior with a minimal model for networks on lattice geometries with simple Hookian elements with disordered spring constant. By developing a mean-field approach to calculate the differential elastic bulk modulus for the macroscopic network response of such networks under large isotropic deformations, we provide insight into the origins of the strain stiffening and softening behavior of these systems. We find that the nonlinear mechanics depends only weakly on the lattice geometry and is governed by the average network connectivity. In particular, the nonlinear response is controlled by the isostatic connectivity, which depends strongly on the applied strain. Our predictions for the strain dependence of the isostatic point as well as the strain-dependent differential bulk modulus agree well with numerical results in both two and three dimensions. In addition, by using a mapping between the disordered network and a regular network with random forces, we calculate the nonaffine fluctuations of the deformation field and compare them to the numerical results. Finally, we discuss the limitations and implications of the developed theory.