https://www.nature.com/articles/s41598-021-04178-8 Skip to main content Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript. Advertisement Advertisement Scientific Reports * View all journals * Search * My Account Login * Explore content * About the journal * Publish with us * Sign up for alerts * RSS feed 1. nature 2. scientific reports 3. articles 4. article Non-Hermitian topology in rock-paper-scissors games Download PDF * Article * Open Access * Published: 12 January 2022 Non-Hermitian topology in rock-paper-scissors games * Tsuneya Yoshida^1, * Tomonari Mizoguchi^1 & * Yasuhiro Hatsugai^1 Scientific Reports volume 12, Article number: 560 (2022) Cite this article * 692 Accesses * 2 Citations * Metrics details Subjects * Biological physics * Evolution * Topological insulators Abstract Non-Hermitian topology is a recent hot topic in condensed matters. In this paper, we propose a novel platform drawing interdisciplinary attention: rock-paper-scissors (RPS) cycles described by the evolutionary game theory. Specifically, we demonstrate the emergence of an exceptional point and a skin effect by analyzing topological properties of their payoff matrix. Furthermore, we discover striking dynamical properties in an RPS chain: the directive propagation of the population density in the bulk and the enhancement of the population density only around the right edge. Our results open new avenues of the non-Hermitian topology and the evolutionary game theory. Download PDF Introduction In these decades, the notion of the topology has played a central role in condensed matter physics^1,2,3,4,5. Analysis of topological aspects of condensed matters dates back to integer quantum Hall systems^6 which exhibit robust chiral edge modes induced by the topology of the Hermitian Hamiltonian in the bulk^7,8,9. While these robust edge modes are originally reported for quantum systems, they have been extended to various disciplines of science^10,11,12,13,14, 15,16,17,18,19,20. In particular, the notion of topology has been extended to an interdisciplinary field^21,22; the emergence of topological edge modes has been predicted^22 for networks of rock-paper-scissors (RPS) cycles which are described by the evolutionary game theory^23,24,25,26,27,28. Among the variety of topological phenomena, non-Hermitian topology^29 ,30,31,32,33,34 has become one of the recent hot topics because it results in novel phenomena which do not have Hermitian counterparts. A representative example is the emergence of exceptional points (EPs) ^35,36,37,38 and their symmetry-protected variants^39,40,41,42,43,44, 45 on which eigenvalues of the non-Hermitian Hamiltonian touch both for the real- and imaginary-parts. Another typical example is a non-Hermitian skin effect^46,47,48,49,50,51,52 which is extreme sensitivity of the eigenvalues and eigenstates to the presence/ absence of boundaries. As is the case of the Hermitian topology, the above non-Hermitian phenomena have also been reported in a variety of systems^53,54,55,56,57,58,59,60,61,62 (e.g., photonic systems^53,54, 55,56 and quantum systems^57,58,59,60,63). Despite the above significant progress, topological phenomena of the evolutionary game theory, which attract interdisciplinary attention, are restricted to the Hermitian topology. Highlighting non-Hermitian topology of such systems is significant as it may provide a new insights and may open a new avenue of the evolutionary game theory. In this paper, we report non-Hermitian topological phenomena in the evolutionary game theory: an EP and a skin effect in RPS cycles. The EP in the single RPS cycle is protected by the realness of the payoffs which is mathematically equivalent to parity-time (PT) symmetry. Our linearized replicator equation elucidates that the EP governs dynamics of the RPS cycle. Furthermore, we discover striking dynamical phenomena in an RPS chain induced by the skin effect: the directive propagation of the population density in the bulk and the enhancement of the population density only around the right edge. These dynamical properties are in sharp contrast to those in the Hermitian systems. The above results open new avenues of the non-Hermitian topology and the evolutionary game theory. Results EP in a single RPS cycle Firstly, we demonstrate the emergence of an EP at which two eigenvalues touch both for the real- and imaginary-parts. Consider two players play the RPS game (see Fig. 1a) whose payoff matrix is given by^64,65,66,67 $$\begin{aligned} A(\lambda )= & {} \left( \begin{array}{ccc} 0 &{} -1 &{} 1 \\ 1 &{} 0 &{} -1 \\ -1 &{} 1 &{} 0 \end{array} \right) +\ lambda \left( \begin{array}{ccc} 0 &{} 0 &{} 0 \\ 0 &{} 1 &{} -1 \\ 0 &{} -1 &{} 1 \end{array} \right) , \end{aligned}$$ (1) with a real number \(\lambda\). Here, players choose one of the strategies \((s_1,s_2,s_3)=(``\mathrm {R}",``\mathrm {P}",``\mathrm {S}")\). The payoff of a player is \(A_{IJ}\) when the player chooses the strategy \(s_I\) and the other player chooses \(s_J\) (\(I,J= 1,2,3\)). For \(\lambda =0\), this game is reduced to the standard zero-sum RPS game where the sum of all players' payoff is zero for an arbitrary set of strategies. The above payoff matrix exhibits an EP, which can be deduced by noting the following two facts. (1) The first (second) term is anti-Hermitian (Hermitian). (2) For an arbitrary \(\lambda\), eigenvalues \(\epsilon _n\) (\(n=1,2,3\)) form a pair \(\epsilon _n=\ epsilon ^*_{n'}\) (\(n\ne n'\)) or are real numbers \(\epsilon _n \in \mathbb {R}\). This constraint arises from the realness of the payoffs $$\begin{aligned} \mathcal {K} A(\lambda ) \mathcal {K}= & {} A(\ lambda ), \end{aligned}$$ (2) where the operator \(\mathcal {K}\) takes complex conjugate. Equation (2) can be recognized as PT symmetry by regarding \(\lambda \) as a momentum (for more details, see Sect. S1 of Supplemental Material [68]). For \(\lambda =0\), the energy eigenvalues are aligned along the imaginary axis due to the anti-Hermiticity of \(A(\ lambda =0)\). Increasing \(\lambda\), two eigenvalues touch at a critical value \(\lambda _\mathrm {c}\) so that they become real numbers when the second term is dominant. At \(\lambda =\lambda _{\ mathrm {c}}\), the EP emerges. The emergence of the EP is supported by Fig. 1b. For \(\lambda =0\), the energy eigenvalues are pure imaginary due to anti-Hermiticity of \(A(\lambda =0)\). As \(\lambda\) is turned on, two eigenvalues approach each other, and the EP emerges at \(\lambda =\lambda _c=\ sqrt{3}\). One can also characterize the topology of this EP by computing the \(\mathbb {Z}_2\)-invariant \(\nu\)^45 $$\begin{aligned} \nu (\lambda )= & {} \mathrm {sgn}\left[ \mathrm {Disc}_E P(E,\lambda ) \right] , \end{aligned}$$ (3) with \(P(E,\lambda )=\mathrm {det}[A(\lambda )-E\mathbbm{1}]\). Here, \(\mathrm {sgn}(x)\) takes 1 (\(-1\)) for \(x>0\) (\(x<0\)), and \(\ mathrm {Disc}_E P(E,\lambda )\) denotes discriminant of \(P(E,\lambda )\). For the \(3\times 3\)-matrix \(A(\lambda )\), \(\mathrm {Disc}_E P(E,\lambda )=(\epsilon _1-\epsilon _2)^2(\epsilon _1-\epsilon _3)^2 (\epsilon _2-\epsilon _3)^2\) holds. Figure 1c plots the \(\mathbb {Z}_2\)-invariant as a function of \(\lambda\). Corresponding to the emergence of the EP, the \(\mathbb {Z}_2\)-invariant jumps at \(\ lambda =\lambda _{\mathrm {c}}\), elucidating the topological protection of the EP. Figure 1 figure 1 (a) Sketch of the RPS cycle. The arrows denote dominance relationship between the strategies for \(\lambda =0\); "R" beats "S"; "S" beats "P"; "P" beats "R". The second terms proportional to \(\lambda\) are introduced between sites connected by the double line. (b) The real- and imaginary-parts of eigenvalues as functions of \(\lambda\). An eigenvalue of \(A(\lambda )\) is zero for an arbitrary \(\lambda\) (see horizontal dashed lines). (c) The \(\mathbb {Z}_2\)-invariant characterizing the EP. The dashed vertical lines in panels (b,c) denote \(\lambda =\lambda _{\mathrm {c}}\). Full size image Dynamical properties and the EP Suppose that a large number of players repeat the above game whose dynamics is described by the replicator equation [see Eq. (4)]. In this case, the EP governs the dynamical behaviors when the population density slightly deviates from a fixed point; the population density shows an oscillatory behavior for \(0\le \lambda <\lambda _{\mathrm {c}}\), while such a behavior disappears for \(\lambda _{\mathrm {c}} <\lambda\) due to the emergence of the EP. Firstly, we linearize the replicator equation. When a large number of players repeat the game, the time-evolution is described by the replicator equation $$\begin{aligned} \partial _t \varvec{e}^T_I \cdot \varvec{x}= \ varvec{e}^T_I \cdot \varvec{x} \left( \varvec{e}^T_I A \varvec{x} -\ varvec{x}^TA\varvec{x} \right) , \end{aligned}$$ (4) where \(\varvec{x}\) denotes the population density (\(\sum _I x_I=1 \)), and \(\varvec{e}_I\) is a unit vector whose I-th element is unity. The second term vanishes when the payoff matrix A is anti-symmetric. However, in order to access the non-Hermitian topology, the second term is inevitable which makes the argument in Ref.^22 unavailable. Nevertheless, we can still obtain the following linearized equation $$\begin{aligned} \partial _t \delta \varvec{x}= \frac{1}{N_0} A \ delta \varvec{x}, \end{aligned}$$ (5) which is mathematically equivalent to the Schrodinger equation. This mathematical equivalence reveals that the dynamics of the RPS cycle (i.e., a classical system) can be understood in terms of quantum physics. Here, \(\delta \varvec{x}\) is defined as \(\delta \varvec {x}=\varvec{x}-\varvec{c}\) with \(\varvec{c}=(1,1,\ldots ,1)^T/N_0\) , and \(N_0\) denotes dimensions of the matrix A. Key ingredients are the following relations: $$\begin{aligned} \sum _{J}A_{IJ}=0 \quad&\mathrm {and}&\quad \sum _ {J}A_{JI}=0, \end{aligned}$$ (6) for an arbitrary I. This equation guarantees that \(\varvec{c}\) denotes a fixed point. Linearizing the replicator equation around this fixed point, we obtain Eq. (5) as detailed in the "Methods". Figure 2 figure 2 Dynamical properties of the RPS game (1). Panels (a,c,e) display the time-evolution of the population \(\varvec{x}\) for \(\lambda =0\), 0.5, and 2, respectively. The data are obtained by employing a fourth order Runge-Kutta method^68 with discretized time \(t_n=n\Delta t_{\ mathrm {RK}}\) with \(\Delta t_{\mathrm {RK}}=0.05\). We set the initial condition as \(\varvec{x}(t=0)=(1-\delta _0,1+\delta _0/2,1+\ delta _0/2)/3\) with \(\delta _0=0.1\). The time-evolution is computed up to \(t_n=120\). Panels (b,d,f) are phase plots for \(\ varvec{x}\). The arrows denote the direction of the dynamics \((\ Delta x_1,\Delta x_2,\Delta x_3)\) which is computed by \(\Delta x_I= x_I (\varvec{e}^T_IA\varvec{x}-\varvec{x}^TA\varvec{x})\Delta t\) with \(\Delta t=0.1\). The black crosses denote the fixed point specified by \(\varvec{c}=(1,1,1)/3\). In panels (b,d,f), the blue line denotes the population density \(\varvec{x}(t)\) plotted in panels (a,c,e), respectively. Full size image The linearized replicator equation indicates that the dynamics is governed by the spectrum of A. In particular, the imaginary-part of the eigenvalues results in oscillatory behaviors. To verify the above statement, we numerically solve the replicator Eq. (4). Figure 2 plots the time-evolution of the population density for several values of \(\lambda\). For \(\lambda =0\), the matrix \(A (\lambda )\) is reduced to the payoff matrix of the standard zero-sum RPS game. In this case, the eigenvalues of A are pure imaginary, which results in the oscillatory behavior (see Fig. 2a). This oscillatory behavior is also observed in the phase plot (see Fig. 2 b), where the orbit forms a closed loop. The above oscillatory behavior is also observed for \(\lambda =0.5\) (see Fig. 2c,d), which is due to the imaginary-part of the eigenvalues. We also note that the deviation \(\delta \varvec{x}\) is enhanced as t increases, which is because the real-part of the eigenvalues is positive (see Fig. 1 b). Further increasing \(\lambda\) induces the EP (see Fig. 1b), and the eigenvalues become real. Correspondingly, the above oscillatory behavior is not observed for \(\lambda =2\) (see Fig. 2e,f), The above numerical data verify that the EP governs the dynamics around the fixed point \(\varvec{c}\); the oscillatory behavior of the RPS cycle disappears as the EP emerges. We note that for \(\lambda < 0\), the fixed point corresponds to evolutionary stable strategy. A similar EP emerges also in this case which governs the dynamics. We close this part with three remarks. Firstly, so far, we have seen the emergence of the EP in the RPS game by changing \(\lambda\). We note that symmetry-protected exceptional rings are also observed by introducing an additional parameter (see Sect. S2 of Supplemental Material [68]), whose topology is also characterized by the \(\mathbb {Z}_2\)-invariant \(\nu\). Secondly, although Refs.^69,70,71,72 discuss topology of interaction networks among strategies (i.e., interaction topology), it differs from the topology discussed in this paper. Thirdly, we note that EPs are also reported for active matters ^73,74. We would like to stress, however, that significance of this paper is to reveal the emergence of the EPs and how they affect the dynamics in systems of the evolutionary game theory which describes the population density of biological systems and human societies. Skin effect in an RPS chain Now, we discuss a one-dimensional system showing the skin effect whose origin is the non-trivial topology characterized by the winding number^48,49. Because of this non-trivial topology, switching from the periodic boundary condition (PBC) to the open boundary condition (OBC) significantly changes the spectrum. Correspondingly, almost of all right eigenstates are localized around the right edge which are called skin modes. Figure 3 figure 3 (a) Sketch of the RPS chain for \(L_x=4\) with an additional site at \(I=9\). The arrows describes payoffs. The terms proportional to \(\ lambda\) are introduced between sites connected by double lines. As shown in this panel, I takes \(I=1,2,3,\ldots\). For the PBC, \ (I+2L_x=I\) holds. Dashed line denotes the unit cell. For \(J=1,2,\ ldots\), \(R_{2I-1}= R_{2I}\) holds. (b) Spectrum of the RPS chain for \(\lambda =-0.5\). The data colored with blue are obtained under the PBC. Here, the shade of color denotes \(k/\pi\). The data denoted by red dots are obtained for \(L_x=15\) and the OBC. (c) Amplitude of the right eigenvectors \(\phi ^{\mathrm {R}}_{In}\) (\(n=1,2,\ldots ,2L_x\)) as functions of I for \(L_x=15\). Red (gray) lines denote data for the OBC (PBC). Full size image The above skin effect can be observed in an RPS chain illustrated in Fig. 3a. Applying the Fourier transformation, the payoff matrix is written as $$\begin{aligned} A(k)= & {} \left( \begin{array}{cc} 0 &{} 1-e^{-ik} \\ -1+e^{ik} &{} -2i\sin k \end{array} \right) + \lambda \left( \ begin{array}{cc} -2 &{} 1+e^{-ik} \\ 1+e^{ik} &{} -2 \end{array} \ right) , \end{aligned}$$ (7) under the PBC. Here, \(\varvec{x}_{k}\) is defined as \(\varvec{x}^T_ {k}=(x_{k\mathrm {A}},x_{k\mathrm {B}})\) with \(x_{k\alpha }=\frac {1}{L_x} \sum _{R_I} e^{ik R_I} x_{R_I\alpha }\) and \(\alpha =\ mathrm {A},\mathrm {B}\). Sets of \(R_I\) and \(\alpha _{I}\) are specified by I (\(x_I=x_{R_I\alpha _I}\)). For the explicit form of the payoff matrix in the real-space, see Sect. S3 of Supplemental Material [68]. Figure 3b plots the spectrum of the payoff matrix for \(\lambda =-0.5 \). When the PBC is imposed, eigenvalues form a loop structure as denoted by blue lines in Fig. 3b. Accordingly, the winding number^31, 48,49 defined as $$\begin{aligned} W= & {} \oint \frac{dk}{2\pi i} \partial _k \log \ mathrm {det}[A(k)-\epsilon _{\mathrm {ref}}\mathbbm{1} ], \end {aligned}$$ (8) takes \(-1\) for \(\epsilon _{\mathrm {ref}}=1\), which implies the skin effect. Indeed, imposing the OBC significantly changes the spectrum (see red dots in Fig. 3b). Correspondingly, almost all of the right eigenvectors are localized around the edges, meaning the emergence of skin modes (see Fig. 3c). Figure 4 figure 4 Time-evolution of the population density \(\delta \varvec{x}(t)=\ varvec{x}(t)-\varvec{c}\) for \(\lambda =-0.5\) and \(L_x=15\). The horizontal axis denotes I. (a-f) The time-evolution under the PBC [OBC']. The data in panels (a,e) are multiplied by 5. The data are obtained by employing the fourth order Runge-Kutta method^68 with discretized time \(t_n=n\Delta t_{\mathrm {RK}}\) with \(\Delta t_{\ mathrm {RK}}=0.05\). We set the initial condition so that \(N_0 x_{I} (0)\) takes \(1+\delta _0\) and \(1-\delta _0\) with \(\delta _0=0.1 \) for \(I=25\) and \(I=26\), respectively. For the other I, \(N_0 x_ {I}(0)\) takes 1. Here, \(N_0\) is chosen so that \(\sum _Ix_I(0)=1\) holds. Full size image The above data (Fig. 3b,c) indicate that the skin effect is observed in the RPS chain. Here, we note that the spectrum and skin modes are almost unchanged under the following perturbation to the edges (see Sect. S3 of Supplemental Material [68]): attaching site at \(I=2L_x+1 \) and tuning diagonal elements to \(A_{II}=\lambda\) only for \(I= 1,2L_x+1\) so that Eq. (6) holds (see Fig. 3a). We refer to this boundary condition as OBC'. The above skin effect results in striking dynamical properties: the directive propagation of the population density in the bulk and the enhancement of the population density only around the right edge. The directive propagation is observed by imposing the PBC. Figure 4a-d indicate that the propagation only to the right direction is enhanced. This phenomenon can be understood by noting the following facts as well as linearized approximation: (1) as t increases, each mode is enhanced corresponding to \(\mathrm {Re}\epsilon _n\); (2) the group velocity of each mode is proportional to \(-\partial _k \ mathrm {Im} \epsilon _n\). Because of the loop structure resulting in \(W=-1\), the modes propagating to the right are more enhanced than the ones propagating to the left. The enhancement of the population density at the right edge is observed by imposing OBC'. We note that Eq. (6) is satisfied for OBC' while it is not for OBC. As shown in Fig. 4e-h, the population density around the right edge is enhanced as t increases. This phenomenon is due to the skin modes whose eigenvalues satisfy \(\ mathrm {Re}\epsilon _n>0\). The above dynamical behaviors are unique to non-Hermitian systems; turning off \(\lambda\), \(iA(\lambda =0)\) becomes Hermitian and the above behaviors disappear (see Sect. S3 of Supplemental Material [68]). We close this part with two remarks. Firstly, we note that a similar behaviors can be observed for another RPS chain, implying ubiquity of the skin effect (see Sect. S4 of Supplemental Material [68]). Secondly, Ref.^21 has analyzed an RPS chain whose payoff matrix is Hermitian up to the imaginary unit i. We would like to stress, however, that our aim is to observe non-Hermitian topological phenomena which are not accessible with the model in Ref.^21. In this non-Hermitian case, we need to take into account the second term of Eq. (4). Discussion We have proposed a new platform of non-Hermitian topology in an interdisciplinary field, i.e., RPS cycles described by the evolutionary game theory. Specifically, by analyzing the payoff matrix, we have demonstrated the emergence of the EP and the skin effect which are representative non-Hermitian topological phenomena. In addition, our linearized replicator equation has revealed that the EP governs the dynamics of the population density around the fixed point. Furthermore, we have discovered the striking dynamical phenomenon in the RPS chain: the directive propagation of the population density in the bulk and the enhancement of the population density only around the right edge which are induced by the skin effect. Our results pose several future directions which we discuss below. The experimental observation is one of the central issues. In particular, our results provide a first step toward the observation of the non-Hermitian topology beyond natural science because the game theory describes a wide variety of systems from biological systems^25 ,75 to human societies^26,27,28,76,77; for instance, dynamics of human cooperation has been discussed in Refs.^76,77. The experimental observation in such a system is considered to be a significant step to the application of topological phenomena beyond natural science. In addition, as is the case of equatorial wave^18, our result may provide a novel perspective of well-known phenomena for biological systems such as bacteria^25 and side-blotched lizards^75; dynamics of such systems may be understood in terms of exceptional points. We also note that the topological classification for systems described by the game theory is also a crucial issue to be addressed. Methods Derivation of the linearized replicator equation Here, we derive Eq. (5). We start with noting Eq. (6) results in the relations \(A\varvec{c}=0\) and \(\varvec{c}^TA=0\) which mean that \ (\varvec{c}=(1,1,1,\ldots ,1)^T/N_0\) is a fixed point. By making use of the above relations we have $$\begin{aligned}&\partial _t \varvec{e}^T_I \cdot \delta \varvec{x} \nonumber \\&\quad = \varvec{e}^T_I \cdot (\varvec{c}+\delta \varvec {x}) \left[ \varvec{e}^T_I A (\varvec{c}+\delta \varvec{x}) -(\varvec {c}+\delta \varvec{x})^T A (\varvec{c}+\delta \varvec{x}) \right] \ nonumber \\&\quad = \varvec{e}^T_I \cdot (\varvec{c}+\delta \varvec {x}) \left( \varvec{e}^T_I A \delta \varvec{x} -\delta \varvec{x}^TA \delta \varvec{x} \right) \nonumber \\&\quad \sim (\varvec{e}^T_I \ cdot \varvec{c}) \varvec{e}^T_I A \delta \varvec{x}. \end{aligned}$$ (9) From the second to the third line, we have used the relations \(A\ varvec{c}=0\) and \(\varvec{c}A=0\). In the last line, we have discarded the second and third order terms of \(\delta \varvec{x}\). Because \((\varvec{e}^T_I \cdot \varvec{c})=1/N_0\) for an arbitrary I, we have $$\begin{aligned} \partial _t \delta \varvec{e}^T_I\cdot \varvec{x}\ sim & {} \frac{1}{N_0} \varvec{e}^T_I A \delta \varvec{x}, \end {aligned}$$ (10) which is equivalent to Eq. (5). References 1. Kane, C. L. & Mele, E. J. \({Z}_{2}\) topological order and the quantum spin hall effect. Phys. Rev. Lett. 95, 146802 (2005). ADS CAS PubMed Google Scholar 2. Bernevig, B. A., Hughes, T. L. & Zhang, S.-C. Quantum spin hall effect and topological phase transition in hgte quantum wells. Science 314, 1757-1761 (2006). ADS CAS PubMed Google Scholar 3. Konig, M. et al. Quantum spin hall insulator state in hgte quantum wells. Science 318, 766-770 (2007). ADS PubMed Google Scholar 4. Hasan, M. Z. & Kane, C. L. Colloquium. Rev. Mod. Phys. 82, 3045-3067 (2010). ADS CAS Google Scholar 5. Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057-1110 (2011). ADS CAS Google Scholar 6. Klitzing, K. V., Dorda, G. & Pepper, M. New method for high-accuracy determination of the fine-structure constant based on quantized hall resistance. Phys. Rev. Lett. 45, 494-497 (1980). ADS Google Scholar 7. Thouless, D. J., Kohmoto, M., Nightingale, M. P. & Nijs, M. D. Quantized hall conductance in a two-dimensional periodic potential. Phys. Rev. Lett. 49, 405-408 (1982). ADS CAS Google Scholar 8. Halperin, B. I. Quantized hall conductance, current-carrying edge states, and the existence of extended states in a two-dimensional disordered potential. Phys. Rev. B 25, 2185-2190 (1982). ADS Google Scholar 9. Hatsugai, Y. Chern number and edge states in the integer quantum hall effect. Phys. Rev. Lett. 71, 3697-3700 (1993). ADS MathSciNet CAS PubMed MATH Google Scholar 10. Haldane, F. D. M. & Raghu, S. Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry. Phys. Rev. Lett. 100, 013904 (2008). ADS CAS PubMed Google Scholar 11. Wang, Z., Chong, Y., Joannopoulos, J. D. & Soljacic, M. Observation of unidirectional backscattering-immune topological electromagnetic states. Nature 461, 772 EP (2009). ADS Google Scholar 12. Fu, J.-X., Liu, R.-J. & Li, Z.-Y. Robust one-way modes in gyromagnetic photonic crystal waveguides with different interfaces. Appl. Phys. Lett. 97, 041112 (2010). ADS Google Scholar 13. Ozawa, T. et al. Topological photonics. Rev. Mod. Phys. 91, 015006 (2019). ADS MathSciNet CAS Google Scholar 14. Prodan, E. & Prodan, C. Topological phonon modes and their role in dynamic instability of microtubules. Phys. Rev. Lett. 103, 248101 (2009). ADS PubMed Google Scholar 15. Kane, C. L. & Lubensky, T. C. Topological boundary modes in isostatic lattices. Nat. Phys. 10, 39 EP (2013). Google Scholar 16. Kariyado, T. & Hatsugai, Y. Manipulation of dirac cones in mechanical graphene. Sci. Rep. 5, 18107 (2015). ADS CAS PubMed PubMed Central Google Scholar 17. Susstrunk, R. & Huber, S. D. Observation of phononic helical edge states in a mechanical topological insulator. Science 349, 47-50 (2015). ADS PubMed Google Scholar 18. Delplace, P., Marston, J. B. & Venaille, A. Topological origin of equatorial waves. Science 358, 1075-1077 (2017). ADS MathSciNet CAS PubMed MATH Google Scholar 19. Sone, K. & Ashida, Y. Anomalous topological active matter. Phys. Rev. Lett. 123, 205502 (2019). ADS CAS PubMed Google Scholar 20. Yoshida, T. & Hatsugai, Y. Bulk-edge correspondence of classical diffusion phenomena. Sci. Rep. 11, 888 (2021). CAS PubMed PubMed Central Google Scholar 21. Knebel, J., Geiger, P. M. & Frey, E. Topological phase transition in coupled rock-paper-scissors cycles. Phys. Rev. Lett. 125, 258301 (2020). ADS MathSciNet CAS PubMed Google Scholar 22. Yoshida, T., Mizoguchi, T. & Hatsugai, Y. Chiral edge modes in evolutionary game theory: A kagome network of rock-paper-scissors cycles. Phys. Rev. E 104, 025003 (2021). ADS MathSciNet CAS PubMed Google Scholar 23. Weibull, J. W. Evolutionary Game Theory (MIT press, 1997). 24. Sigmund, K. Evolutionary Game Dynamics Vol. 69 (American Mathematical Society, 2011). 25. Kirkup, B. C. & Riley, M. A. Antibiotic-mediated antagonism leads to a bacterial game of rock-pape-scissors in vivo. Nature 428, 412-414 (2004). ADS CAS PubMed Google Scholar 26. Rosas, A. Evolutionary game theory meets social science: Is there a unifying rule for human cooperation?. J. Theor. Biol. 264, 450-456 (2010). ADS MathSciNet PubMed MATH Google Scholar 27. Wang, Z., Kokubo, S., Tanimoto, J., Fukuda, E. & Shigaki, K. Insight into the so-called spatial reciprocity. Phys. Rev. E 88, 042145 (2013). ADS Google Scholar 28. Perc, M. et al. Statistical physics of human cooperation. Phys. Rep. 687, 1-51 (2017). ADS MathSciNet MATH Google Scholar 29. Hatano, N. & Nelson, D. R. Localization transitions in non-hermitian quantum mechanics. Phys. Rev. Lett. 77, 570-573 (1996). ADS CAS PubMed Google Scholar 30. Bender, C. M. & Boettcher, S. Real spectra in non-hermitian hamiltonians having \(P{\rm T}\) symmetry. Phys. Rev. Lett. 80, 5243-5246 (1998). ADS MathSciNet CAS MATH Google Scholar 31. Gong, Z. et al. Topological phases of non-hermitian systems. Phys. Rev. X 8, 031079 (2018). CAS Google Scholar 32. Kawabata, K., Shiozaki, K., Ueda, M. & Sato, M. Symmetry and topology in non-hermitian physics. Phys. Rev. X 9, 041015 (2019). CAS Google Scholar 33. Bergholtz, E. J., Budich, J. C. & Kunst, F. K. Exceptional topology of non-hermitian systems. Rev. Mod. Phys. 93, 015005 (2021). ADS MathSciNet Google Scholar 34. Ashida, Y., Gong, Z. & Ueda, M. Non-hermitian physics. Adv. Phys. 69, 249-435 (2020). ADS Google Scholar 35. Rotter, I. A non-hermitian hamilton operator and the physics of open quantum systems. J. Phys. A Math. Theor. 42, 153001 (2009). ADS MathSciNet MATH Google Scholar 36. Shen, H., Zhen, B. & Fu, L. Topological band theory for non-hermitian hamiltonians. Phys. Rev. Lett. 120, 146402 (2018). ADS MathSciNet CAS PubMed Google Scholar 37. Kozii, V. & Fu, L. Non-hermitian topological theory of finite-lifetime quasiparticles: Prediction of bulk fermi arc due to exceptional point. Preprint arXiv:1708.05841 (2017). 38. Yoshida, T., Peters, R. & Kawakami, N. Non-hermitian perspective of the band structure in heavy-fermion systems. Phys. Rev. B 98, 035141 (2018). ADS CAS Google Scholar 39. Budich, J. C., Carlstrom, J., Kunst, F. K. & Bergholtz, E. J. Symmetry-protected nodal phases in non-hermitian systems. Phys. Rev. B 99, 041406 (2019). ADS CAS Google Scholar 40. Okugawa, R. & Yokoyama, T. Topological exceptional surfaces in non-hermitian systems with parity-time and parity-particle-hole symmetries. Phys. Rev. B 99, 041202 (2019). ADS CAS Google Scholar 41. Yoshida, T., Peters, R., Kawakami, N. & Hatsugai, Y. Symmetry-protected exceptional rings in two-dimensional correlated systems with chiral symmetry. Phys. Rev. B 99, 121101 (2019). ADS CAS Google Scholar 42. Zhou, H., Lee, J. Y., Liu, S. & Zhen, B. Exceptional surfaces in pt-symmetric non-hermitian photonic systems. Optica 6, 190-193 (2019). ADS CAS Google Scholar 43. Kimura, K., Yoshida, T. & Kawakami, N. Chiral-symmetry protected exceptional torus in correlated nodal-line semimetals. Phys. Rev. B 100, 115124 (2019). ADS CAS Google Scholar 44. Yoshida, T. & Hatsugai, Y. Exceptional rings protected by emergent symmetry for mechanical systems. Phys. Rev. B 100, 054109 (2019). ADS CAS Google Scholar 45. Delplace, P., Yoshida, T. & Hatsugai, Y. Symmetry-protected higher-order exceptional points and their topological characterization. Preprint arXiv:2103.08232 (2021). 46. Yao, S. & Wang, Z. Edge states and topological invariants of non-hermitian systems. Phys. Rev. Lett. 121, 086803 (2018). ADS CAS PubMed Google Scholar 47. Lee, C. H. & Thomale, R. Anatomy of skin modes and topology in non-hermitian systems. Phys. Rev. B 99, 201103 (2019). ADS CAS Google Scholar 48. Zhang, K., Yang, Z. & Fang, C. Correspondence between winding numbers and skin modes in non-hermitian systems. Phys. Rev. Lett. 125, 126402 (2020). ADS MathSciNet CAS PubMed Google Scholar 49. Okuma, N., Kawabata, K., Shiozaki, K. & Sato, M. Topological origin of non-hermitian skin effects. Phys. Rev. Lett. 124, 086801 (2020). ADS MathSciNet CAS PubMed Google Scholar 50. Lee, J. Y., Ahn, J., Zhou, H. & Vishwanath, A. Topological correspondence between hermitian and non-hermitian systems: Anomalous dynamics. Phys. Rev. Lett. 123, 206404 (2019). ADS MathSciNet CAS PubMed Google Scholar 51. Borgnia, D. S., Kruchkov, A. J. & Slager, R.-J. Non-hermitian boundary modes and topology. Phys. Rev. Lett. 124, 056802 (2020). ADS MathSciNet CAS PubMed Google Scholar 52. Yoshida, T., Mizoguchi, T. & Hatsugai, Y. Mirror skin effect and its electric circuit simulation. Phys. Rev. Res. 2, 022062 (2020). CAS Google Scholar 53. Guo, A. et al. Observation of \(\cal{P}\cal{T}\)-symmetry breaking in complex optical potentials. Phys. Rev. Lett. 103, 093902 (2009). ADS CAS PubMed Google Scholar 54. Zhen, B. et al. Spawning rings of exceptional points out of dirac cones. Nature 525, 354 EP (2015). ADS Google Scholar 55. Hassan, A. U., Zhen, B., Soljacic, M., Khajavikhan, M. & Christodoulides, D. N. Dynamically encircling exceptional points: Exact evolution and polarization state conversion. Phys. Rev. Lett. 118, 093002 (2017). ADS PubMed Google Scholar 56. Zhou, H. et al. Observation of bulk fermi arc and polarization half charge from paired exceptional points. Science 359, 1009-1012 (2018). ADS MathSciNet CAS PubMed MATH Google Scholar 57. Zyuzin, A. A. & Zyuzin, A. Y. Flat band in disorder-driven non-hermitian weyl semimetals. Phys. Rev. B 97, 041203 (2018). ADS CAS Google Scholar 58. Papaj, M., Isobe, H. & Fu, L. Nodal arc of disordered dirac fermions and non-hermitian band theory. Phys. Rev. B 99, 201107 (2019). ADS CAS Google Scholar 59. Shen, H. & Fu, L. Quantum oscillation from in-gap states and a non-hermitian landau level problem. Phys. Rev. Lett. 121, 026403 (2018). ADS CAS PubMed Google Scholar 60. Matsushita, T., Nagai, Y. & Fujimoto, S. Disorder-induced exceptional and hybrid point rings in weyl/dirac semimetals. Phys. Rev. B 100, 245205 (2019). ADS CAS Google Scholar 61. Hofmann, T. et al. Reciprocal skin effect and its realization in a topolectrical circuit. Phys. Rev. Res. 2, 023265 (2020). CAS Google Scholar 62. Helbig, T. et al. Generalized bulk oundary correspondence in non-hermitian topolectrical circuits. Nat. Phys. 16, 747-750 (2020). CAS Google Scholar 63. Matsushita, T., Nagai, Y. & Fujimoto, S. Spectrum collapse of disordered dirac landau levels as topological non-hermitian physics. J. Phys. Soc. Jpn. 90, 074703 (2021). ADS Google Scholar 64. Tainaka, K. Paradoxical effect in a three-candidate voter model. Phys. Lett. A 176, 303-306 (1993). ADS Google Scholar 65. Juul, J., Sneppen, K. & Mathiesen, J. Clonal selection prevents tragedy of the commons when neighbors compete in a rock-paper-scissors game. Phys. Rev. E 85, 061924 (2012). ADS Google Scholar 66. Szolnoki, A., Vukov, J. & Perc, M. C. V. From pairwise to group interactions in games of cyclic dominance. Phys. Rev. E 89, 062125 (2014). ADS Google Scholar 67. Szolnoki, A. & Perc, M. Biodiversity in models of cyclic dominance is preserved by heterogeneity in site-specific invasion rates. Sci. Rep. 6, 38608 (2016). ADS CAS PubMed PubMed Central Google Scholar 68. Suli, E. & Mayers, D. F. An Introduction to Numerical Analysis (Cambridge University Press, 2003). 69. Szabo, G. & Fath, G. Evolutionary games on graphs. Phys. Rep. 446 , 97-216 (2007). ADS MathSciNet Google Scholar 70. Szolnoki, A. et al. Cyclic dominance in evolutionary games: A review. J. R. Soc. Interface 11, 20140735 (2014). PubMed PubMed Central Google Scholar 71. Dobramysl, U., Mobilia, M., Pleimling, M. & Tauber, U. C. Stochastic population dynamics in spatially extended predator-prey systems. J. Phys. A: Math. Theor. 51, 063001 (2018). ADS MathSciNet MATH Google Scholar 72. Szolnoki, A., Oliveira, B. F. D. & Bazeia, D. Pattern formations driven by cyclic interactions: A brief review of recent developments. Europhys. Lett. 131, 68001 (2020). ADS CAS Google Scholar 73. Sone, K., Ashida, Y. & Sagawa, T. Exceptional non-hermitian topological edge mode and its application to active matter. Nat. Commun. 11, 5745 (2020). ADS CAS PubMed PubMed Central Google Scholar 74. Tang, E., Agudo-Canalejo, J. & Golestanian, R. Topology protects chiral edge currents in stochastic systems. Phys. Rev. X 11, 031015 (2021). CAS Google Scholar 75. Sinervo, B. & Lively, C. M. The rock-paper-scissors game and the evolution of alternative male strategies. Nature 380, 240-243 (1996). ADS CAS Google Scholar 76. Szolnoki, A., Perc, M. C. V. & Szabo, G. Defense mechanisms of empathetic players in the spatial ultimatum game. Phys. Rev. Lett. 109, 078701 (2012). ADS PubMed Google Scholar 77. Szolnoki, A. & Perc, M. C. V. Correlation of positive and negative reciprocity fails to confer an evolutionary advantage: Phase transitions to elementary strategies. Phys. Rev. X 3, 041021 (2013). Google Scholar Download references Acknowledgements This work is supported by JSPS Grant-in-Aid for Scientific Research on Innovative Areas "Discrete Geometric Analysis for Materials Design": Grant No. JP20H04627. This work is also supported by JSPS KAKENHI Grants Nos. JP17H06138, JP20K14371, and JP21K13850. Author information Affiliations 1. Department of Physics, University of Tsukuba, Tsukuba, Ibaraki, 305-8571, Japan Tsuneya Yoshida, Tomonari Mizoguchi & Yasuhiro Hatsugai Authors 1. Tsuneya Yoshida View author publications You can also search for this author in PubMed Google Scholar 2. Tomonari Mizoguchi View author publications You can also search for this author in PubMed Google Scholar 3. Yasuhiro Hatsugai View author publications You can also search for this author in PubMed Google Scholar Contributions T.Y., T.M., and Y.H. planned the project. T.Y. performed the calculations. T.Y., and T.M., and Y.H. were involved in the discussion of the materials and the preparation of the manuscript. Corresponding author Correspondence to Tsuneya Yoshida. Ethics declarations Competing interests The authors declare no competing interests. Additional information Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Information Supplementary Information. Rights and permissions Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and Permissions About this article Verify currency and authenticity via CrossMark Cite this article Yoshida, T., Mizoguchi, T. & Hatsugai, Y. Non-Hermitian topology in rock-paper-scissors games. Sci Rep 12, 560 (2022). https://doi.org/ 10.1038/s41598-021-04178-8 Download citation * Received: 07 October 2021 * Accepted: 16 December 2021 * Published: 12 January 2022 * DOI: https://doi.org/10.1038/s41598-021-04178-8 Share this article Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative Comments By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Download PDF Advertisement Advertisement Explore content * Research articles * News & Comment * Collections * Subjects * Follow us on Facebook * Follow us on Twitter * Sign up for alerts * RSS feed About the journal * About Scientific Reports * Journal policies * Guide to referees * Contact * Calls for Papers * Editor's Choice * Guest Edited Collections * Scientific Reports Top 100 2019 * Scientific Reports Top 100 2018 * Scientific Reports Top 10 2018 * Scientific Reports Top 100 2017 * Editorial Board Highlights * Author Highlights * Announcements * 10th Anniversary Editorial Board Interviews * Scientific Reports Top 100 2020 Publish with us * For authors * Submit manuscript Search Search articles by subject, keyword or author [ ] Show results from [All journals] Search Advanced search Quick links * Explore articles by subject * Find a job * Guide to authors * Editorial policies Scientific Reports (Sci Rep) ISSN 2045-2322 (online) nature.com sitemap Nature portfolio * About us * Press releases * Press office * Contact us * * * Discover content * Journals A-Z * Articles by subject * Nano * Protocol Exchange * Nature Index Publishing policies * Nature portfolio policies * Open access Author & Researcher services * Reprints & permissions * Research data * Language editing * Scientific editing * Nature Masterclasses * Nature Research Academies Libraries & institutions * Librarian service & tools * Librarian portal * Open research * Recommend to library Advertising & partnerships * Advertising * Partnerships & Services * Media kits * Branded content Career development * Nature Careers * Nature Conferences * Nature events Regional websites * Nature Africa * Nature China * Nature India * Nature Italy * Nature Japan * Nature Korea * Nature Middle East Legal & Privacy * Privacy Policy * Use of cookies * Manage cookies/Do not sell my data * Legal notice * Accessibility statement * Terms & Conditions * California Privacy Statement Springer Nature (c) 2022 Springer Nature Limited Close Nature Briefing Sign up for the Nature Briefing newsletter -- what matters in science, free to your inbox daily. Email address [ ] Sign up [ ] I agree my information will be processed in accordance with the Nature and Springer Nature Limited Privacy Policy. Close Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing * *