# Topological mechanics of gyroscopic metamaterials

See allHide authors and affiliations

Edited by David A. Weitz, Harvard University, Cambridge, MA, and approved October 13, 2015 (received for review April 17, 2015)

## Significance

We have built a new type of mechanical metamaterial: a “gyroscopic metamaterial” composed of rapidly spinning objects that are coupled to each other. At the edges of these materials, we find sound waves that are topologically protected (i.e. they cannot be scattered backward or into the bulk). These waves, which propagate in one direction only, are directly analogous to edge currents in quantum Hall systems. Through a mathematical model, we interpret the robustness of these edge waves in light of the subtle topological character of the bulk material. Crucially, these edge motions can be controlled by distorting the metamaterial lattice, opening new avenues for the control of sound in matter.

## Abstract

Topological mechanical metamaterials are artificial structures whose unusual properties are protected very much like their electronic and optical counterparts. Here, we present an experimental and theoretical study of an active metamaterial—composed of coupled gyroscopes on a lattice—that breaks time-reversal symmetry. The vibrational spectrum displays a sonic gap populated by topologically protected edge modes that propagate in only one direction and are unaffected by disorder. We present a mathematical model that explains how the edge mode chirality can be switched via controlled distortions of the underlying lattice. This effect allows the direction of the edge current to be determined on demand. We demonstrate this functionality in experiment and envision applications of these edge modes to the design of one-way acoustic waveguides.

A vast range of mechanical structures, including bridges, covalent glasses, and conventional metamaterials, can be ultimately modeled as networks of masses connected by springs (1⇓⇓⇓⇓–6). Recent studies have revealed that despite its apparent simplicity, this minimal setup is sufficient to construct topologically protected mechanical states (7⇓⇓⇓–11) that mimic the properties of their quantum analogs (12). This follows from the fact that, irrespective of its classic or quantum nature, a periodic material with a gapped spectrum of excitations can display topological behavior as a result of the nontrivial topology of its band structure (13⇓⇓⇓⇓⇓⇓⇓–21).

All such mechanical systems, however, are invariant under time reversal because their dynamics are governed by Newton’s second law, which, unlike the Schrödinger equation, is second order in time. If time-reversal symmetry is broken, as in recently suggested acoustic structures containing circulating fluids (16), theoretical work (13) has suggested that phononic chiral topological edge states that act as unidirectional waveguides resistant to scattering off impurities could be supported. In this paper, we show that by creating a coupled system of gyroscopes, a “gyroscopic metamaterial,” we can produce an effective material with intrinsic time-reversal symmetry breaking. As a result, our gyroscopic metamaterials support topological mechanical modes analogous to quantum Hall systems, which have robust chiral edge states (22⇓–24). We demonstrate these effects by building a real system of gyroscopes coupled in a honeycomb lattice. Our experiments show long-lived, unidirectional transport along the edge, even in the presence of significant defects. Moreover, our theoretical analysis indicates that direction of edge propagation is controlled both by the gyroscope spin and the geometry of the underlying lattice. As a result, deforming the lattice of gyroscopes allows one to control the edge mode direction, offering unique opportunities for engineering novel materials.

Much of the counterintuitive behavior of rapidly spinning objects originates from their large angular momentum, which endows the axis of spin with a resistance to change. If we fix one end of a gyroscope and apply a force, *ω* is the gyroscope angular frequency and *I* is its rotational inertia. The behavior of a gyroscope differs from that of a simple mass in two important ways: (*i*) It moves perpendicular to applied forces and (*ii*) its response is first order in time. The canonical example of this unusual behavior is precession: A spinning top does not simply fall over, but rather its free end orbits around the contact point (precesses) with a constant period,

