Observation of non-Hermitian topology and its bulk-edge correspondence in an active mechanical metamaterial
Ananya Ghatak Martin Brandenbourger Jasper van Wezel Corentin Coulais
cond-mat.mes-hall cond-mat.str-el quant-ph






























    ∂̣|̌ł⟨⟩̊/1/2𝒫𝒯𝒫𝒯𝒞𝒞𝒫𝒯University of Amsterdam, Institute of Physics, Science Park 904, 1098 XH, Amsterdam, the NetherlandsUniversity of Amsterdam, Institute of Physics, Science Park 904, 1098 XH, Amsterdam, the NetherlandsUniversity of Amsterdam, Institute of Physics, Science Park 904, 1098 XH, Amsterdam, the NetherlandsUniversity of Amsterdam, Institute of Physics, Science Park 904, 1098 XH, Amsterdam, the Netherlands
        Observation of non-Hermitian topology and its bulk-edge correspondence in an active mechanical metamaterial
            Corentin Coulais
            December 30, 2023
        ===========================================================================================================
        
        
        Topological edge modes are excitations that are localized at the materials' edges and yet are characterized by a topological invariant defined in the bulk. Such bulk-edge correspondence has enabled the creation of robust electronic, electromagnetic and mechanical transport properties across a wide range of systems, from cold atoms to metamaterials, active matter and geophysical flows. Recently, the advent of non-Hermitian topological systems—wherein energy is not conserved—has sparked considerable theoretical advances. In particular, novel topological phases that can only exist in non-Hermitian systems have been introduced. However, whether such phases can be experimentally observed, and what their properties are, have remained open questions. Here, we identify and observe a novel form of bulk-edge correspondence for 
        a particular non-Hermitian topological phase. We find that a change in the bulk non-Hermitian topological invariant leads to a change of topological edge mode localisation together with peculiar purely non-Hermitian properties. Using a quantum-to-classical analogy, we create a mechanical metamaterial with non-reciprocal interactions, in which we observe experimentally the predicted bulk-edge correspondence, demonstrating its robustness. Our results open new avenues for the field of non-Hermitian topology and for manipulating waves in unprecedented fashions.Keywords. topological insulators, broken Hermiticity, mechanical metamaterials
        
        Significance statement.
        In recent years, the mathematical concept of topology has been used to predict and harness the propagation of waves such as light or sound in materials. However, these advances have so far been realized in idealised scenarios, where waves do not attenuate. In this research, we demonstrate that topological properties of a mechanical system can predict the localisation of waves in realistic settings where the energy can grow and/or decay. These findings may lead to novel strategies to manipulate waves in unprecedented fashions, for applications in vibration damping, energy harvesting and sensing technologies.
        
        
        
        § INTRODUCTION
        
        The inclusion of non-Hermitian features in topological insulators has recently seen an explosion of activity. Exciting developments include tunable wave guides that are robust to disorder <cit.>, structure-free systems <cit.>, and topological lasers and pumping <cit.>. In these systems, active components are introduced to typically: (i) break time-reversal symmetry to create topological insulators with unidirectional edge modes <cit.>; (ii) pump topologically protected edge modes, thus harnessing Hermitian topology in non-Hermitian settings <cit.>.
        In parallel, extensive theoretical efforts have generalized the concept of a topological insulator to truly non-Hermitian phases that cannot be realised in Hermitian materials <cit.>.
        However, such non-Hermitian topology and its bulk-edge correspondence remain a matter of intense debate. 
        On the one hand it has been argued that the usual bulk-edge correspondence breaks down in non-Hermitian settings, while on the other hand new topological invariants specific to non-Hermitian systems have been proposed to capture particular properties of their edge modes <cit.>.
        
        
        
        
        Here, focussing on a non-Hermitian version of the Su-Schrieffer-Heeger (SSH) model <cit.> with an odd number of sites (Fig. <ref>a), we find that a change in the bulk non-Hermitian topological invariant is accompanied by a localization change in the zero-energy edge modes. This finding suggests the existence of a
         bulk-edge correspondence for this new type of truly non-Hermitian topology. We further construct a mechanical analogue of the non-Hermitian quantum model (Fig. <ref>b), and create a mechanical metamaterial (Fig. <ref>c) in which we observe the predicted correspondence between the non-Hermitian topological invariant and the topological edge mode. In particular, we report that the edge-mode in the non-Hermitian topological phase has a peculiar nature, as it is localized on the rigid rather than the floppy side of the mechanical metamaterial.
        
        
        
        
        § NON-HERMITIAN WINDING NUMBER
        
        The one-dimensional model depicted schematically in Fig. <ref>a, is described by the quantum mechanical Bloch Hamiltonian
        
        H(k) =([    0 Q(k); R(k)    0;      ]),
        
        where k is the wave vector. The coefficients Q(k)= a_1+b_2 e^-ik and R(k)=a_2+b_1 e^ik allow electrons to hop between neighbouring sites within the unit cell (a_1 and a_2), as well as between unit cells (b_1 and b_2). If the amplitudes for hopping to the left (a_1 and b_1) are different from the corresponding amplitudes for hopping to the right (a_2 and b_2), the Hamiltonian is non-Hermitian, with complex eigenvalues E_±(k)=±√(Q(k) R(k)) that come in pairs related by reflection in the point E=0. Thus Eq. (<ref>) has a chiral symmetry,  (σ_z^-1H(k) σ_z = -H(k)) and falls in symmetry class AIII <cit.>. 
        
        A non-Hermitian Hamiltonian such as Eq. (<ref>) may host two different types of topological invariants, corresponding either to a winding of the phase of their eigenvectors as the wave vector k is varied across the Brillouin zone <cit.> (see Eq. (<ref>) of the Materials and Methods), or to the complex energies winding around one another in the complex energy plane <cit.> (see Eq. (<ref>) of the Materials and Methods). The former type of topology exists both for Hermitian and non-Hermitian systems, while the latter is exclusive to non-Hermitian systems, has not been observed yet, and is the focus of the present work.
        
        
        
        § MAPPING BETWEEN NON-HERMITIAN QUANTUM AND CLASSICAL MODELS
        
        The non-Hermitian topology contained in the model of Eq. (<ref>) stems from the non-reciprocity of its hopping parameters. This renders a direct implementation within a quantum material challenging, but recent advances on non-reciprocal mechanical metamaterials <cit.> suggest that such non-reciprocal interactions can be realized within a mechanical platform. In particular, inspired by the works of Kane and Lubensky <cit.> and Brandenbourger et al. <cit.>, we introduce the one-dimensional mechanical system (Fig. <ref>b), which is described by the dynamical matrix:
        
        D(k)=(-a+b e^-i k)(-a(1-ε )+b(1+ε )e^i k).
        
        Here, a=(p+2rsinθ)/√(p^2+4r^2cos^2θ) and b=(p-2rsinθ)/√(p^2+4r^2cos^2θ) are geometrical parameters that depend on the length r, the initial angle θ of the red rotors and the lattice spacing p (see Fig. <ref>b and Materials and Methods). The parameter ε modifies the stiffness of the blue springs in a non-reciprocal way, so that a strain in the spring causes a larger torque on the left rotor than on the right.  This non-reciprocal interaction is created locally for each robotic unit cell by an active-control loop: the motor of each unit cell applies a torque that depends on the strain of its neighbouring springs (see Materials and Methods).
        
        The equations of motion imposed by the dynamical matrix D(k) on the displacements and their time derivatives may be combined into a Schrödinger-like equation, as proposed by Kane and Lubensky <cit.>. The matrix taking the place of the Hamiltonian in this formulation has precisely the same form as Eq. (<ref>), with Q(k)=-a+be^-ik and R(k)=-a(1-ε)+b(1+ε )e^ik but the eigenvalues represent frequencies ω_±(k), rather than the energies (see Materials and Methods). This generalises the formal mapping between the dynamical matrix D(k) and the Hamiltonian H(k) introduced by Kane and Lubensky to a non-Hermitian setting.
        
        
        
        § NON-HERMITIAN BULK-EDGE CORRESPONDENCE
        
        In the following, we restrict our attention to a particular model with parameter values a_1=a, a_2=a(1-ε), b_1=-b(1+ε) and b_2=-b. 
        In the reciprocal and non-Hermitian limit of ε=0, the two bands of the model lie entirely on the real axis. As shown in Fig. <ref>, increasing ε leads to the bands developing imaginary components, and eventually touching at the exceptional point E=0 and k=0 before coalescing. The non-Hermitian winding invariant ν then becomesone. At even larger values of ε another exceptional point is encountered at k=π and the bands separate again into a non-winding phase.
        
        
        We show the full phase diagram of this system in Fig. <ref>a, as a function of the hopping parameters a/b and the non-reciprocal parameter ε. In the hatched pink region, the complex energies (frequencies) wind, and the non-Hermitian topological invariant has the value ν=1. The other region has non-winding non-Hermitian topology with ν=0, in accordance with the fact that the energies (frequencies) form disconnected bands. The phase boundaries correspond to the parameter values at which the bands coalesce at the exceptional point E=0 <cit.>, and are given by ε_c=(a/b± 1)/(a/b∓ 1). The ε=0 axis represents the Hermitian limit, in which the non-Hermitian invariant is always zero and the two branches of exceptional points combine into a Dirac point at a/b=1. For the classical system, this axis corresponds precisely to the Kane-Lubensky model <cit.>. The fact that the regions of non-Hermitian topology span large parts of parameter space suggests that they may be realized experimentally.
        
        
        The hatched pink regions of the phase diagram Fig. <ref>a are based on the behaviour of bulk topological invariants, calculated in a system with periodic boundary conditions. The non-Hermitian topology, however, is expected to be most visible experimentally in the emergence or suppression of edge modes localised at the edges of the chain. The edge modes can be found for the quantum (classical) model by solving Schrödinger's (Newton's) equation for zero-energy modes (see Materials and Methods). We focus in the following on a SSH chain with an odd number of sites (Fig. <ref>a) and on the mechanical Kane-Lubensky chain (Fig. <ref>b), for which the bulk-edge correspondences are strictly equivalent. Namely, we investigate a SSH (Kane-Lubensky) chain with N A-sites (rotors) and N-1 B-sites (springs). There is a vast literature on even-sized SSH chains and the choice of an odd-sized SSH chain can appear as less conventional. However, the topological nature of edge modes is the same for odd and even chains. Even in the Hermitian limit and for any given value of the topological invariant, it is the termination of the chain rather than the distinction between an odd and even number of atoms that determines the presence or absence of an edge mode. This can be clearly seen when considering for example the half-infinite chain, which is neither odd nor even and whose edge mode can be predicted to be present or absent based on knowledge of the invariant and the termination <cit.>. Last but not least, choosing the odd chain is necessary to obtain a formal mapping with the mechanical system.
        
        In the Hermitian limit ε=0, both the quantum and the classical chain always have a single zero mode, which is localized to the right (left) edge for a>b (a<b) (green contours in Fig. <ref>a). In the non-Hermitian case ε≠ 0, the zero mode changes sides precisely at the critical lines ε=ε_c of the bulk, periodic system (Fig. <ref>a). In all cases, we find that the tails of the edge modes become oscillatory for |ε|>1 (Fig. <ref>bc), as a consequence of imaginary contributions to their eigenvectors. Other choices of parameters will lead to a qualitatively similar correspondence between edge mode localisation and bulk winding (see Materials and Methods). Finally, we find that perturbations of the ideal model considered here, such as the inclusion of on-site potentials or mechanical bending interactions, progressively gap the system and suppress the zero-modes (see Materials and Methods).
        
        The coincidence between the change of the non-Hermitian winding number and the change of localisation of the zero modes demonstrates that these zero modes are topological and that a change of localisation corresponds to a topological transition.These topological zero modes have several peculiar properties, that are only possible because Hermiticity is broken:
        
        (i) increasing the non-reciprocity at a fixed value of the ratio a/b can be seen in Fig. <ref>a to cause two consecutive changes in the edge mode location,  one of which goes against the direction of the non-reciprocal bias; (ii) in the case of the quantum system, the shape of the phase diagram can not be explained by a simple argument involving the shifting of unit cells, as can be done in the Hermitian limit; (iii) in the mechanical system, as the topological edge mode in the winding region occurs where the mechanical degrees of freedom are constrained—this is a zero-energy mode, yet it involves stretching of the springs. 
        
        
        
        The bulk-edge correspondence shown in Fig. <ref>a differs from but is complementary to recent results on even non-Hermitian SSH chains, where the topological modes appear or disappear 
        at the values ε_s=((a/b)^2± 1)/((a/b)^2∓ 1)at which the gap of the open chain is closed <cit.>. That the gap closings of the open and closed chain do not coincide is a manifestation of the non-Hermitian skin-effect, which also causes all modes in the system except the topological zero mode to localise on one end of the chain
        
        
        (see Materials and Methods) <cit.>.Recent results show that taking into account the non-Hermitian skin-effect allows the definition of a non-Bloch-topological invariant,  which  switches  value  at ε_s <cit.> A  physically  compelling picture  thus  emerges  for  non-Hermitian  topological  phases: while the winding topological invariant predicts the edge at which the topological zero modes are localized,  the non-Bloch invariant predicts the existence of the topological modes and the location of the gap closing.
        
        § NON-HERMITIAN BULK-EDGE CORRESPONDENCE IN AN ACTIVE MECHANICAL METAMATERIAL
        
        To demonstrate the non-Hermitian bulk-edge correspondence described above we provide an experimental realization. To this end, we build an active mechanical metamaterial (Figs. <ref>c and <ref>a), which consists of nine robotic unit cells and in which a combination of geometry and active control is used to implement D(k), as defined in Eq. (<ref>). While the geometry allows us to obtain suitable values of a=1.00 and b=0.73, active control makes it possible to tune the non-reciprocal parameter ε (see Materials and Methods). We selectively access properties of the periodic (bulk) or open (edged) system by either including or omitting a rigid connection between the first and last unit cells (Fig. <ref>a and Materials and Methods).
        
        
        In this setup, we first perform relaxation experiments on the periodic metamaterial to quantify directly the bulk eigenfrequencies for the wave vectors k=0,π and hence the bulk topological invariant ν (Fig. <ref>a-d). We find that the non-Hermitian topological invariant jumps from zero—where the eigenfrequencies are disconnected as in Fig. <ref>b—to one—where the eigenfrequencies wind as in Fig. <ref>c—for a non-reciprocal parameter ε_c^periodic= 0.12 (Fig. <ref>d). Second, we probe the signature of the zero modes of the open chain (Fig. <ref>e), by applying a low-frequency excitation at the central unit cell. We observe a right-to-left (left-to-right) decaying displacement field for small (large) values of the non-reciprocal parameter ε (Fig. <ref>f and Supplementary Video). We find that the amplification factor |δθ_9/δθ_1| crosses the value one at a critical value of ε_c^open=0.12 (Fig. <ref>g). Remarkably, the correspondence ε_c^periodic=ε_c^open agrees precisely with the theoretically predicted non-Hermitian bulk-edge correspondence. It shows that the experimentally observable switching of edge state localization in the open chain coincides with the changing value of the non-Hermitian topological invariant in the (bulk) system with periodic boundary conditions. Moreover, it proves the robustness of both the bulk-boundary correspondence and the non-Hermitian topology to inherent deviations from the ideal model such as geometric and motor non-linearities, spring bending, time delays and noise in the micro-controllers, friction, and geometric irregularities. 
        
        To show more clearly the connection between the topological transition and the behaviour of the edge modes, we also create a domain wall in the metamaterial, with the leftmost part remaining reciprocal (ε=0) and the non-reciprocal parameter ε being tuned away from zero in the rightmost part (Fig. <ref>c) and vice-versa (Fig. <ref>d).  As expected, beyond the threshold value, the localization of the displacement field changes from the right edge to the domain boundary at the center (Fig. <ref>c) or the displacement field localises at both edges away from the domain boundary (Fig. <ref>d).
        
        
        
        
        § DISCUSSION AND OUTLOOK
        
        To conclude, we discovered and experimentally observed a novel type of bulk-edge correspondence for the non-Hermitian topological phase of a mechanical metamaterial with non-reciprocal interactions. This particular form of non-Hermitian bulk-edge correspondence, connected to energy winding, exhibits marked differences with the recently proposed non-Hermitian bulk-edge correspondence connected to a bi-orthogonal expectation value <cit.>. First, the correspondence based on energy winding reported here is unaffected by the non-Hermitian skin-effect: despite the complete reorganisation of the spectrum between a periodic and an open system, the energy winding of the periodic system predicts changes in the edge modes of the open system. 
        
        Second, the energy winding and the bi-orthogonal condition both predict the emergence of zero modes. However, while the bi-orthogonal condition predicts the existence of edge modes, the energy winding additionally predicts the side of the chain at which the topological mode appears. 
        
        
        These differences call for further investigation and generalization of the bulk-edge correspondence based on energy winding, beyond the particular system considered here.
        
        Further, we envision the study of nonlinearity, robustness to disorder, different interactions, higher spatial dimensions and other strategies to achieve non-Hermiticity—to be exciting future research directions. We believe that our work provides conceptual and technological advances, opening up avenues for the topological design of tunable wave phenomena. 
        
        
        
        
        
        
        
        
        § MATERIALS AND METHODS
        Hermitian and non-Hermitian topology of the non-reciprocal SSH model.
        The Hamiltonian of Eq. (<ref>) of the Main Text may acquire topological character either from the winding of the Berry connection determined by its eigenfunctions |ψ^±(k)⟩, or from the direct winding of its eigenenergies E^±(k) in the complex plane. While the former type of winding corresponds to conventional Hermitian topology <cit.>, the latter corresponds to a unique form of non-Hermitian topology, that can only exist when the eigenenergies are complex <cit.>.
        
        The Hamiltonian in Eq. (<ref>) is pseudo-Hermitian with respect to a positive definite metric operator <cit.> and can be diagonalized to find its eigenenergies E^±(k)=±√(R(k)Q(k)) as well as its left and right eigenmodes,
        ⟨ψ_L^+(k)|=(
        1 , √(Q(k)/R(k))), 
        ⟨ψ_L^-(k)|=(-√(R(k)/Q(k)),
        1 ), 
        |ψ_R^+(k)⟩=(1,√(R(k)/Q(k)))^T 
        and |ψ_R^-(k)⟩=(-√(Q(k)/R(k)),1 )^T. 
        These eigenmodes obey the bi-orthonormality condition and the Hamiltonian preserves a generalized unitarity condition <cit.>.
        
        To compute the conventional Hermitian topological invariant in the non-Hermitian setting, we can define the Berry connection as <cit.>𝒜^±(k)=-i⟨ψ_L^±(k)|∂_k ψ_R^±(k)⟩,
        from which the topological invariant is then calculated to be
        γ^± =1/2 π∫_0^2πdk A^±(k). 
        
        Choosing different left and right eigenvectors in the definition of the Berry connection, 
        or introducing a modified inner product, does not yield additional invariants <cit.>. The invariants γ^± are zero (when a/b >1) or integer (when a/b <1) for the regions with non-winding non-Hermitian 
        topology (non-hatched regions in Fig. <ref>a of the Main Text). In the winding non-Hermitian topological regions (pink hatched regions in Fig. <ref>a), the two bands coalesce into a 4π-periodic structure. In that case the Berry connection winds an integer number of times after integration over the full 4π period, and we find 1/2 π∫_0^4πdk A^±(k)=1 for the single, coalesced band.
        
        The truly non-Hermitian topology can be defined in terms of the winding of eigenenergies around an exceptional point, as shown in Fig. <ref> of the Main Text. The corresponding winding number can be calculated to be ν=-1/2 π∫_0^4πdk ∂/∂ k(  [E^+(k)-E^-(k)]).
        
        The value of this non-Hermitian topological invariant is shown in Fig. <ref>a of the Main Text. Notice that the definition of Eq. (<ref>) differs by a factor of two from the convention used in some other works <cit.>.
        
        
        Edge modes of the non-reciprocal SSH model with an odd number of sites.
        In Eq. (<ref>) of the Main Text, a non-reciprocal version of the SSH model is defined in reciprocal space. Here, we use the corresponding real space formulation to identify the edge modes of a finite non-reciprocal SSH chain with open boundary conditions. Specifically, we consider N sites of type A, and N-1 sites of type B, which is strictly analogous to the mechanical Kane-Lubensky chain <cit.>. The Hamiltonian is given in real space by:
        
            H=(
            [   0 a_1   0       …       0; a_2   0 b_1                ;   0 b_2   0 a_1   ⋱       ⋮;         a_2   0   ⋱        ;   ⋮       ⋱   ⋱   ⋱ a_1   0;                 a_2   0 b_1;   0       …       0 b_2   0;     ]),
        
        where a_1, a_2, b_1, and b_2 are hopping coefficients, which we assume to be real. In the Hermitian case, with a_1=a_2 and b_1=b_2, the Hamiltonian has a unique zero mode satisfying the equation H |ψ⟩= 0, with the eigenmode written as |ψ⟩=(ψ^A_1,ψ^B_1,…,ψ^A_N-1,ψ^B_N-1,ψ^A_N)^T. However, in the generic non-Hermitian case considered here, with a_1≠ a_2 and b_1≠ b_2, the equations H |ψ_R⟩= 0 and ⟨ψ_L|  H= 0 respectively for the right and left eigenvectors may yield distinct zero-energy modes. We solve these two equations and find (ψ^A_n)_R/(ψ^A_1)_R =(-a_2/b_1)^n-1, (ψ^B_n)_R=0 for the right eigenmode, and (ψ^A_n)_L /(ψ^A_1)_L =(-a_1/b_2)^n-1, (ψ^B_n)_L=0 for the left eigenmode. In Fig. <ref>b of the Main Text, we only plot the right eigenmodes for a_1=a, a_2=a(1-ε), b_1=-b(1+ε) and b_2=-b.
        
        
        Non-Hermitian skin effect. In the non-Hermitian case, boundary conditions have a significant effect on the shape of the entire spectrum, namely open boundaries shift modes at all energies (frequencies) towards one side of the chain, in what is known as the non-Hermitian skin effect <cit.>. This effect is not related to topology, and was recently observed in both a non-reciprocal mechanical metamaterial <cit.> and a non-reciprocal electronic circuit <cit.>. 
        Theoretically, it has been shown that in a non-Hermitian SSH chain, the closing of the bandgap appears at parameter values that are different for open and periodic boundary conditions, in an apparent breakdown of the bulk-edge correspondence <cit.>(indicated by the grey region in Fig. <ref>b). 
        We find consistent results (indicated by the grey region in Fig. <ref>a), but we report in addition a clear correspondence between the topology of the bulk spectrum computed with closed boundary-conditions, and the zero-energy edge modes obtained with open boundary conditions: (i) a zero-energy edge mode always exists, as calculated analytically in the section above and confirmed numerically (Fig. <ref>a); (ii) this edge mode in the chain with open boundaries changes its localization at the exceptional points of the periodic—bulk— model (See Fig. <ref>ab of the Main Text); (iii) surprisingly, this edge mode is unaffected by the gap closing ε=(a^2/b^2± 1)/(a^2/b^2∓ 1) of the open chain <cit.> (dashed lines in Fig. <ref>a of the Main Text and Fig. <ref>). These results are consistent with and complementary to recent results in the case of the even non-Hermitian SSH chain <cit.>
        
        
        
        Notice that in our system, a signature of the non-Hermitian skin effect can be seen in the response to a local excitation. While the localization of the response at low frequency (Fig. <ref>a) is essentially the same as that of the zero-mode (Fig. 3 of the main text), changing localization edge at the topological phase boundaries, the localization at large frequencies (Fig. <ref>b) solely depends on the non-reciprocal parameter, which is a direct signature of the non-Hermitian skin effect. For extensive portions of parameter space, the localization of the topological zero mode dominating the low-frequency response is opposite to that induced by the non-Hermitian skin effect at high frequencies.
        
        
        
        Non-reciprocal Kane-Lubensky chain.
        The classical analogue of the non-Hermitian SSH chain is a non-reciprocal version of the Kane-Lubensky chain <cit.>, as shown in Fig. <ref>b of the Main Text. In this system, N rotors of length r, with an initial tilt angle θ and a staggered orientation, are connected by N-1 springs between subsequent rotors, with lattice spacing p. 
        To construct the equation of motion for such a system, we first write the relation between the angular displacements of the rotors, |δθ⟩=(δθ_1,…,δθ_N)^T, and the length change of the springs (positive for stretching, negative for compression), |δℓ⟩=(δℓ_1,…,δℓ_N-1)^T. This is given by 
        |δℓ⟩= R |δθ⟩, with the compatibility matrix <cit.>:
        
            R = r cosθ([ -a  b  0  …  0;  0 -a  b     ⋮;  ⋮     ⋱  ⋱  0;  0  …  0 -a  b ]).
        
        The coefficients a= (p- 2rsinθ)/√(p^2+4r^2cos^2θ) and b=(p+ 2rsinθ)/√(p^2+4r^2cos^2θ) are geometric parameters
         <cit.>. We can similarly write the relation between the torque on each rotor, |τ⟩=(τ_1,…,τ_N)^T, and the tension in the springs, |t ⟩=(t_1,…,t_N-1)^T, in the form |τ⟩= Q |t ⟩, where Q is the equilibrium matrix. In the Hermitian case R and Q are transposes of each other <cit.>. The compatibility and equilibrium matrices can be multiplied to compute the so-called dynamical matrix D^Hermitian=k/JQR, where k is the spring constant of the elastic link between subsequent rotors and J is the rotational moment of inertia of the rotors. 
        
        In the active mechanical metamaterial, each unit cell n has a local control system that can apply a local torque τ^control loop_n that depends on the angular displacement of its rotor δθ_n and that of its neighbours δθ_n-1 and δθ_n+1. We choose to apply the following torque
        τ^control loop_n= ε k r^2cos^2θ (b (bδθ_n-aδθ_n-1) -  a (aδθ_n-bδθ_n+1)). Since the added torques are proportional to the angular displacements, the system can still be described by an effective dynamical matrix D=D^Hermitian+D^control loop of the form
         D= k r^2cos^2θ/J([          a^2 (1-ε )           -ab(1+ε )                   0                   0                   -                   -                   0;           -ab(1-ε ) a^2(1-ε )+b^2(1+ε )           -ab(1+ε )                   0                   -                   -                   0;                   0           -ab(1-ε ) a^2(1-ε )+b^2(1+ε )           -ab(1+ε )                   -                   -                   0;                                                           ⋱                   ⋱                   ⋱                                        ;                   0                   0                   -                   -                   0           -ab(1-ε )           b^2(1+ε ) ])
        ,
        
        which is non-symmetric for ε≠ 0.
        In the Main Text, we assume that the ratio k r^2cos^2θ/J=1 without loss of generality. 
         
        We compute the right and left zero modes by solving D|δθ_R ⟩= 0 and ⟨δθ_L |  D= 0, and we find (δθ_n)_R/(δθ_1)_R=(a(1-ε)/b(1+ε))^n-1 and (δθ_n)_L/(δθ_1)_L=(a/b)^n-1, respectively. We show and discuss only the right zero modes in Fig. <ref>bc and in the Main Text, because  the right zero modes dominate the observed angular displacement profile <cit.>.
        
        With periodic boundary conditions, the Fourier transform of Eq. (<ref>) becomes
        D(k)=a^2(1-ε )+b^2(1+ε) -ab(1+ε )e^ik-ab(1-ε )e^-ik, which in its factored form D(k)=Q(k) R(k), with Q(k)=(-a+b e^-ik), R(k)=(-a(1-ε)+b(1+ε) e^ik), coincides with Eq. (<ref>) of the Main Text. This factored mathematical expression allows us to construct the mapping between quantum and classical systems, following References <cit.>, where the quantum Hamiltonian is written as in Eq. (<ref>) of the Main Text.
        Notice that the physical meaning of the Fourier equilibrium and compatibility matrices Q(k) and R(k) is ill-defined in the non-Hermitian case.
        
        
        Role of perturbations
        The computational model introduced in Eq. (<ref>) and discussed in the Main Text is necessarily an idealisation. In the actual mechanical, integrated system in the experimental setup there are unavoidable small effects of bending in each of the elastomeric bands, on top of other essential effects from frictional forces, geometrical and electromechanical non-linearities in the chain, time delays and noise from micro-controllers and geometric irregularities. That we nevertheless see the theoretically predicted bulk-edge correspondence in our experimental results is thus witness to the robustness of the non-Hermitian topology described by the numerical predictions.
        
        To test the limits of the topological robustness, we explicitly probe the role of two types of perturbations in the numerical model:  (i) bending interactions inherently present in the elastomeric bands connecting the nearest neighbour rotors; (ii) an on-site potential.
        With these, the Fourier-transformed dynamical matrix becomes D(k)=a^2(1-ε +ρ+g)+b^2(1+ε+ρ+g) -ab(1+ε -ρ)e^ik-ab(1-ε -ρ)e^-ik], where ρ is the relative bending stiffness of each rubber band and g is the on-site potential. For the system with periodic boundary conditions, described by D(k), the phase boundaries between non-winding non-Hermitian topology and winding non-Hermitian topology become
        
            ε_c=a/b± 1/a/b∓ 1+ρa/b∓ 1/a/b± 1+g(a/b)^2+ 1/(a/b)^2- 1.
        
        Fig. <ref>ab shows that the presence of bending gaps the spectrum of the model with periodic boundary conditions around ε=0 for all values of the ratio a/b. At higher values of non-reciprocal parameter ε, however, the transition into a phase with non-Hermitian topology survives. Likewise, for non-zero on-site potential the gapping of the spectrum in the periodic chain causes a breakdown of the bulk-boundary correspondence around the point a/b=1 and ε=0, while for higher values of ε it survives (Fig. <ref>cd). Both of these observations are a testament to the robustness of the non-Hermitian topology and its bulk-boundary correspondence, which remain intact even in the presence of perturbations, for sufficiently large values of the non-reciprocity. These considerations have been taken into account in the design of the experiments described below, wherein a specific shape of the rubber band is chosen to minimize the bending.
        Experimental platform
        To perform the experiments, we followed Brandenbourger et al. <cit.>, and created a one-dimensional, non-reciprocal active mechanical metamaterial (Fig. <ref>c) consisting of nine unit cells, each of which has a single rotational degree of freedom θ_n, where n is the unit cell index. 
        The unit cells are mechanically connected in a specific geometry initially proposed in the paper of Kane and Lubensky <cit.> (Fig. <ref>b). Each unit cell consists of a rigid rotor of length 36 mm, and initial angle (-1)^n(θ-π/2) with respect to the horizontal axis, where θ=π/12. Each rotor is connected to its neigbours at its midpoint r=18 mm by a flexible, laser-cut elastomeric band of thickness 4 mm, whose shape has been optimized such that the band can easily hinge at its anchoring points and such that in-plane compression and stretching of the band dominates its elastic response (Fig. <ref>c). The lattice spacing between subsequent rotors is p=60 mm. Assuming that the bending of the band and friction can be neglected, to linear order, the deformations of the elastomeric bands induce a torque on rotor n given by τ^elastic_n= kr^2cos^2θ(  a (aδθ_n-bδθ_n+1) -  b (bδθ_n-aδθ_n-1)), where δθ_n is the angular displacement with respect to the initial angle in the reference state, k=7.8× 10^-3 N.m is the spring constant of the rubber band and a=1.00 and b=0.73 are the geometric parameters defined above.  
        
        In addition, each robotic unit cell is made of a mechanical oscillator, an angular encoder (Broadcom HEDR-55L2-BY09), a DC coreless motor (Motraxx CL1628), a microcontroller (Espressif ESP 32) and a custom-made electronic board that connects these components and ensures power conversion and communication between neighbouring unit cells. Each robotic unit cell uses a local active control loop to apply a strain-dependent torque on rotor n given by τ^control loop_n= ε k r^2cos^2θ (b (bδθ_n-aδθ_n-1) -  a (aδθ_n-bδθ_n+1)), where ε is a tunable dimensionless non-reciprocal gain parameter. The total torque on the rotors then realizes effective non-reciprocal interactions,
        
            τ_n=k r^2cos^2θ(   -ab(1-ε)δθ_n-1
               +(a^2(1-ε)+b^2(1+ε))δθ_n
               -ab(1+ε)δθ_n+1),
         
        where τ_n=τ^elastic_n+τ^control loop_n. 
        This system precisely realizes the dynamical matrix D defined in Eq. (<ref>), and therefore its mechanical response exhibits the bulk-edge correspondence shown by the right zero modes of D. We record the rotors' instantaneous positions δθ_n(t) via the angular encoders at a resolution 4.4× 10^-4 rad and sampling frequency 100 Hz.
        
        
        Measurements of the non-Hermitian winding number
        In order to measure the winding of the bands in the spectrum of the system with periodic boundary conditions, we connect the first and the last rotors with a rigid bar and ball-bearing hinges. To ensure homogeneity of the moments of inertia throughout the system, we add small masses at the end of each rotor.
        
        A rigid pin is attached to each rotor. We impose the initial position of each rotor δθ_n(t=0) away from their equilibrium position as follows: for each wave vector k, we manufacture drilled rigid plates with precisely positioned holes, in which the rigid pins can be inserted and that set the initial condition. The initial condition δθ_n(t=0)=δ̅θ̅cos k n is chosen such that the overall configuration of the chain is modified from equilibrium according to a specific wave vector (k=0 or k=π). In order to stay in the linear regime as well as within the limit of angular resolution, we impose δ̅θ̅ = 0.21 rad and 0.04 rad, respectively, for the measurements on the wave vectors k=0 and k=π. We remove rapidly the drilled plate to let the system freely relax. For each experiment, every unit cell is observed to relax the same way, except for their phase. For each wave vector k, we fit the displacement overtime of the far left unit cell to the equation (A_1 exp(λ_k t)+ A_2 exp(λ'_k t) ) cos(ω_k t) to deduce the real (imaginary) parts of the eigenfrequencies ω_k (λ_k and λ'_k), see Figs. <ref>bc of the Main Text. From there, we use a discretized version of Eq. (<ref>)
        
        ν=-2/π(arctanλ_π-λ'_π/2ω_π-arctanλ_0-λ'_0/2ω_0)
        
        to compute the winding invariant as a function of ε, as shown in Fig. <ref>d of the Main Text.
        
        
        Measurements of the edges modes
        We excite the metamaterial at the center rotor (n=5) by applying sinusoidal oscillations of amplitude 0.028 rad and frequency 0.05 Hz over 5 periods of oscillations. We extract the magnitude of the oscillation of each rotor via a Fourier series analysis to produce the data shown in Fig. <ref>f-e of the Main Text.
        
        Acknowledgments.—We thank D. Giesen, T. Walstra and T. Weijers for their skilful technical assistance. We are grateful to C. Bender, L. Fu, T. Lubensky, and D. Z. Rocklin, R. Thomale and Z. Wang for insightful discussions and to M.S. Golden for critical review of the manuscript. J.v.W. acknowledge funding from the a Netherlands Organization for Scientific Research (NWO) VIDI grant. C.C. acknowledges funding from the European Research Council grant ERC-StG-Coulais-852587-Extr3Me.
        
        
        pnas-new
        10 url<#>1 urlprefixURL Nash_PNAS2015authorNash, L. M.et al.titleTopological mechanics of gyroscopic metamaterials.
        journalProc. Natl. Ac. Sc. U. S. A.volume112, pages14495
          (year2015).
        
        Mitchell_NatPhys2018authorMitchell, N. P., authorNash, L. M.,
          authorHexner, D., authorTurner, A. M.&authorIrvine, W. T. M.titleAmorphous topological insulators constructed from
          random point sets.
        journalNat. Phys.volume14,
          pages380 (year2018).
        
        Khanikaev_NatComm2015authorKhanikaev, A. B., authorFleury, R.,
          authorMousavi, S. H.&authorAlù, A.titleTopologically robust sound propagation in an
          angular-momentum-biased graphene-like resonator lattice.
        journalNat. Comm.volume6,
          pages8260 (year2015).
        
        Souslov_PRL2019authorSouslov, A., authorDasbiswas, K.,
          authorFruchart, M., authorVaikuntanathan, S.&authorVitelli, V.titleTopological waves in fluids with odd viscosity.
        journalPhys. Rev. Lett.volume122, pages128001
          (year2019).
        
        Delplace_Science2017authorDelplace, P., authorMarston, J. B.&authorVenaille, A.titleTopological origin of equatorial waves.
        journalSciencevolume358,
          pages1075 (year2017).
        
        Bandres_Science2018authorBandres, M. A.et al.titleTopological insulator laser: Experiments.
        journalSciencevolume359,
          pageseaar4005 (year2018).
        
        Zilberberg_Nature2018authorZilberberg, O.et al.titlePhotonic topological boundary pumping as a probe of
          4d quantum hall physics.
        journalNaturevolume553,
          pages59 (year2018).
        
        Kraus_PRL2012authorKraus, Y. E., authorLahini, Y.,
          authorRingel, Z., authorVerbin, M.&authorZilberberg, O.titleTopological states and adiabatic pumping in
          quasicrystals.
        journalPhys. Rev. Lett.volume109, pages106402
          (year2012).
        
        Lohse_NatPhys2015authorLohse, M., authorSchweizer, C.,
          authorZilberberg, O., authorAidelsburger, M.&authorBloch, I.titleA thouless quantum pump with ultracold bosonic atoms
          in an optical superlattice.
        journalNat. Phys.volume12,
          pages350 (year2015).
        
        Pedro_PRL2019authorPedro, R. P., authorPaulose, J.,
          authorSouslov, A., authorDresselhaus, M.&authorVitelli, V.titleTopological protection can arise from thermal
          fluctuations and interactions.
        journalPhys. Rev. Lett.volume122, pages118001
          (year2019).
        
        Zeuner_PRL2015authorZeuner, J. M.et al.titleObservation of a topological transition in the bulk
          of a non-hermitian system.
        journalPhys. Rev. Lett.volume115, pages040402
          (year2015).
        
        Gong_PRX2018authorGong, Z.et al.titleTopological phases of non-hermitian systems.
        journalPhys. Rev. Xvolume8,
          pages031079 (year2018).
        
        Shen_PRL2018authorShen, H., authorZhen, B.&authorFu,
          L.titleTopological band theory for non-hermitian
          hamiltonians.
        journalPhys. Rev. Lett.volume120, pages146402
          (year2018).
        
        Ghatak_JPhysCondMat2019authorGhatak, A.&authorDas, T.titleNew topological invariants in non-hermitian systems.
        journalJ. Phys.: Cond. Mattervolume31, pages263001
          (year2019).
        
        Yao_PRL2018authorYao, S.&authorWang, Z.titleEdge states and topological invariants of
          non-hermitian systems.
        journalPhys. Rev. Lett.volume121, pages086803
          (year2018).
        
        Kunst_PRL2018authorKunst, F. K., authorEdvardsson, E.,
          authorBudich, J. C.&authorBergholtz, E. J.titleBiorthogonal bulk-boundary correspondence in
          non-hermitian systems.
        journalPhys. Rev. Lett.volume121, pages026808
          (year2018).
        
        Lee_PRB2019authorLee, C. H.&authorThomale, R.titleAnatomy of skin modes and topology in non-hermitian
          systems.
        journalPhys. Rev. Bvolume99,
          pages201103(R) (year2019).
        
        Helbig_2019authorHelbig, T.et al.titleGeneralized bulk–boundary correspondence in non-Hermitian topolectrical circuits. journalNat. Phys.volume16,
          pages747–750 (year2020).
        
        Bergholtz_review2019authorBergholtz, E. J., authorBudich, J. C.&authorKunst, F. K.titleExceptional topology of non-hermitian systems.
        journalarXiv:1912.10048  (year2019).
        
        Xiao_NatPhys2020authorXiao, L.et al.titleNon-hermitian bulk–boundary correspondence in
          quantum dynamics. journalNat. Phys.volume16,
          pages761–766 (year2020).
        
        MartinezAlvarez_PRB2018authorMartinez Alvarez, V. M., authorBarrios Vargas,
          J. E.&authorFoa Torres, L. E. F.titleNon-hermitian robust edge states in one dimension:
          Anomalous localization and eigenspace condensation at exceptional points.
        journalPhysical Review Bvolume97 (year2018).
        
        Yoshida_arxiv2019authorYoshida, T.&authorHatsugai, Y.titleExceptional rings protected by emergent symmetry for
          mechanical systems.
          journalPhys. Rev. Bvolume100,
          pages054109 (year2019).
        
        Hasan_RMP2010authorHasan, M. Z.&authorKane, C. L.titleColloquium: Topological insulators.
        journalRev. Mod. Phys.volume82, pages3045 (year2010).
        
        Fleury_Science2014authorFleury, R., authorSounas, D. L.,
          authorSieck, C. F., authorHaberman, M. R.&authorAlù, A.titleSound isolation and giant linear nonreciprocity in a
          compact acoustic circulator.
        journalSciencevolume343,
          pages516 (year2014).
        
        Coulais_Nature2017authorCoulais, C., authorSounas, D.&authorAlù, A.titleStatic non-reciprocity in mechanical metamaterials.
        journalNaturevolume542,
          pages461 (year2017).
        
        Wang_PRL2018authorWang, Y.et al.titleObservation of nonreciprocal wave propagation in a
          dynamic phononic lattice.
        journalPhys. Rev. Lett.volume121, pages194301
          (year2018).
        
        Kane_NatPhys2014authorKane, C. L.&authorLubensky, T. C.titleTopological boundary modes in isostatic lattices.
        journalNat. Phys.volume10,
          pages39 (year2014).
          
        Brandenbourger_arxiv2019authorBrandenbourger, M., authorLocsin, X.,
          authorLerner, E.&authorCoulais, C.titleNon-reciprocal robotic metamaterials.
        journalNat. Commun.volume10,
          pages4608 (year2019).
        
        
        Lubensky_RepProgPhys2015authorLubensky, T. C., authorKane, C. L.,
          authorMao, X., authorSouslov, A.&authorSun, K.titlePhonons and elasticity in critically coordinated
          lattices.
        journalRep. Prog. Phys.volume78, pages073901
          (year2015).
        
        Huber_NatPhys2016authorHuber, S. D.titleTopological mechanics.
        journalNat. Phys.volume12,
          pages621 (year2016).
        
        Susstrunk_PNAS2016authorSusstrunk, R.&authorHuber, S. D.titleClassification of topological phonons in linear
          mechanical metamaterials.
        journalProc. Natl. Ac. Sc. U. S. A.volume113, pagesE4767
          (year2016).
          
        history_topo
        
        
        
        
        
        
        
        
        
        
        
        
        
        
        
            authorZak, J.titleSymmetry criterion for surface states in solids.
        journalPhys. Rev. Bvolume32,
          pages2218 (year1985).
          
        McDonald_PRX2018authorMcDonald, A., authorPereg-Barnea, T.&authorClerk, A. A.titlePhase-dependent chiral transport and effective
          non-hermitian dynamics in a bosonic kitaev-majorana chain.
        journalPhys. Rev. Xvolume8,
          pages041031 (year2018).
        
        Mostafazadeh_IJGMMP2010authorMostafazadeh, A.titlePseudo-hermitian representation of quantum
          mechanics.
        journalInt. J. Geom. Meth. Mod. Phys.volume07, pages1191 (year2010).
        
        Borgnia_arxiv2019authorBorgnia, D. S., authorKruchkov, A. J.&authorSlager, R.-J.titleNon-Hermitian Boundary Modes and Topology.
        journalPhys. Rev. Lett.volume124pages056802  (year2020).
        
        SchomerusauthorSchomerus, H.titleNonreciprocal response theory of non-hermitian
          mechanical metamaterials: Response phase transition from the skin effect of
          zero modes.
        journalPhys. Rev. Res.volume2pages013058 (year2020).
