Quantum Dynamical Tunneling Breaks Classical Conserved Quantities
Abstract
We discover that quantum dynamical tunneling, occurring between phase space regions in a classically forbidden way, can break conserved quantities in pseudointegrable systems. We rigorously prove that a conserved quantity in a class of typical pseudointegrable systems can be broken quantum mechanically. We then numerically compute the uncertainties of this broken conserved quantity, which remain non-zero for up to eigenstates and exhibit universal distributions similar to energy level statistics. Furthermore, all the eigenstates with large uncertainties show the superpositions of regular orbits with different values of the conserved quantity, showing definitive manifestation of dynamical tunneling. A random matrix model is constructed to successfully reproduce the level statistics in pseudointegrable systems.
I Introduction
Dynamical tunneling is a fundamental quantum phenomenon which refers to classically forbidden transitions between phase-space regions, even without an energy barrier in the middle [1]. It has been well stuided in both classically integrable and mixed systems (where integrable and chaotic phase-space regions coexist) within various contexts [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 2]. In such systems, dynamical tunneling was not observed to impact the overall integrability associated with the conserved quantities governing the entire system. In this paper, we use pseudointegrable systems to show that dynamical tunneling can break classical conserved quantities other than energy, giving rise to non-integrable behaviors. The pseudointegrable systems were introduced by Richens and Berry in [14]. They are defined as classical Hamiltonian systems with equal degrees of freedom and conserved quantities, where the phase trajectories are restricted to an invariant surface featuring a multi-handled sphere with genus larger than , in contrast to the genus- torus found in integrable systems. Typical examples of such systems are rational polygon billiards (see Fig. 1(a)), whose conserved quantity other than energy is given by
(1) |
where is the angle between the momentum and the horizontal axis and is the least common multiple of for all vertices with rational angles [15]. The conservation arises from the fact that the directions of each classical trajectory follow the orbits of a dihedral group .
The dynamical properties of pseudointegrable systems have been studied over decades both classically and quantum mechanically but from very different perspectives. In classical mechanics, people are concerned about how the unremovable singularities on split the beams of trajectories [14, 15]. In the quantum counterpart, the exploration has primarily revolved around their energy level statistics [14, 16, 17, 18, 19]. Following this line, Bogomolny identified a subset of them characterized by the semi-Poisson distribution [20, 21, 22, 21, 24, 25, 23], which had been proposed as a universal intermediate bridging the gap between Poisson statistics and Wigner-Dyson statistics [26, 27]. Pseudointegrable systems present two advantages for our study: i) Singularities in phase space were suggested as an important element to provoke dynamical tunneling [3]; ii) Absence of a classical chaotic sea makes the non-integrable behaviors more attributable to the effects of dynamical tunneling. Here we focus on the rational right triangle billiards, as depicted in Fig. 1(b). This subclass has all the features of pseudointegrable systems. We analytically prove the absence of conserved quantity in their quantum counterparts. The entire proof unfolds in two steps. First, assuming as the quantization of is conserved, we can expand the eigenstates by plane waves with the directions following the orbits of group . Second, imposing Dirichlet boundary conditions on the sides, we find that they cannot be satisfied, unless the billiards are completely integrable (see Appendix A for details). To quantify the quantum mechanical breaking of the conserved quantity , we calculate the uncertainty of operator for each eigenstate, defined by
(2) |
where denotes the -th eigenstate. We observe that the decreasing trend of with rising energy cannot be described by elementary functions. Instead, we examine the probability density functions (PDFs) of within an energy shell, revealing nontrivial distributions that can be accurately fitted by the Brody distribution [28, 29]. The eigenstates with sufficiently small are collected, and their unfolded level spacing statistics conform to the Poisson distribution. While, the eigenstates with large showcase a superposition with classical regular orbits with different , regardless of whether they are periodic or not. It is the consequence of dynamical tunneling, which results in avoided crossings among regular states. Inspired by the statistics of , we develop a random matrix model for dynamical tunneling and successfully reproduce the level spacing distributions of pseudointegrable systems.
II conserved quantity
Each classical trajectory has directions being the orbits of group , which fold into an interval . This corresponds to a quotient map , functioning as a conserved quantity. For simplicity, we use the form (1). Consequently, the trajectory in phase space can be projected onto . It is simple to prove is conserved along a classical trajectory, i.e., for a certain reflection at one side, for . Additionally, can be expanded as a polynomial in , which can be quantized as the operator , where is the magnitude of the momentum and is the momentum operator along direction.
III Uncertainty
We calculate the quantum right triangle billiards with one interior angle being , , , , , (the angles will be used to represent the corresponding billiards). The genus of is given by the formula . Consequently, each pair of billiards corresponds to with genus-2, 3, 4 . It is noted that the energy level statistics of and has been found to be semi-Poisson distribution [20]. For each billiard, eigenstates are calculated using a hybrid of scaling method and decomposition method [30, 31]. We use the Weyl’s level to indicate the energy of eigenstates[32], which helps compare different systems with a unified energy scale. The direct results of are depicted in Fig. 2(a). It demonstrates that most points cluster around the small values of . However, there are still eigenstates with large values at higher energy levels. Since does not exhibit a clear decreasing trend with , we analyze statistically. If were conserved, we would expect the PDF of uncertainties, denoted as , to be a delta distribution at . However, the statistical results show that follow various nontrivial distributions. This is exemplified by and billiards in Fig. 2(b1) and Fig. 2(b2). Both have right tails but different peaks. With energy rising, they exhibit a trend to converge to delta distribution. It indicates that in the classical limit, the conserved quantity will recover.
To gain further insights, we consider the unfolded uncertainty, defined as , whose domain is extended from to . The PDFs of , denoted as , are shown in Fig. 3(a). All presented here can be well fitted by Brody distributions , with a single parameter . In particular, for and billiards, whose level statistics are semi-Poissonian, closely resemble the Poisson distribution . The differences in cumulative distribution functions (CDFs) from the best-fitted model, denoted by where the superscript represents “optimal fitting”, are shown in the inset, which give the errors less than . This result shows that for genus-3, 4 billiards have a polynomial repulsion at small and a super-exponential tail at large . Roughly, exhibits an elevated with genus increasing, as shown in Fig. 3(a). This may be understood from Diophantine approximations. Consider a right triangle billiard, as depicted in Fig. 1(b), where the ratio is a good Diophantine approximation of some irrational number. Then, is typically large, resulting in a high genus. Since irrational triangle billiards lack the conserved quantity , a with higher genus naturally implies a more chaotic system.
To confirm our expectation that the level statistics of eigenstates with small enough follows a Poisson distribution, we scan the values of to determine a small enough threshold , such that the unfolded level spacing distribution of eigenstates with aligns with the Poisson distribution with minimal fitting error. This can be formulated as where and is the unfolded level spacing. We call the eigenstates with as “pseudoregular” states. The results for billiard are depicted in Fig. 3(b), where the differences in CDFs of unfolded level spacing with different thresholds from the Poisson distribution are calculated, defined as . It shows that the level statistics of pseudoregular states approximates to , implying an absence of level repulsions in the spectrum if is conserved. It is established that when all the eigenstates are involved (), the level spacing statistics for the whole spectrum follows a semi-Poisson distribution, i.e., . By increasing the threshold to , the level statistics occupies an intermediate position between the Poisson and semi-Poisson distribution. This demonstrates that by controlling the uncertainty threshold , the level statistics can be changed between these two distributions. We will see that the eigenstates with are characterized by the superposition of classical trajectories with distinct values of , which is caused by dynamical tunneling.
IV Dynamical Tunneling
The direct approach to investigate dynamical tunneling is mapping the eigenstates into phase space [35, 33, 34]. Fortunately, in pseudointegrable systems, each classical orbit has a specific value of . As a result, each eigenstate that may be a superposition of different classical orbits can be converted into a distribution of different values . This distribution is expressed by the boundary function of , where is the included angle between the outer normal vector of sides and the momentum as shown in Fig. 1(b). This method follows Ref.[36] and its details are given in Appendix B. With , the conserved quantity can be rewritten as , where the absolute value arises at the hypotenuse . In this way, becomes the connection between and . The results of billiard are shown in Fig. 4. For this billiard, there are only two types of classical orbits: periodic orbits (POs) and uniformly distributed orbits (UDOs) that cover the entire billiard table [37]. It has been established that, in rational right triangle billiards, nearly every periodic orbit contains segments that are perpendicular to a specific side. Consequently, almost all POs has , while UDOs have various values of . Quantum mechanically, there exists “superscar” states, named by Bogomolny [22, 38, 39], exhibiting a superposition of spatially parallel POs. cannot be used to distinguish different superscar states as they all equal to . In both Fig. 4.(a) and Fig. 4.(b), exhibit a single peak, illustrating the features of superscar and UDO, respectively. Both of them are pseudoregular states, characterized by uncertainties smaller than . However, in Fig. 4.(c1) and Fig. 4.(d1), orbits with different values of are superposed (one superscar plus one UDO in (c1), two UDOs in (d1)), causing their relatively large . Notably, eigenstates with small wavenumbers may distribute with a considerable width on the values of , leading to an elevated . But the presence of separated peaks here eliminates this concern. The eigestates neighboring to Fig. 4.(c1) and Fig. 4.(d1) are shown in Fig. 4.(c2) and Fig. 4.(d2). And each pair of them shares common peak positions, signifying common tunneled components. Additionally, their individual level splittings are much smaller than the mean level spacing . This is another indication of dynamical tunneling. In Fig. 4.(c2), there emerges sub-peaks at small comparing to Fig. 4.(c1), which may illustrate more resonant orbits. Therefore, the results indicate that dynamical tunneling destroys the conservation of by superposing trajectories with different directions.
V Avoided Crossing
The dynamical tunneling can be modeled by a random matrix Hamiltonian
(3) |
where the diagonals are energy levels of different classical orbits, and the off-diagonal represent tunneling rate between classical orbits with different values of . is real due to the time-reversal symmetry. As different classical orbits are not correlated, it is reasonable to assume that has a Poisson distribution, i.e., .
The numerical results in Fig. 2(b1) and Fig. 3(a) suggest that the distribution of have the same shape as , that is,
(4) |
Once the two distributions and are known, we can compute the level spacing distribution of (3), which is compared with in Fig. 5, demonstrating a perfect fit. Therefore, combing the results in Fig. 4(b) and Fig. 5, we can conclude that dynamical tunneling is responsible for distorting the level spacing distribution from a Poisson distribution to the semi-Poisson distribution. We expect this mechanism applies to all pseudointegrable systems with diverse level spacing distributions.
VI Discussion
The dynamical tunneling observed in pseudointegrable systems shows different effects compared to those in integrable and mixed systems. We break down them into three aspects. i) Conserved quantity. In mixed systems, one can define a “local” conserved quantity for the trajectories within symmetry-related integrable islands. Dynamical tunneling among these integrable islands or from them into the chaotic sea can destroy this “local” conserved quantity [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. However, here, dynamical tunneling has the potential to break the conserved quantity governing the entire systems by superposing classical orbits with different values of (as shown in Fig. 2 and Fig. 4). ii) Level statistics. For mixed systems, dynamical tunneling can only introduce sub-leading corrections to the Berry-Robnik level statistics, since it weakly couples regular and chaotic states, thus increasing small distances between the corresponding levels [40, 4, 6, 7, 5, 8]. The level repulsion in pseudointegrable systems is completely induced by dynamical tunneling (as shown in Fig. 3.(b)). iii) Tunneling rate. In principle, the tunneling rate between two resonant “double-well states” decreases exponentially with , which can be expressed as where is a parameter related to systems and resonant orbits (for a double-well system , is the imaginary part of action over the energy barrier in the middle) [4, 3]. In the presence of more resonant states, the tunneling rates between the original pair can be enhanced or suppressed by several orders of magnitude at certain values of [9, 10, 11, 12, 13]. We formulate this qualitatively as where is the average value over a shell of , is the enhanced (or suppressed) point and is the corresponding strength. In the limit of strong tunneling, the distribution of becomes dense on the shell, allowing us to treat as a random variable centered around . This is consistent with our semi-empirical distribution of tunneling rate (4). However, a rigorous quantitative derivation of tunneling rate remains an open problem.
ACKNOWLEDGMENTS
LQK and BW are supported by National Natural Science Foundation of China (Grant No. 11921005 & 92365202), and Shanghai Municipal Science and Technology Major Project (Grant No.2019SHZDZX01). ZG is supported by The University of Tokyo Excellent Young Researcher Program.
Appendix A A proof of conserved quantity breaking in rational right triangle billiards
Consider a rational right triangle billiard described in Fig. 6. The direction can undergo three kinds of reflections by sides: (reflection at ), (reflection at ) and (reflection at ). These can generate a dihedral group () with dimension , yielding a conserved quantity independent with the energy. can be quantized after expanding to the polynomial of where is the momentum along -direction and is the norm of total momentum.
The conserved quantity originates from the existence of quotient projection mapping each classical trajectory on the surface , where is the unit circle. It means that for a classical trajectory, the possible directions, denoted as , can be related to a reference so that
(A5) |
If the quantized operator remains conserved in quantum mechanics, the eigenstate should also follow the orbits of group , just as each classical trajectory does. In this manner, the eigenstate, denoted as , can be expanded by plane waves with direction ,
(A6) |
where are the superposing parameters. Then, the Dirichlet boundary conditions at , and can be evaluated sequentially:
(A7) | |||
(A8) | |||
(A9) |
where . Suppose the subscript of has a cyclic relation for , Eq. (A8) can be rewritten as
(A10) |
We can prove that Eq. (A7) and Eq. (A10) are equivalent to
(A11) |
We take the -th derivative of the both sides of Eq. (A7) and make equal to zero. In this way, we obtain homogeneous linear equations by defining ,
(A12) |
Its determinant is known as Vandermonde determinant, yielding
(A13) |
If let , we obtain which cannot be confined in the domain . As a result, this determinant is nonzero for . In this way, the solution to Eq. (A12) is a null vector, i.e.,
(A14) |
When , the billiard has directions corresponding to degenerate cases with . We redefine the superposing parameter as , so that . For , the degenerate version of boundary conditions Eq. (A7) and Eq. (A8) become
(A15) | |||
(A16) |
These two relations are equivalent to , yielding . Consequently, all must equal to . For , we have
(A17) | |||
(A18) |
They are equivalent to , yielding . Consequently, all must equal to . So far, we have proved an absence of common eigenstates for the degenerate cases of . Using the same tricks on Eq. (A10), we can obtain . Applying the relation Eq. (A11) on the eigenstates , we get
(A19) |
One can verify is invariant under the action of dihedral group , for . is a irreducible representation of group . We apply the boundary condition Eq. (A9) to Eq. (A19), yielding
(A20) |
We can prove that Eq. (A20) is equivalent to
(A21) |
where . This proof is similar to the proof of Eq. (A11). We take the -th deviatives of both side of Eq. (A15) and make equal to zero. We obtain homogeneous linear equations by defining ,
(A22) |
Its determinant is nonzero, so the solution is , yielding . The relation Eq. (A21) can be used as the quantization condition. We divide it by itself with a different subscript, and obtain
(A23) |
where for . By using the identity
(A24) |
and setting to be , we find that
(A25) |
illustrating the l.h.s. is rational. It is known that if the cosine of a rational angle is a rational number, the only possible cosine values are and (A proof was provided in [41]). They correspond to the right triangle billiards with and , which are the only two complete integrable cases. For others, Eq. (A25) are not satisfied. Therefore, we conclude that in rational right triangle billiards, is not conserved unless the billiards are integrable.
Appendix B Poincaré-Husimi representation
Poincaré-Husimi representation is used to compare the eigenstates with classical trajectories [36, 42]. In this representation, the Poincaré section is typically chosen at the boundary of billiards, yielding the boundary function defined as the normal derivative of eigenfunction, i.e.,
(B26) |
where is a point on the boundary parameterized by the arc-length and denotes the outer normal unit vector at . Conversely, the wave function can be obtained by boundary integral
(B27) |
where is the billiard boundary curve, is the point inside the billiard table and is the 2-dimensional free particle Green’s function. In phase space, the conjugate variable to is , where represents the included angle between the outer normal vector of sides and the velocity of billiard flow, as shown in Fig. 6. The conjugate function is expressed as
(B28) |
where is the wavenumber. It should be emphasized that the relations between and have three kind depending on the sides: (for ), (for ), (for ). Hence, remains the classical conserved quantity. Without ambiguity, we also use the symbol to denote . The values of can be used to distinguish different trajectories. For example, in billiard, there are only two types of classical orbits: periodic orbits (POs) and uniformly distributed orbits (UDOs) that cover the entire billiard table [37]. It has been established that, almost each PO possesses segments perpendicular to side of right triangle. Consequently, almost all POs has , while UDOs have various values of . Quantum mechanically, there exists “superscar” states, exhibiting a superposition of spatially parallel POs [22, 38, 39]. cannot be used to distinguish different superscar states as they all equal to . But, the complete details of trajectories can be obtained by projecting the eigenstate using Husimi functions on the phase space expanded by and . The Husimi function with , is defined as
(B29) |
where is the normalization parameter. is the standard coherent state, i.e.,
(B30) | ||||
where is the length of boundary and the sum of ensures the coherent states being periodic with a period of .
The results for billiard are shown in Fig. 7. And the Poincaré-Husimi representations and the plot of the conjugate boundary function with respect to the conserved quantity are compared to illustrate the trajectories of typical eigenstates with different (the uncertainty of ). The eigenstate in Fig. 7.(a) is a superscar state, superposed by two parallel POs with segments perpendicular to the hypotenuse. The eigenstate in Fig. 7.(b) illustrates a UDO. Both of these are pseudoregular states, each exhibiting a single peak at the corresponding value of and thus possessing the conserved quantity . The eigenstate in Fig. 7.(c1) is a superposition of one superscar (including four parallel POs) and one UDO from the observation on . Consequently, it has two separated peaks on the distribution of , causing an elevated . The eigenstate in Fig. 7.(c2), which is next to the one in Fig. 7.(c1), includes the trajectories in Fig. 7.(c1) plus an additional UDO (this UDO also exists in Fig. 7.(c1) but with a minimal peak). The eigenstates in Fig. 7.(d1) and Fig. 7.(d2) are another pair with the energy splitting much smaller than the mean level spacing and exhibit a dynamical tunneling between two UDOs regular state.
References
- [1] M. J. Davis and E. J. Heller, Quantum dynamical tunneling in bound states, J. Chem. Phys. 75, 246 (1981).
- [2] N. T. Maitra and E. J. Heller, Semiclassical perturbation approach to quantum reflection, Phys. Rev. A 54, 4763 (1997).
- [3] S. Keshavamurthy and P. Schlagheck, Dynamical Tunneling: Theory and Experiment (CRC Press, Boca Raton, FL, 2011).
- [4] S. Tomsovic and D. Ullmo, Chaos-assisted tunneling, Phys. Rev. E 50, 145 (1994).
- [5] V. A. Podolskiy and E. E. Narimanov, Semiclassical Description of Chaos-Assisted Tunneling, Phys. Rev. Lett. 91, 263601 (2003).
- [6] A. Bäcker, R. Ketzmerick, S. Löck and L. Schilling, Regular-to-Chaotic Tunneling Rates Using a Fictitious Integrable System, Phys. Rev. Lett. 100, 104101 (2008).
- [7] A. Bäcker, R. Ketzmerick, S. Löck and N. Mertig, Fractional-Power-Law Level Statistics Due to Dynamical Tunneling, Phys. Rev. Lett. 106, 024101 (2011).
- [8] A. Relaño, Chaos-assisted tunneling and spectral fluctuations in the order-chaos transition, Phys. Rev. Lett. 100, 224101 (2008).
- [9] D. A. Steck, W. H. Oskay and M. G. Raizen, Observation of chaos-assisted tunneling between islands of stability, Science 293, 274 (2001).
- [10] S. D. Frischat and E. Doron, Dynamical tunneling in mixed systems, Phys. Rev. E 57, 1421 (1998).
- [11] W. A. Lin and L. E. Ballentine, Quantum tunneling and chaos in a driven anharmonic oscillator, Phys. Rev. Lett. 65, 2927 (1990).
- [12] E. Doron and S. D. Frischat, Semiclassical Description of Tunneling in Mixed Systems: Case of the Annular Billiard, Phys. Rev. Lett. 75, 3661 (1995).
- [13] E. Doron and S. D. Frischat, Resonance-Assisted Tunneling in Near-Integrable Systems, Phys. Rev. Lett. 87, 064101 (2001).
- [14] P. J. Richens and M. V. Berry, Pseudointegrable systems in classical and quantum mechanics, Physica D 2, 495 (1981).
- [15] E. Gutkin, Billiards in polygons, Physica D 19, 311 (1986).
- [16] T. Cheon, T. D. Cohen, Quantum level statistics of pseudointegrable billiards, Phys. Rev. Lett. 62, 2769 (1989).
- [17] A. Shudo, Y. Shimizu, P. Šeba, J. Stein, H. J. Stöckmann and K. Życzkowski, Statistical properties of spectra of pseudointegrable systems, Phys. Rev. E 49, 3748 (1994).
- [18] T. Shigehara, Conditions for the ayyearance of wave chaos in quantum singular systems with a pointlike scatterer, Phys. Rev. E 50, 4357 (1994).
- [19] J. Wiersig, Spectral properties of quantized barrier billiards, Phys. Rev. E 65, 046217 (2002).
- [20] E. B. Bogomolny, U. Gerland, C. Schmit, Models of intermediate spectral statistics, Phys. Rev. E 59, R1315 (1999).
- [21] E. B. Bogomolny, O. Giraud, C. Schmit, Periodic orbits contribution to the 2-point correlation form factor for pseudo-integrable systems, Comm. Math. Phys. 222, 327 (2001).
- [22] E. B. Bogomolny and C. Schmit, Structure of wave functions of pseudointegrable billiards, Phys. Rev. Lett. 92, 244102 (2004).
- [23] E. B. Bogomolny, Barrier billiard and random matrices, J. Phys. A 55, 024001 (2021).
- [24] E. B. Bogomolny and C. Schmit, Spectral statistics of a quantum interval-exchange map, Phys. Rev. Lett 93, 254102 (2004).
- [25] E. B. Bogomolny, R. Dubertrand, and C. Schmit, Spectral statistics of a pseudo-integrable map: the general case, Nonlinearity 22, 2101 (2009).
- [26] B. I. Shklovskii, B. Shapiro, B. R. Sears, P. Lambrianides and H. B. Shore, Statistics of spectra of disordered systems near the metal-insulator transition, Phys. Rev. B 47, 11487 (1993).
- [27] V. E. Kravtsov, I. V. Lerner, B. L. Altshuler and A. G. Aronov, Universal spectral correlations at the mobility edge, Phys. Rev. Lett. 72, 888 (1994).
- [28] T. A. Brody, A statistical measure for the repulsion of energy levels, Lett. Nuovo Cimento 7, 482 (1973).
- [29] F. M. Izrailev, Simple models of quantum chaos: Spectrum and eigenfunctions, Phys. Rep. 196, 299 (1990).
- [30] E. Vergini and M. Saraceno, Calculation by scaling of highly excited states of billiards, Phys. Rev. E 52, 2204 (1995).
- [31] A. Barnett, Dissipation in deforming chaotic billiards, Ph.D. thesis, Harvard University, 2001.
- [32] H. P. Baltes and E. R. Hilf, Spectra of Finite Systems (BIWissenschaftsverlag, Mannheim, 1976).
- [33] K. Husimi, Some formal properties of the density matrix, Proc. Phys. Math. Soc. Jpn. 22, 264 (1940).
- [34] Z. Wang, J. Feng and B. Wu, Microscope for quantum dynamics with Planck cell resolution, Phys. Rev. Research 3, 033239 (2021).
- [35] E.Wigner, On the quantum correction for thermodynamic equilibrium, Phys. Rev. 40, 749 (1932).
- [36] A. Bäcker, S. Fürstberger and R. Schubert, Poincaré Husimi representation of eigenstates in quantum billiards, Phys. Rev. E 70, 036204 (2004).
- [37] W. A. Veech, Teichmüller curves in moduli space, Eisenstein series and an application to triangular billiards, Invent. Math 97, 553 (1989).
- [38] E. B. Bogomolny, B. Dietz, T. Friedrich, M. Miski-Oglu, A. Richter, F. Schäfer and C. Schmit, First experimental observation of superscars in a pseudointegrable barrier billiard, Phys. Rev. Lett. 97, 254102 (2006).
- [39] S. Åberg, T. Guhr, M. Miski-Oglu, and A. Richter, Superscars in Billiards: A Model for Doorway States in Quantum Spectra, Phys. Rev. Lett. 100, 204101 (2008).
- [40] M. V. Berry and M. Robnik, Semiclassical level spacings when regular and chaotic orbits coexist, J. Phys. A 17, 2413 (1984).
- [41] J. Jahnel, When is the (co) sine of a rational angle equal to a rational number?, arXiv:1006.2938[quant-ph], (2010).
- [42] Č. Lozej, G. Casati, and T. Prosen, Quantum chaos in triangular billiards, Phys. Rev. Res. 4, 013138 (2022).