What happens if we replace the masses in a conventional mechanical metamaterial with gyroscopes? A first glimpse is provided by a normal mode analysis of honeycomb lattices composed of mass–spring and gyroscope–spring networks. The density of states of these two systems (Fig. 1*A*) shows qualitatively similar features: Each is characterized by two bands, a lower “acoustic” band (where neighboring sites move in phase) and an upper “optical” band (where neighboring sites move out of phase). The connections between these two bands, however, show key differences: In the mass–spring system the two bands touch at a Dirac point, whereas in the gyroscopic system a gap opens up between the bands. Crucially, this gap is not empty, but populated by nearly equally spaced modes and the number of these edge modes scales with the length of the edge. Examination of these gap modes reveals them to be confined to the edge and to be chiral: The phases always rotate in the same direction as one moves around the lattice (Fig. 1*B*). As we show below, these edge modes are topologically robust and can therefore serve as unidirectional waveguides.

It is far from obvious that in a real system the phonon spectrum would be resistant to the presence of both disorder (lattice imperfections, gyroscopic nonuniformity, etc.) and mixed-order dynamics (e.g., nutation). However, an appealing feature of topological states is that they are often resistant to disorder, suggesting that they may be useful for acoustic applications and observable under a wide range of experimental conditions.

To explore the relevance of these effects, we have assembled a prototype system of 54 interacting gyroscopes on a honeycomb lattice (Fig. 2*A* and Fig. S2). Our gyroscopes consist of small dc motors spinning cylindrical masses at ∼300 Hz (with *B*), producing an individual precession frequency of

The magnetically coupled system is conceptually equivalent to the system of gyrosocopes connected by springs discussed earlier; the linearized magnetic coupling differs, however, in detail from the coupling given by springs because the equilibrium results from the cancellation of opposing forces instead of the absence of forces. As detailed in *Supporting Information*, this results in a mode spectrum that is shifted to lower frequencies (Fig. 2*C*), compared with a spring-coupled gyroscope system (Fig. 1*A*). However, the topological character of the band structure is not affected and acoustic and optical bands are still apparent with chiral edge modes in between.

To test the mechanical response of the gyroscopic metamaterial, we excite it with periodic bursts of air through a small nozzle and follow the resulting disturbance. We probe the normal modes by weakly exciting a single gyroscope at a fixed frequency for many (>100) periods and recording the resulting motion of the network. All excitations were kept to small amplitudes (<10% of the lattice spacing) to avoid the nonlinearities associated with coupling the gyroscopes magnetically (Figs. S3–S5).

The effect of disorder inherent to our experiment (e.g., variation in motor speed and gyroscope pivot position) can be clearly seen in the comparison between the structure of bulk modes as shown in Fig. 2 *D* and *F* for the idealized (*Left*) and experimental system (*Right*). There is little overall agreement between calculated and measured modes, though acoustic modes show approximate in-phase oscillation of adjacent gyroscopes and optical modes show approximate out-of-phase oscillation of adjacent gyroscopes. This is characteristic of the effects of disorder (25), which produces the same effect in numerically evaluated modes with comparable disorder (Fig. S6).

Despite these experimental imperfections, exciting a mode in the gap between the acoustic and optical bands produces clean excitations along the edge (Fig. 2*E*). The orientations and relative orbit sizes of these modes closely match the modes numerically computed for an idealized model. A comparison between calculated and measured edge mode frequencies is shown in Fig. S5. The robustness of these modes against disorder is characteristic of their topological character.

To demonstrate that our experimental metamaterial functions as a unidirectional waveguide, we excite a single edge gyroscope for five periods at a frequency in the gap. As shown in Fig. 3*A* and Movie S1, the resulting excitation propagates in only one direction around the edge of the lattice. The motion of this wave packet around the edge is persistent, circumnavigating the boundary several times. As expected, short excitations at a frequency not in the band gap do not produce a similar robust edge excitation (Movie S2). Crucially, this indicates that the chiral edge modes are topologically protected from coupling to the bulk modes, functioning as an efficient one-directional waveguide.

We further demonstrate the robustness of these edge modes by intentionally introducing disorder in the lattice, for example by removing three gyroscopes. As shown in Movie S3 and Fig. 3*B*, even this significant disturbance does not destroy the chiral edge modes. An excitation on the edge is seen to move around this disturbance—in the same direction as before—and emerge undisturbed on the other side. Remarkably, the excitation traverses the defect region without scattering backward or into the bulk. As before, the resilience of the edge modes suggests these edge states are topological in character.

