Abstract
Mechanical metamaterials are artificial structures with unusual properties, such as negative Poisson ratio, bistability or tunable vibrational properties, that originate in the geometry of their unit cell1,2,3,4,5. Often at the heart of such unusual behaviour is a soft mode: a motion that does not significantly stretch or compress the links between constituent elements. When activated by motors or external fields, soft modes become the building blocks of robots and smart materials. Here, we demonstrate the existence of topological soft modes that can be positioned at desired locations in a metamaterial while being robust against a wide range of structural deformations or changes in material parameters6,7,8,9,10. These protected modes, localized at dislocations in deformed kagome and square lattices, are the mechanical analogue of topological states bound to defects in electronic systems11,12,13,14. We create physical realizations of the topological modes in prototypes of kagome lattices built out of rigid triangular plates. We show mathematically that they originate from the interplay between two Berry phases: the Burgers vector of the dislocation and the topological polarization of the lattice. Our work paves the way towards engineering topologically protected nanomechanical structures for molecular robotics or information storage and read-out.
Similar content being viewed by others
Main
Central to our approach is a simple insight: mechanical structures on length scales ranging from the molecular to the architectural can often be viewed as networks of nodes connected by links15. Whether the linking components are chemical bonds or metal beams, mechanical stability depends crucially on the number of constraints relative to the degrees of freedom. When the degrees of freedom exceed the constraints, the structure exhibits excess zero (potential) energy modes. Conversely, when the constraints exceed the degrees of freedom, there are excess states of self-stressâbalanced combinations of tensions and compressions of the links, with no resultant force on the nodes. The generalized Maxwell relation16 stipulates that the index ν, given by the difference between the number of zero modes, nm, and the number of states of self-stress, nss, is equal to the difference between the number of degrees of freedom, Ndf, and the number of constraints, Nc,
A trivial way to position a zero-energy mode in the interior of a generic rigid lattice is to remove some bonds, locally reducing the number of constraints. Consider, instead, a network that satisfies everywhere the local isostatic condition Ndf = Nc (which precludes bond removal). In this case, zero modes can be present only in conjunction with an equal number of states of self-stress, invisible partners from the perspective of motion. Isostaticity by itself, however, does not dictate how the modes are distributed spatially. Kane and Lubensky6 recently introduced a special class of isostatic lattices that possesses a further feature called topological polarization. Much as an electrically polarized material can localize positive and negative charges at opposite boundaries, a topologically polarized lattice can harbour zero modes or states of self-stress at sample edges (whose outward normal is respectively aligned or anti-aligned with the polarization). The edge mode distribution is biased even though both boundaries are indistinguishable on the basis of local constraint counting. Furthermore, this bias is insensitive to local variations in bond lengths, spring constants, or node masses, provided no bonds are cut and the lattice remains rigid in the bulk6.
In this Letter, we harness the topological polarization to place zero modes in the interior of an isostatic lattice where topological defects called dislocations are positioned. Dislocations are termination points of incomplete lattice rows that have an edge-like character. They are characterized by a topological charge called the Burgers vector, b, which measures the deficit in any circuit surrounding the dislocation, see Fig. 1a, c. A dislocation is composed of a dipole of under-coordinated (green) and over-coordinated (orange) points (Fig. 1a) or plaquettes (Fig. 1bâd), whose orientation is obtained on rotating b by Ï/2. The dipole moment, d, a vector connecting the under-coordinated point/plaquette to the over-coordinated one, points outward from the added strip of material that terminates at the dislocation and has a magnitude equal to the strip width. Therefore, d quantifies the orientation and size of the effective âedgeâ created by the dislocation.
To localize topological modes at these effective edges, we need to incorporate the dislocations into polarized lattices without modifying the local constraint count (that is, trivial zero modes must be excluded). We demonstrate this construction for two polarized lattices: a deformed kagome lattice introduced in ref. 6 and a deformed square lattice. As shown in Fig. 1, both lattices are obtained by decorating a 2D crystal lattice (regular hexagonal and square, with primitive vectors {a1, a2} indicated in Fig. 1a, c, respectively) with a multi-atom basis at each unit cell. In the absence of defects, the unit cell determines the topological polarization PT of the bulk lattice. In the Supplementary Information we show that PT = a1 â a2 for the deformed square lattice; it has been shown in ref. 6 that PT = a1 for the deformed kagome lattice. Dislocations in the undecorated lattice carry over to the polarized lattice and, when appropriately chosen, produce a lattice that is four-coordinated everywhere (Fig. 1b, d). See the Supplementary Information for more details of the construction.
Having constructed polarized lattices with dislocations, we numerically compute their vibrational spectrum by treating each bond as a harmonic spring. Results are shown in Fig. 2a, b for deformed kagome and square networks, respectively. We use periodic boundary conditions, which preserves isostaticity everywhere but requires a net Burgers vector of zero. As a result, dislocations appear in pairs with equal and opposite dipole moments. In both networks, the dipole moments of the dislocation on the left (dL = a1 â a2) and on the right (dR = â dL) are aligned, respectively, with and against the lattice polarization. The left dislocation has an associated soft mode in both cases, labelled by arrows, whose energy decreases with system size. The opposite dislocation is associated with an approximate state of self-stress, labelled by coloured and thickened bonds following ref. 6. See Supplementary Information for computational details. These observations are consistent with the intuitive interpretation of dislocations as edges oriented by their dipole moment.
To assess whether these modes can be observed in metamaterials with realistic bonds, joints and boundary conditions, we have built prototypes of the deformed kagome lattice, composed of rigid triangles laser-cut out of 3 mm thick PMMA sheets. The corner of the triangles are connected by plastic bolts that act as hinges. The boundary points are pinned to a flat base by screws, but can pivot freely. The design ensures that each internal vertex has as many constraints as degrees of freedom, satisfying the local isostatic condition away from the boundary. Figure 2c shows such a prototype mimicking the dislocation configuration of the computer model from Fig. 2a. Theoretically, the boundary pinning and the use of rigid triangles push the phonon gap to infinity, so that only zero modes can be observed. In practice, the prototype has some compliance and mechanical play at the pivots. Nonetheless, it is rigid in the bulk, as can be verified by unsuccessfully trying to move the (white) triangles far from the dislocations, see Supplementary Movie 1.
Despite the differences between the two systems, the soft mode observed in the simulated harmonic network survives in the real-world prototype as a collective motion of points near the left dislocation. The motion is easily activated by pushing the hinge joints of the triangles that make up the dislocation (inset to Fig. 2c and Supplementary Movie 2). The motion is not a strict zero mode because it interacts with the pinned boundary of the finite system, but the structural compliance is sufficient for the remnant soft motion to be observed. In contrast, the dislocation on the right does not admit displacements in its vicinity and remains rigid (Supplementary Movie 3), consistent with the simulations.
To quantify the number and type of modes associated with a dislocation, an electrostatic analogy proves useful. Once the connectivity of a locally isostatic lattice is fixed, the index ν in equation (1) can be viewed as a topological charge, invariant under smooth deformations of the lattice. Just as the Gauss law yields the net charge enclosed in a region from the flux of the electric polarization through its boundary, the net value of ν in an arbitrary portion of an isostatic lattice is given by the flux of the topological polarization through its boundary6. On evaluating the flux of PT on a contour encircling an isolated dislocation, we obtain
where Vcell is the unit cell area. In the Supplementary Information, we present a detailed derivation of equation (2) that accounts for the elastic strains around the dislocation. Here, we simply comment on its physical interpretation.
The topologically protected modes arise from a delicate interplay between a Berry phase associated with cycles in the Brillouin zone, embedded in PT (ref. 6), and the Berry phase of a topological defect in real space, represented by its Burgers vector (or dipole d). A similar interplay dictates the existence of localized electronic modes at dislocations in conventional topological insulators12,13. Equation (2) gives ν = +1â(â1) for the left (right) dislocation in the deformed kagome lattice, and ν = +2â(â2) in the deformed square lattice. The sign of ν distinguishes zero modes (+) from states of self-stress (â), while its magnitude gives their numbers. For instance, we correctly predict that the square lattice of Fig. 2b admits two soft modes localized to the left dislocation (as verified in Supplementary Fig. 1).
The soft modes investigated here have unusual localization properties, as shown in Fig. 3. The mode amplitude falls off exponentially along most rays originating at the core of the left dislocation, but the decay length depends on the direction of the ray relative to the underlying lattice. There are two special directions in each lattice (one of which is highlighted by red circles in both Fig. 3a, b) along which the localization is weak. In all other directions, the mode decays over much shorter length scales, of the order of a few lattice constants (green circles and symbols in Fig. 3). In the Supplementary Information, we show that the weak localization directions track lines in momentum space along which the acoustic modes vary quadratically, rather than linearly, with the momentum. The localization of the approximate state of self-stress behaves similarly, with weak directions that are the opposite of those for the soft mode.
The topologically protected modes we have identified could have applications across a wide range of systems and length scales. At macroscopic scales, isostatic origami structures exist whose deformations are restricted to rotations of hinged triangles, much as in the kagome lattice17. At the microscale, dislocations could be used for robust information storage, with a bit encoded by the presence (+) or absence (â) of a topological soft mode, in turn controlled by the orientation of the Burgers vector. Such protected bits could be hard-wired into microscopic âpunch cardsâ that could be read out mechanically by probing the region around dislocations for soft motions (or lack thereof). We also envisage molecular robots and smart metamaterials that could exploit the protected modes as activated mechanisms.
References
Babaee, S., Shim, J., Weaver, J. C., Patel, N. & Bertoldi, K. 3D soft metamaterials with negative Poissonâs ratio. Adv. Mater. 25, 5044â5049 (2013).
Schenk, M. & Guest, S. D. Geometry of Miura-folded metamaterials. Proc. Natl Acad. Sci. USA 110, 3276â3281 (2013).
Wei, Z. Y., Guo, Z. V., Dudte, L., Liang, H. Y. & Mahadevan, L. Geometric mechanics of periodic pleated origami. Phys. Rev. Lett. 110, 215501 (2013).
Sun, K., Souslov, A., Mao, X. & Lubensky, T. C. Surface phonons, elastic response, and conformal invariance in twisted kagome lattices. Proc. Natl Acad. Sci. USA 109, 12369â12374 (2012).
Shan, S. et al. Harnessing multiple folding mechanisms in soft periodic structures for tunable control of elastic waves. Adv. Funct. Mater. 24, 4935â4942 (2014).
Kane, C. L. & Lubensky, T. C. Topological boundary modes in isostatic lattices. Nature Phys. 10, 39â45 (2014).
Prodan, E. & Prodan, C. Topological phonon modes and their role in dynamic instability of microtubules. Phys. Rev. Lett. 103, 248101 (2009).
Vitelli, V. Topological soft matter: Kagome lattices with a twist. Proc. Natl Acad. Sci. USA 109, 12266â12267 (2012).
Chen, B. G., Upadhyaya, N. & Vitelli, V. Nonlinear conduction via solitons in a topological mechanical insulator. Proc. Natl Acad. Sci. USA 111, 13004â13009 (2014).
Vitelli, V., Upadhyaya, N. & Chen, B. G. Topological mechanisms as classical spinor fields. Preprint at http://arXiv.org/abs/1407.2890 (2014).
Stern, A. Anyons and the quantum Hall effect: A pedagogical review. Ann. Phys. 323, 204â249 (2008).
Ran, Y., Zhang, Y. & Vishwanath, A. One-dimensional topologically protected modes in topological insulators with lattice dislocations. Nature Phys. 5, 298â303 (2009).
Teo, J. C. Y. & Kane, C. L. Topological defects and gapless modes in insulators and superconductors. Phys. Rev. B 82, 115120 (2010).
JuriÄiÄ, V., Mesaros, A., Slager, R-J. & Zaanen, J. Universal probes of two-dimensional topological insulators: Dislocation and Ï flux. Phys. Rev. Lett. 108, 106403 (2012).
Maxwell, J. C. On the calculation of the equilibrium and stiffness of frames. Phil. Mag. 27, 294â299 (1864).
Calladine, C. Buckminster Fullerâs âTensegrityâ structures and Clerk Maxwellâs rules for the construction of stiff frames. Int. J. Solids Struct. 14, 161â172 (1978).
Tachi, T. Designing freeform origami tessellations by generalizing Reschâs patterns. J. Mech. Design 135, 111006 (2013).
Acknowledgements
We thank J. C. Y. Teo, R-J. Slager and A. Turner for helpful discussions, and the Leiden University Fine Mechanics Department for technical support. This work was funded by FOM and by the D-ITP consortium, a program of the Netherlands Organisation for Scientific Research (NWO), funded by the Dutch Ministry of Education, Culture and Science (OCW).
Author information
Authors and Affiliations
Contributions
J.P., B.G.C. and V.V. designed the research; J.P. and B.G.C. conducted the research and interpreted the results; V.V. supervised the research and interpreted the results; J.P., B.G.C. and V.V. prepared the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Information (PDF 1426 kb)
Supplementary Movie
Supplementary Movie 1 (MP4 1762 kb)
Supplementary Movie
Supplementary Movie 2 (MP4 1719 kb)
Supplementary Movie
Supplementary Movie 3 (MP4 1030 kb)
Rights and permissions
About this article
Cite this article
Paulose, J., Chen, Bg. & Vitelli, V. Topological modes bound to dislocations in mechanical metamaterials. Nature Phys 11, 153â156 (2015). https://doi.org/10.1038/nphys3185
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys3185
This article is cited by
-
Breakdown of conventional winding number calculation in one-dimensional lattices with interactions beyond nearest neighbors
Communications Physics (2023)
-
Topological phenomena at defects in acoustic, photonic and solid-state lattices
Nature Reviews Physics (2023)
-
Strongly nonlinear topological phases of cascaded topoelectrical circuits
Frontiers of Physics (2023)
-
Conformal elasticity of mechanism-based metamaterials
Nature Communications (2022)
-
Observation of topological action potentials in engineered tissues
Nature Physics (2022)