To analyze the origin of these effects, we return to an ideal coupled gyroscope model. For simplicity, we represent the displacement of the tip of the gyroscope from equilibrium as **1** is *i*, arises from the cross-product. Accordingly, the linearized equation of motion for each site in the gyroscopic metamaterial is*p* is the site label, *q* the neighboring sites, *p*, *k* is the spring constant.

Symmetries often play a fundamental role in characterizing a system’s topological behavior; in the case of the gyroscopic materials, broken time-reversal symmetry is a natural starting point. We note that the linearized equation of motion bears remarkable similarity to the Schrödinger equation for the wavefunction of an electron in a tight-binding model. Thus, by analogy, we may analyze the breaking of temporal symmetry using the “time-reversal” operation in quantum mechanics: *ψ* mirrors their displacement in the *y* axis; applying the complete time-reversal operation to a single gyroscope leaves the equation of motion unchanged. Similarly, for a network of gyroscopes Eq. **2** is invariant under this operation only if the coefficient

The geometric origin of the time-reversal symmetry breaking can also be seen in the case of gyroscopes connected by springs, by considering the energy of two connected gyroscopes. In the linearized limit, the stretching/compression of the spring is given by *y* axis (*A*). When considering an entire lattice, we see that for arbitrary displacements the bond energy will be conserved under time reversal if (and only if) we are able to choose a global mirror axis to which all bonds are either perpendicular or parallel. As a result, time-reversal invariance is only guaranteed for lattices composed of square or rectangular building blocks.

It is instructive to note that, in the limit that gyroscopes are coupled by weak springs, *Supporting Information* and Figs. S7 and S8 for details) (26). In the Haldane model, time-reversal symmetry is broken by a staggered magnetic field. This field can be varied, resulting in a change in the topological character of the modes as quantified by the Chern number (12, 27, 28). Accordingly, depending on the strength of the field and asymmetry between the two sites in the unit cell, the Chern number of the bottom band is

In gyroscopic metamaterials, the analog to changing the magnetic field is to geometrically distort the lattice. In either case, the relevant operation produces a phase shift in the hopping between neighboring sites; in the gyroscope system this phase shift is determined by the bond angles,

In a honeycomb lattice, it is possible to distort the constituent hexagons without changing the bond length (Fig. 4 *B*–*D*), allowing us to change the gyroscopic phase between neighboring sites without changing the network connectivity. As predicted by the time-reversal analysis above, the band-gap and chiral edge modes disappear when the bonds fall on a rectangular grid (in which case *C*); this is topologically equivalent to the dispersion relationship of a honeycomb network of springs and masses. Continuing to distort the lattice past this point restores the band gap, but the edge modes now have opposite chirality, as reflected in an inversion of the bands and hence of the Chern number;

We have presented an experimental proof of concept and theoretical analysis of a topologically protected unidirectional waveguide in a real mechanical metamaterial. The origin of our topological edge modes is due to time-reversal symmetry breaking; our analysis indicates this arises from the combination of the chiral nature of the gyroscopes and the geometry of the underlying lattice. Because the direction of the edge modes can be changed discontinuously with geometric distortions, in principle small displacements should be capable of inverting the edge mode direction. This mechanism may have practical applications for creating direction-tunable materials, but it also suggests interesting nonlinear effects should occur in the regime near the mechanically induced topological phase transition.

The prototypical gyroscopic solids we have developed here are examples of active metamaterials: Their design relies on the presence of internal motors that keep each gyroscope in a fast spinning state. An open challenge is to construct scalable gyroscopic metamaterials using nanofabrication techniques (e.g., microelectromechanical systems) or active molecules that convert chemical energy into rotation very much like the motors powering each gyroscope (30, 31). Such an implementation would pave the way toward realizing materials that support, at a microscopic scale, robust topological acoustic modes.

### Note Added in Proof.

In the concluding stages of the present work, we became aware of a parallel independent effort in which a class of topological gyroscopic metamaterials was theoretically analyzed (32).

## Equation of Motion

If a gyroscope is spinning with a large constant angular frequency, *ω*, around its principal axis, *z*-component of

When a gyroscope is subject to a restoring force (e.g., gravity),

## Linearized Equation of Motion

To find the linearized equation of motion, we consider first an interaction of two gyroscopes *p* and *q*, with *p* being at the origin and *q* a distance *a* from *p* at an angle

The effective spring constant for the interaction is given by the gradient of the force. In general, these gradients may be asymmetric in the components that are parallel and perpendicular to vector from the bond, which is defined as the line connecting points *p* and *q*. In the case of a linear spring with constant

As an illustration, we now consider the linearized equation of motion for gyroscopes *p* and *q* interacting via this linear spring. In this case, we need only to find the displacement parallel to the bond (*p* to *q* as *x* axis, taking the real part of **S1** has resulted in the imaginary coefficient on the left hand side of Eq. **S4**.

For a general radial interaction force, *a* is the separation between the two lattice sites.

In our experiment, the gyroscopes are coupled through small magnets. The force can be approximated by treating each gyroscope as a magnetic dipole with strength *M*; this produces an *a*) to pendulum length (

In a honeycomb lattice, the symmetry of the lattice allows for the equations of motion to be simplified to

The equation of motion for the magnets, Eq. **S7**, is nearly equivalent to the simple linear spring case, excepting the slight asymmetry between the *ψ* and **S4**.) The asymmetry of these terms has an interesting effect when viewed from the perspective of time-reversal symmetry breaking: The

## Experimental Details

The gyroscopic metamaterial was constructed from 54 hanging gyroscopes, as shown in Fig. S2. Each individual gyroscope consisted of a 3D printed cylindrical rotor (radius 10.4 mm) heat fitted to the shaft of a small (20 mm long) brushed dc pager motor. An N40 neodymium magnet (9-mm diameter, 4-mm thickness) was embedded with its moment aligned with the *z* axis in each rotor. The spinning mass of each gyroscope rotor was 6.1 g. The dipole moment of the magnets was determined from direct measurement of the force between attracting magnets (with dipole moments aligned) to be ^{2}. Each spinning mass was enclosed in a housing that was glued to the motor shaft and suspended by a weak spring. When the end of the gyroscope is displaced, the bending of the spring gives an effective pendulum length of

The spacing between gyroscopes in the experimental lattice was 30.5 mm, corresponding to an effective spring constant that is roughly equivalent to the magnitude of the effective gravitational pinning spring,

Fixed magnets (N52, 10-mm o.d., 5.5-mm i.d., 3 mm thick) were placed at the perimeter of the system in the position of each edge gyroscope’s missing neighbor. These fixed magnets were approximately two-thirds the strength of the inner magnets. We note that this detail is not crucial to the topological nature of the system. The strength of the boundary magnets does not change the frequencies of the gap edges (which depend on the bulk properties), and only slightly affects frequency distribution of modes within the gap.

We obtained the experimental normal modes for the frequency range

## Linearity of the Magnetically Coupled System

It has been shown experimentally that lattices of interacting magnets can exhibit nonlinear behavior when excited at sufficiently large amplitudes (33). To ensure that our experiments were performed in a regime where linear analysis can be applied, we tested for the presence of nonlinear effects with three experiments that are summarized here and listed in detail below. (*i*) We characterized the shift in the normal mode frequencies of a pair of interacting gyroscopes as a function of the excitation amplitude. (*ii*) We verified that the frequency of an edge mode in our system of 54 interacting gyroscopes does not shift significantly as the amplitude of excitation was changed. (*iii*) We measured speed of the wave packet shown in Movie S1 and verified that the change in amplitude does not significantly affect the average speed of the wave packet. Our measurements of the effects of nonlinearities (detailed below) cover a range of excitation amplitude beyond the maximum excitation amplitude in the experiments of the main text.

*i*) Fig. S3 shows the measured normal mode frequencies of a pair of interacting gyroscopes, each with two neighboring boundary magnets in a honeycomb configuration (as described in*Experimental Details*). Two normal mode frequencies are expected: a lower frequency mode in which the gyroscopes precess in phase (red dots in Fig. S3) and a higher frequency mode in which the gyroscopes precess out of phase (green dots in Fig. S3). The frequencies of both modes were tracked as the amplitude of an excitation damped during an interaction. The nonlinearity of the interaction is apparent because the mode frequencies are not constant with amplitude. However, at the upper limit of excitation amplitude of the experiments presented in the main text (3 mm), the deviations of the frequency are not appreciable. At this amplitude, a frequency shift of only less than 2% was observed for either of the modes. We note that the frequency shift reaches a maximum of 5% for the antisymmetric mode when the excitation amplitude reaches 4 mm.*ii*) The frequency shift of a gap mode was measured at varying excitation amplitudes. For each test a region of 0.06 Hz was tested near a known peak, and each frequency was excited for 80 s before recording and subsequently recorded for 100 s, as described previously. We determined that for changes in excitation amplitude ranging from 1.2 to 4.0 mm, the shift in frequency was less than 0.02 Hz. An example of the effect is shown in Fig. S5*B*. The frequencies for the fourth mode number are plotted for three amplitudes in the range indicated above.*iii*) We tracked the position of the edge mode wave packet shown in Fig. S4 and Movie S1 as the amplitude of oscillation decayed. The position of the wave packet was found in each frame using a center-of-mass method, in which each lattice point was weighted by the squared displacement of its gyroscope. As shown in Fig. S4, the average speed of the wave packet is unaffected by the gradual decay in amplitude.

Furthermore, we do not observe mode mixing at the amplitudes of the data presented in this work.

## Normal Mode Analysis

For second-order mass and spring systems with *n* lattice sites in two dimensions, we find the normal modes by considering the system of equations*x* and *y* displacements of each mass in the network and **K** is a coupling matrix between sites. The normal modes are found by finding the eigenvalues, **S5**. As with the spring and mass system, these equations of motion can be expressed as a matrix:*Experimental Details*. The energy minimization gave the equilibrium points that were used in the finite-system normal mode calculation. The spring constant between gyroscopes was given by *r* was the separation between interacting gyroscopes. The mode frequencies were calculated by assuming coupling of individual gyroscopes with all gyroscopes in the lattice as well as with only nearest-neighbor coupling.

A comparison between the expected gap modes from the model and the measured edge modes in the experiment is shown in Fig. S5. The shaded regions indicate the range of possible gap mode frequencies calculated with the measured experimental values ^{2} and

We note that disorder (numerical and experimental) can greatly alter mode profiles. We find that experimentally observed mode profiles are qualitatively similar to numerically calculated mode profiles in a system with random disorder, as shown in Fig. S6.

## Time Domain Simulation

We simulate a 2D gyromaterial in the time domain only for the case of spring interactions with free boundary conditions. We numerically integrate Eq. **S1** considering a spring interaction and gravity. We consider only the *x* and *y* displacements and integrate using a fourth order Runge–Kutta method.

## Band Structure and Chern Number Calculation

We find the band structure for both the spring and magnetically coupled systems on a honeycomb lattice using the linearized equations of motion, and assuming the solutions*a* and *b* refer to the two sites in each unit cell.

The resulting equations can be expressed as a *n* is the number of lattice sites per unit cell.) The four eigenvalues of this matrix give the values of the four dispersion bands at a particular value of

The Chern number of each band is given by an integral of the Berry curvature

## Mapping to the Haldane Model

Our experiment is performed in the regime where

In the weak spring limit, the motion of each gyroscope is approximately circular, so there is only one degree of freedom per gyroscope. This can be seen by splitting the displacement of the gyroscope two polarizations,

Substituting this form of *ψ* into the linearized equation of motion and matching the coefficients of the exponentials gives*u*’s alone, which is equivalent to Haldane’s model of the quantum Hall effect. Because **S17** implies **S16** gives*m* and *n*, and the last sum is over all neighbors *m* of *n* and neighbors *l* of *m*.

If one expresses the right side as **S18**, the bond angles lead to a phase shift between next-nearest neighbors on the lattice, similar to the phase shift on hopping terms from moving in a magnetic field. For the honeycomb lattice this differs from Haldane’s model only in that the second neighbor term has a real part.

The topological character of Haldane’s system can be quantified by calculation of the Chern number, which is an integral of the Berry curvature over the Brillouin zone of the lattice. A nonzero Chern number indicates a topologically nontrivial state and implies the existence of chiral edge currents. Systems with time-reversal symmetry must have a Chern number equal to zero, because time-reversal symmetry implies zero Berry curvature. However, not all systems with broken time-reversal symmetry must have a nonzero Chern number. In general, two bands separated by a finite gap will not acquire nonzero Chern number because of infinitesimal perturbations.

However, an infinitesimal perturbation can produce a large change in Berry curvature at Dirac points. Therefore, even the small complex phase terms in the Haldane model and in Eq. **S18** can open a gap and induce a nonzero Chern number.

To see this effect mathematically, one can expand in powers of displacement from the Dirac point, *T*, can be written in terms of Pauli matrices. For a system without next-nearest-neighbor coupling, we find that after rotating the wavevector,

For the honeycomb lattice there are only two sites per unit cell, so there are two modes with each wave number. A basis can be obtained by defining *a* is the edge length of the hexagon; Fig. S7 shows the two basis states near

Let us focus on just the leading terms of Eq. **S18** and the terms that arise from hopping along the diagonals (next-nearest neighbors,

Both wave functions pick up the same phase under translation, but they transform oppositely under rotation. Because the phase differences are different, the energies (i.e., frequencies of the normal modes in the gyroscope system) of the states are different; the one whose phase shifts match the phase of the hopping better has the lower energy (i.e., the energies are *r* next-nearest neighbors).

If we include the energy splitting and the linear terms in

The total Berry curvature for each of the two Dirac points is

In general, distorting the honeycomb lattice produces different phase shifts along different diagonals (Fig. S8). For a distorted honeycomb lattice, there are four angles of the hexagon equal to *α* and two equal to

## Acknowledgments

We acknowledge the Materials Research and Engineering Centers (MRSEC) Shared Facilities at The University of Chicago for the use of their instruments. This work was supported by the National Science Foundation MRSEC Program at The University of Chicago (Grant DMR-1420709). W.T.M.I. acknowledges support from the A. P. Sloan Foundation through a Sloan Fellowship and the Packard Foundation through a Packard Fellowship.

## Footnotes

- ↵
^{1}To whom correspondence may be addressed. Email: aturne26{at}jhu.edu or wtmirvine{at}uchicago.edu.

Author contributions: L.M.N., D.K., V.V., A.M.T., and W.T.M.I. designed research; L.M.N., D.K., A.R., A.M.T., and W.T.M.I. performed research; L.M.N., D.K., A.R., A.M.T., and W.T.M.I. analyzed data; and L.M.N., D.K., V.V., A.M.T., and W.T.M.I. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1507413112/-/DCSupplemental.

## References

- ↵.
- Maxwell JC

- ↵
- ↵.
- Thorpe MF

- ↵
- ↵
- ↵
- ↵
- ↵Paulose J, Chen BG-g, Vitelli V (2015) Topological modes bound to dislocations in mechanical metamaterials..
*Nat Phys*11(2):153–156 - ↵.
- Chen BG,
- Upadhyaya N,
- Vitelli V

- ↵.
- Paulose J,
- Meeussen AS,
- Vitelli V

- ↵Vitelli V, Upadhyaya N, Chen BG-g (2014) Topological mechanisms as classical spinor fields. arXiv:1407.2890 [cond-mat.soft].
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Süsstrunk R,
- Huber SD

- ↵.
- Peano V,
- Brendel C,
- Schmidt M,
- Marquardt F

- ↵.
- Ningyuan J,
- Owens C,
- Sommer A,
- Schuster D,
- Simon J

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Palacci J,
- Sacanna S,
- Steinberg AP,
- Pine DJ,
- Chaikin PM

- ↵
- ↵

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Applied Physical Sciences