License: CC BY 4.0
arXiv:2401.06474v1 [nlin.CD] 12 Jan 2024

Quantum Dynamical Tunneling Breaks Classical Conserved Quantities

Lingchii Kong(孔令琦) International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China    Zongping Gong(龚宗平) Department of Applied Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan    Biao Wu(吴飙) [email protected] International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China Wilczek Quantum Center, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
(January 12, 2024)
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 105superscript10510^{5}10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT 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 1981198119811981 [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 \mathcal{R}caligraphic_R featuring a multi-handled sphere with genus larger than 1111, in contrast to the genus-1111 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

T(cosθ)cos(Nθ)𝑇𝜃𝑁𝜃T(\cos\theta)\triangleq\cos\left(N\theta\right)italic_T ( roman_cos italic_θ ) ≜ roman_cos ( italic_N italic_θ ) (1)

where θ𝜃\thetaitalic_θ is the angle between the momentum and the horizontal axis and N𝑁Nitalic_N is the least common multiple of nisubscript𝑛𝑖n_{i}italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT for all vertices with rational angles miπ/nisubscript𝑚𝑖𝜋subscript𝑛𝑖m_{i}\pi/n_{i}italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_π / italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT [15]. The conservation arises from the fact that the directions of each classical trajectory follow the orbits of a dihedral group DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT.

Refer to caption
Figure 1: Pseudointegrable systems. (a) Rational polygon and invariant surface \mathcal{R}caligraphic_R with constant energy and T𝑇Titalic_T. Each angle at the vertex i𝑖iitalic_i of polygon has a rational degree miπ/nisubscript𝑚𝑖𝜋subscript𝑛𝑖m_{i}\pi/n_{i}italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_π / italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT with misubscript𝑚𝑖m_{i}italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT and nisubscript𝑛𝑖n_{i}italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT being coprime integers. At least one vertex has mi>1subscript𝑚𝑖1m_{i}>1italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > 1. (b) Rational right triangle. m1subscript𝑚1m_{1}italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and n1subscript𝑛1n_{1}italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT are coprime integers. The directions of momentum can be described by two alternative variables: the angle θ𝜃\thetaitalic_θ of inclination with respect to the horizontal axis, or the included angle θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG with respect to the outer normal vector of sides.

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 \mathcal{R}caligraphic_R 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 T𝑇Titalic_T in their quantum counterparts. The entire proof unfolds in two steps. First, assuming T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG as the quantization of T𝑇Titalic_T is conserved, we can expand the eigenstates by 2N2𝑁2N2 italic_N plane waves with the directions following the orbits of group DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT. 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 T𝑇Titalic_T, we calculate the uncertainty of operator T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG for each eigenstate, defined by

σ(ψn)ψn|T^2|ψnψn|T^|ψn2,𝜎subscript𝜓𝑛quantum-operator-productsubscript𝜓𝑛superscript^𝑇2subscript𝜓𝑛superscriptquantum-operator-productsubscript𝜓𝑛^𝑇subscript𝜓𝑛2\sigma\left(\psi_{n}\right)\triangleq\langle\psi_{n}|\hat{T}^{2}|\psi_{n}% \rangle-\langle\psi_{n}|\hat{T}|\psi_{n}\rangle^{2},italic_σ ( italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) ≜ ⟨ italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT | over^ start_ARG italic_T end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ⟩ - ⟨ italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT | over^ start_ARG italic_T end_ARG | italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ⟩ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (2)

where ψnsubscript𝜓𝑛\psi_{n}italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT denotes the n𝑛nitalic_n-th eigenstate. We observe that the decreasing trend of σ(ψn)𝜎subscript𝜓𝑛\sigma\left(\psi_{n}\right)italic_σ ( italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) with rising energy cannot be described by elementary functions. Instead, we examine the probability density functions (PDFs) of σ𝜎\sigmaitalic_σ within an energy shell, revealing nontrivial distributions that can be accurately fitted by the Brody distribution [28, 29]. The eigenstates with sufficiently small σ𝜎\sigmaitalic_σ are collected, and their unfolded level spacing statistics conform to the Poisson distribution. While, the eigenstates with large σ𝜎\sigmaitalic_σ showcase a superposition with classical regular orbits with different T𝑇Titalic_T, 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 σ𝜎\sigmaitalic_σ, 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 DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT, which fold 0θ2π0𝜃2𝜋0\leq\theta\leq 2\pi0 ≤ italic_θ ≤ 2 italic_π into an interval [0,π/N]0𝜋𝑁[0,\pi/N][ 0 , italic_π / italic_N ]. This corresponds to a quotient map |θmod(2π/N)|modulo𝜃2𝜋𝑁|\theta\mod(2\pi/N)|| italic_θ roman_mod ( 2 italic_π / italic_N ) |, functioning as a conserved quantity. For simplicity, we use the form (1). Consequently, the trajectory in phase space can be projected onto \mathcal{R}caligraphic_R. It is simple to prove T(cosθ)𝑇𝜃T(\cos\theta)italic_T ( roman_cos italic_θ ) is conserved along a classical trajectory, i.e., for a certain reflection at one side, T(cosθ)=cos(Nθ)cos[N(2jπ/N±θ)]=cos(Nθ)=T(cosθ)𝑇𝜃𝑁𝜃𝑁plus-or-minus2𝑗𝜋𝑁𝜃𝑁𝜃𝑇𝜃T(\cos\theta)=\cos\left(N\theta\right)\to\cos\left[N(2j\pi/N\pm\theta)\right]=% \cos\left(N\theta\right)=T(\cos\theta)italic_T ( roman_cos italic_θ ) = roman_cos ( italic_N italic_θ ) → roman_cos [ italic_N ( 2 italic_j italic_π / italic_N ± italic_θ ) ] = roman_cos ( italic_N italic_θ ) = italic_T ( roman_cos italic_θ ) for j=0,1,,N1𝑗01𝑁1j=0,1,\dots,N-1italic_j = 0 , 1 , … , italic_N - 1. Additionally, T(cosθ)𝑇𝜃T(\cos\theta)italic_T ( roman_cos italic_θ ) can be expanded as a polynomial in cosθ𝜃\cos\thetaroman_cos italic_θ, which can be quantized as the operator P^x/Psubscript^𝑃𝑥𝑃\hat{P}_{x}/Pover^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT / italic_P, where P𝑃Pitalic_P is the magnitude of the momentum and P^xsubscript^𝑃𝑥\hat{P}_{x}over^ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT is the momentum operator along x𝑥xitalic_x direction.

III Uncertainty

We calculate the quantum right triangle billiards with one interior angle being π/8𝜋8\pi/8italic_π / 8, π/5𝜋5\pi/5italic_π / 5, π/7𝜋7\pi/7italic_π / 7, 2π/72𝜋72\pi/72 italic_π / 7, π/9𝜋9\pi/9italic_π / 9, 2π/92𝜋92\pi/92 italic_π / 9 (the angles will be used to represent the corresponding billiards). The genus of \mathcal{R}caligraphic_R is given by the formula 1+Ni(mi1)/(2ni)1𝑁subscript𝑖subscript𝑚𝑖12subscript𝑛𝑖1+N\sum_{i}(m_{i}-1)/(2n_{i})1 + italic_N ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - 1 ) / ( 2 italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ). Consequently, each pair of billiards corresponds to \mathcal{R}caligraphic_R with genus-2, 3, 4 . It is noted that the energy level statistics of π/8𝜋8\pi/8italic_π / 8 and π/5𝜋5\pi/5italic_π / 5 has been found to be semi-Poisson distribution [20]. For each billiard, 105superscript10510^{5}10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT eigenstates are calculated using a hybrid of scaling method and decomposition method [30, 31]. We use the Weyl’s level NWeylsubscript𝑁WeylN_{\operatorname{Weyl}}italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT to indicate the energy of eigenstates[32], which helps compare different systems with a unified energy scale. The direct results of σ𝜎\sigmaitalic_σ are depicted in Fig. 2(a). It demonstrates that most points cluster around the small values of σ𝜎\sigmaitalic_σ. However, there are still eigenstates with large σ𝜎\sigmaitalic_σ values at higher energy levels. Since σ𝜎\sigmaitalic_σ does not exhibit a clear decreasing trend with NWeylsubscript𝑁WeylN_{\operatorname{Weyl}}italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT, we analyze σ𝜎\sigmaitalic_σ statistically. If T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG were conserved, we would expect the PDF of uncertainties, denoted as P(σ)𝑃𝜎P(\sigma)italic_P ( italic_σ ), to be a delta distribution at σ=0𝜎0\sigma=0italic_σ = 0. However, the statistical results show that P(σ)𝑃𝜎P(\sigma)italic_P ( italic_σ ) follow various nontrivial distributions. This is exemplified by π/8𝜋8\pi/8italic_π / 8 and π/7𝜋7\pi/7italic_π / 7 billiards in Fig. 2(b1) and Fig. 2(b2). Both P(σ)𝑃𝜎P(\sigma)italic_P ( italic_σ ) 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 T𝑇Titalic_T will recover.

Refer to caption
Figure 2: (a) Scatter plots of σ𝜎\sigmaitalic_σ in relation to NWeylsubscript𝑁WeylN_{\operatorname{Weyl}}italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT for π/8𝜋8\pi/8italic_π / 8 billiard, where 104superscript10410^{4}10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT points are sampled uniformly. (b1) and (b2) PDFs of σ𝜎\sigmaitalic_σ. The 105superscript10510^{5}10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT eigenstates are categorized into three groups based on their energy levels, with the colors representing NWeylsubscript𝑁WeylN_{\operatorname{Weyl}}italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT within the ranges of 02×10402superscript1040-2\times 10^{4}0 - 2 × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT, 46×10446superscript1044-6\times 10^{4}4 - 6 × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT and 810×104810superscript1048-10\times 10^{4}8 - 10 × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT in blue, green, and red, respectively. The insets provide a closer look at σ𝜎\sigmaitalic_σ smaller than 0.10.10.10.1. Notably, (b1) and (b2) respectively correspond to π/8𝜋8\pi/8italic_π / 8 and π/7𝜋7\pi/7italic_π / 7 billiards.

To gain further insights, we consider the unfolded uncertainty, defined as σ~:=σ/mean(σ)assign~𝜎𝜎mean𝜎\widetilde{\sigma}:=\sigma/\operatorname{mean}(\sigma)over~ start_ARG italic_σ end_ARG := italic_σ / roman_mean ( italic_σ ), whose domain is extended from [0,1]01[0,1][ 0 , 1 ] to [0,]0[0,\infty][ 0 , ∞ ]. The PDFs of σ~~𝜎\widetilde{\sigma}over~ start_ARG italic_σ end_ARG, denoted as P(σ~)𝑃~𝜎P(\widetilde{\sigma})italic_P ( over~ start_ARG italic_σ end_ARG ), are shown in Fig. 3(a). All P(σ~)𝑃~𝜎P(\widetilde{\sigma})italic_P ( over~ start_ARG italic_σ end_ARG ) presented here can be well fitted by Brody distributions PBrody(σ~)=a(q+1)σ~qexp(aσ~q+1)subscript𝑃Brody~𝜎𝑎𝑞1superscript~𝜎𝑞𝑎superscript~𝜎𝑞1P_{\operatorname{Brody}}(\widetilde{\sigma})=a(q+1)\widetilde{\sigma}^{q}\exp% \left(-a\widetilde{\sigma}^{q+1}\right)italic_P start_POSTSUBSCRIPT roman_Brody end_POSTSUBSCRIPT ( over~ start_ARG italic_σ end_ARG ) = italic_a ( italic_q + 1 ) over~ start_ARG italic_σ end_ARG start_POSTSUPERSCRIPT italic_q end_POSTSUPERSCRIPT roman_exp ( - italic_a over~ start_ARG italic_σ end_ARG start_POSTSUPERSCRIPT italic_q + 1 end_POSTSUPERSCRIPT ), a=Γ[(q+2)/(q+1)]q+1𝑎Γsuperscriptdelimited-[]𝑞2𝑞1𝑞1a=\Gamma\left[(q+2)/(q+1)\right]^{q+1}italic_a = roman_Γ [ ( italic_q + 2 ) / ( italic_q + 1 ) ] start_POSTSUPERSCRIPT italic_q + 1 end_POSTSUPERSCRIPT with a single parameter q[0,1]𝑞01q\in[0,1]italic_q ∈ [ 0 , 1 ]. In particular, for π/8𝜋8\pi/8italic_π / 8 and π/5𝜋5\pi/5italic_π / 5 billiards, whose level statistics are semi-Poissonian, P(σ~)𝑃~𝜎P(\widetilde{\sigma})italic_P ( over~ start_ARG italic_σ end_ARG ) closely resemble the q=0𝑞0q=0italic_q = 0 Poisson distribution exp(σ~)~𝜎\exp\left(-\widetilde{\sigma}\right)roman_exp ( - over~ start_ARG italic_σ end_ARG ). The differences in cumulative distribution functions (CDFs) from the best-fitted model, denoted by UB(σ~):=0σ~𝑑σ~1[P(σ~1)PBrody*(σ~1)]assignsubscript𝑈𝐵~𝜎superscriptsubscript0~𝜎differential-dsubscript~𝜎1delimited-[]𝑃subscript~𝜎1superscriptsubscript𝑃Brodysubscript~𝜎1U_{B}(\widetilde{\sigma}):=\int_{0}^{\widetilde{\sigma}}d\widetilde{\sigma}_{1% }\left[P(\widetilde{\sigma}_{1})-P_{\operatorname{Brody}}^{*}(\widetilde{% \sigma}_{1})\right]italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( over~ start_ARG italic_σ end_ARG ) := ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over~ start_ARG italic_σ end_ARG end_POSTSUPERSCRIPT italic_d over~ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_P ( over~ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) - italic_P start_POSTSUBSCRIPT roman_Brody end_POSTSUBSCRIPT start_POSTSUPERSCRIPT * end_POSTSUPERSCRIPT ( over~ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ] where the superscript *** represents “optimal fitting”, are shown in the inset, which give the errors less than 2%percent22\%2 %. This result shows that P(σ~)𝑃~𝜎P(\widetilde{\sigma})italic_P ( over~ start_ARG italic_σ end_ARG ) for genus-3, 4 billiards have a polynomial repulsion at small σ~~𝜎\widetilde{\sigma}over~ start_ARG italic_σ end_ARG and a super-exponential tail at large σ~~𝜎\widetilde{\sigma}over~ start_ARG italic_σ end_ARG. Roughly, P(σ~)𝑃~𝜎P(\widetilde{\sigma})italic_P ( over~ start_ARG italic_σ end_ARG ) exhibits an elevated q𝑞qitalic_q 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 m1/n1subscript𝑚1subscript𝑛1m_{1}/n_{1}italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is a good Diophantine approximation of some irrational number. Then, N𝑁Nitalic_N is typically large, resulting in a high genus. Since irrational triangle billiards lack the conserved quantity T𝑇Titalic_T, a \mathcal{R}caligraphic_R with higher genus naturally implies a more chaotic system.

Refer to caption
Figure 3: (a) PDFs of σ~~𝜎\widetilde{\sigma}over~ start_ARG italic_σ end_ARG for 105superscript10510^{5}10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT eigenstates. The shaded band represents Poisson distribution, i.e., exp(σ~)~𝜎\exp\left(-\widetilde{\sigma}\right)roman_exp ( - over~ start_ARG italic_σ end_ARG ). These distributions are fitted by the Brody distribution. The inset displays the differences in CDFs from their optimal fitting distributions, denoted as UB(σ~)subscript𝑈𝐵~𝜎U_{B}(\widetilde{\sigma})italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( over~ start_ARG italic_σ end_ARG ). The optimal fitting parameters q𝑞qitalic_q are approximately 1.01.01.01.0, 1.01.01.01.0, 1.41.41.41.4, 1.51.51.51.5, 1.51.51.51.5, 1.91.91.91.9 for π/8𝜋8\pi/8italic_π / 8, π/5𝜋5\pi/5italic_π / 5, π/7𝜋7\pi/7italic_π / 7, 2π/72𝜋72\pi/72 italic_π / 7, π/9𝜋9\pi/9italic_π / 9, 2π/92𝜋92\pi/92 italic_π / 9 billiards, respectively. (b) Differences in CDFs of unfolded level spacing with different uncertainty thresholds from the Poisson distribution, denoted as UP(s~)subscript𝑈𝑃~𝑠U_{P}(\widetilde{s})italic_U start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ). σsuperscript𝜎\sigma^{\prime}italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT is the selected threshold value. For π/8𝜋8\pi/8italic_π / 8 billiard, σcsubscript𝜎𝑐\sigma_{c}italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is chosen to be 0.01890.01890.01890.0189 with the number of pseudoregular states amounting to 11481114811148111481. Notably, in the absence of uncertainty threshold (σ=1superscript𝜎1\sigma^{\prime}=1italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1), the level spacing statistics for the whole spectrum follows PsemiPoisson(s~)subscript𝑃semiPoisson~𝑠P_{\operatorname{semi-Poisson}}(\widetilde{s})italic_P start_POSTSUBSCRIPT roman_semi - roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ), as indicated by the red dashed line. Additionally, the level spacing distribution for the uncertainties less than 0.140.140.140.14 is shown at the intermediate position between Poisson and semi-Poisson distribution.

To confirm our expectation that the level statistics of eigenstates with small enough σ𝜎\sigmaitalic_σ follows a Poisson distribution, we scan the values of σ𝜎\sigmaitalic_σ to determine a small enough threshold σcsubscript𝜎𝑐\sigma_{c}italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, such that the unfolded level spacing distribution of eigenstates with σ<σc𝜎subscript𝜎𝑐\sigma<\sigma_{c}italic_σ < italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT aligns with the Poisson distribution with minimal fitting error. This can be formulated as P(s~|{ψn,σ(ψn)<σc})PPoisson(s~)similar-to-or-equals𝑃conditional~𝑠subscript𝜓𝑛𝜎subscript𝜓𝑛subscript𝜎𝑐subscript𝑃Poisson~𝑠P\left(\widetilde{s}\;|\left\{\psi_{n},\sigma(\psi_{n})<\sigma_{c}\right\}% \right)\simeq P_{\operatorname{Poisson}}\left(\widetilde{s}\right)italic_P ( over~ start_ARG italic_s end_ARG | { italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_σ ( italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) < italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT } ) ≃ italic_P start_POSTSUBSCRIPT roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ) where PPoisson(s~)=exp(s~)subscript𝑃Poisson~𝑠~𝑠P_{\operatorname{Poisson}}(\widetilde{s})=\exp\left(-\widetilde{s}\right)italic_P start_POSTSUBSCRIPT roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ) = roman_exp ( - over~ start_ARG italic_s end_ARG ) and s~~𝑠\widetilde{s}over~ start_ARG italic_s end_ARG is the unfolded level spacing. We call the eigenstates with σ<σc𝜎subscript𝜎𝑐\sigma<\sigma_{c}italic_σ < italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT as “pseudoregular” states. The results for π/8𝜋8\pi/8italic_π / 8 billiard are depicted in Fig. 3(b), where the differences in CDFs of unfolded level spacing with different thresholds σsuperscript𝜎\sigma^{\prime}italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT from the Poisson distribution are calculated, defined as UP(s~):=0s~𝑑s~1[P(s~1|{ψn,σ(ψn)<σ})PPoisson(s~1)]assignsubscript𝑈𝑃~𝑠superscriptsubscript0~𝑠differential-dsubscript~𝑠1delimited-[]𝑃conditionalsubscript~𝑠1subscript𝜓𝑛𝜎subscript𝜓𝑛superscript𝜎subscript𝑃Poissonsubscript~𝑠1U_{P}(\widetilde{s}):=\int_{0}^{\widetilde{s}}d\widetilde{s}_{1}\left[P\left(% \widetilde{s}_{1}\;|\left\{\psi_{n},\sigma(\psi_{n})<\sigma^{\prime}\right\}% \right)-P_{\operatorname{Poisson}}(\widetilde{s}_{1})\right]italic_U start_POSTSUBSCRIPT italic_P end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ) := ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over~ start_ARG italic_s end_ARG end_POSTSUPERSCRIPT italic_d over~ start_ARG italic_s end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_P ( over~ start_ARG italic_s end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | { italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_σ ( italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) < italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT } ) - italic_P start_POSTSUBSCRIPT roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ]. It shows that the level statistics of pseudoregular states approximates to PPoisson(s~)subscript𝑃Poisson~𝑠P_{\operatorname{Poisson}}(\widetilde{s})italic_P start_POSTSUBSCRIPT roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ), implying an absence of level repulsions in the spectrum if T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG is conserved. It is established that when all the eigenstates are involved (σ=1superscript𝜎1\sigma^{\prime}=1italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1), the level spacing statistics for the whole spectrum follows a semi-Poisson distribution, i.e., PsemiPoisson(s~)=4s~exp(2s~)subscript𝑃semiPoisson~𝑠4~𝑠2~𝑠P_{\operatorname{semi-Poisson}}(\widetilde{s})=4\widetilde{s}\exp\left(-2% \widetilde{s}\right)italic_P start_POSTSUBSCRIPT roman_semi - roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ) = 4 over~ start_ARG italic_s end_ARG roman_exp ( - 2 over~ start_ARG italic_s end_ARG ). By increasing the threshold to σ=0.14>σcsuperscript𝜎0.14subscript𝜎𝑐\sigma^{\prime}=0.14>\sigma_{c}italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 0.14 > italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, the level statistics occupies an intermediate position between the Poisson and semi-Poisson distribution. This demonstrates that by controlling the uncertainty threshold σsuperscript𝜎\sigma^{\prime}italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, the level statistics can be changed between these two distributions. We will see that the eigenstates with σ>σc𝜎subscript𝜎𝑐\sigma>\sigma_{c}italic_σ > italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT are characterized by the superposition of classical trajectories with distinct values of T𝑇Titalic_T, which is caused by dynamical tunneling.

Refer to caption
Figure 4: Profiles of different types of eigenstates for π/8𝜋8\pi/8italic_π / 8 billiard, illustrated by un(θ¯)subscript𝑢𝑛¯𝜃u_{n}(\bar{\theta})italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) with respect to |T(cosθ¯)|𝑇¯𝜃\left|T(\cos\bar{\theta})\right|| italic_T ( roman_cos over¯ start_ARG italic_θ end_ARG ) |. (a) Superscar state with NWeyl10017.3subscript𝑁Weyl10017.3N_{\operatorname{Weyl}}\approx 10017.3italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 10017.3 and σ0.006𝜎0.006\sigma\approx 0.006italic_σ ≈ 0.006. (b) Spacial uniformly distributed orbit with NWeyl40090.0subscript𝑁Weyl40090.0N_{\operatorname{Weyl}}\approx 40090.0italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 40090.0 and σ0.002𝜎0.002\sigma\approx 0.002italic_σ ≈ 0.002. (c1) Superposition of a single superscar and a single UDO with NWeyl60144.6subscript𝑁Weyl60144.6N_{\operatorname{Weyl}}\approx 60144.6italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 60144.6 and σ0.563𝜎0.563\sigma\approx 0.563italic_σ ≈ 0.563. (c2) is the eigenstate next to (c1) with NWeyl60144.1subscript𝑁Weyl60144.1N_{\operatorname{Weyl}}\approx 60144.1italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 60144.1 and σ0.489𝜎0.489\sigma\approx 0.489italic_σ ≈ 0.489. (d1) Superposition of two UDOs with NWeyl61476.4subscript𝑁Weyl61476.4N_{\operatorname{Weyl}}\approx 61476.4italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 61476.4 and σ0.461𝜎0.461\sigma\approx 0.461italic_σ ≈ 0.461. (d2) is the eigenstate next to (d1) with NWeyl61476.6subscript𝑁Weyl61476.6N_{\operatorname{Weyl}}\approx 61476.6italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 61476.6 and σ0.509𝜎0.509\sigma\approx 0.509italic_σ ≈ 0.509.

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 T𝑇Titalic_T. As a result, each eigenstate that may be a superposition of different classical orbits can be converted into a distribution of different T𝑇Titalic_T values . This distribution is expressed by the boundary function un(θ¯)subscript𝑢𝑛¯𝜃u_{n}(\bar{\theta})italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) of ψnsubscript𝜓𝑛\psi_{n}italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT, where θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG 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 θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG, the conserved quantity T𝑇Titalic_T can be rewritten as T=|T(cosθ¯)|=|cos(Nθ¯)|𝑇𝑇¯𝜃𝑁¯𝜃T=|T(\cos\bar{\theta})|=|\cos\left(N\bar{\theta}\right)|italic_T = | italic_T ( roman_cos over¯ start_ARG italic_θ end_ARG ) | = | roman_cos ( italic_N over¯ start_ARG italic_θ end_ARG ) |, where the absolute value arises at the hypotenuse θ¯=π/2+mπ/nθ¯𝜃𝜋2𝑚𝜋𝑛𝜃\bar{\theta}=\pi/2+m\pi/n-\thetaover¯ start_ARG italic_θ end_ARG = italic_π / 2 + italic_m italic_π / italic_n - italic_θ. In this way, θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG becomes the connection between un(θ¯)subscript𝑢𝑛¯𝜃u_{n}\left(\bar{\theta}\right)italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) and T𝑇Titalic_T. The results of π/8𝜋8\pi/8italic_π / 8 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 T=1𝑇1T=1italic_T = 1, while UDOs have various values of T𝑇Titalic_T. Quantum mechanically, there exists “superscar” states, named by Bogomolny [22, 38, 39], exhibiting a superposition of spatially parallel POs. T𝑇Titalic_T cannot be used to distinguish different superscar states as they all equal to 1111. In both Fig. 4.(a) and Fig. 4.(b), |un(θ¯)|2superscriptsubscript𝑢𝑛¯𝜃2\left|u_{n}(\bar{\theta})\right|^{2}| italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT exhibit a single peak, illustrating the features of superscar and UDO, respectively. Both of them are pseudoregular states, characterized by uncertainties σ𝜎\sigmaitalic_σ smaller than σcsubscript𝜎𝑐\sigma_{c}italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. However, in Fig. 4.(c1) and Fig. 4.(d1), orbits with different values of T𝑇Titalic_T are superposed (one superscar plus one UDO in (c1), two UDOs in (d1)), causing their relatively large σ𝜎\sigmaitalic_σ. Notably, eigenstates with small wavenumbers may distribute with a considerable width on the values of T𝑇Titalic_T, leading to an elevated σ𝜎\sigmaitalic_σ. 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 1111. This is another indication of dynamical tunneling. In Fig. 4.(c2), there emerges sub-peaks at small T𝑇Titalic_T comparing to Fig. 4.(c1), which may illustrate more resonant orbits. Therefore, the results indicate that dynamical tunneling destroys the conservation of T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG by superposing trajectories with different directions.

V Avoided Crossing

The dynamical tunneling can be modeled by a random matrix Hamiltonian

h=(ϵγγϵ),matrixitalic-ϵ𝛾𝛾italic-ϵh=\begin{pmatrix}\epsilon&\gamma\\ \gamma&-\epsilon\end{pmatrix}\,,italic_h = ( start_ARG start_ROW start_CELL italic_ϵ end_CELL start_CELL italic_γ end_CELL end_ROW start_ROW start_CELL italic_γ end_CELL start_CELL - italic_ϵ end_CELL end_ROW end_ARG ) , (3)

where the diagonals ±ϵplus-or-minusitalic-ϵ\pm\epsilon± italic_ϵ are energy levels of different classical orbits, and the off-diagonal γ𝛾\gammaitalic_γ represent tunneling rate between classical orbits with different values of T𝑇Titalic_T. γ𝛾\gammaitalic_γ is real due to the time-reversal symmetry. As different classical orbits are not correlated, it is reasonable to assume that ϵitalic-ϵ\epsilonitalic_ϵ has a Poisson distribution, i.e., P(ϵ)=exp(2|ϵ|)𝑃italic-ϵ2italic-ϵP(\epsilon)=\exp\left(-2|\epsilon|\right)italic_P ( italic_ϵ ) = roman_exp ( - 2 | italic_ϵ | ).

The numerical results in Fig. 2(b1) and Fig. 3(a) suggest that the distribution g(γ)𝑔𝛾g(\gamma)italic_g ( italic_γ ) of γ𝛾\gammaitalic_γ have the same shape as P(σ)𝑃𝜎P(\sigma)italic_P ( italic_σ ), that is,

g(γ)=g(γ)12exp(|γ|).𝑔𝛾𝑔𝛾12𝛾g(\gamma)=g(-\gamma)\approx\frac{1}{2}\exp\left(-{|\gamma|}\right).italic_g ( italic_γ ) = italic_g ( - italic_γ ) ≈ divide start_ARG 1 end_ARG start_ARG 2 end_ARG roman_exp ( - | italic_γ | ) . (4)

Once the two distributions P(ϵ)𝑃italic-ϵP(\epsilon)italic_P ( italic_ϵ ) and g(γ)𝑔𝛾g(\gamma)italic_g ( italic_γ ) are known, we can compute the level spacing distribution of (3), which is compared with PsemiPoisson(s~)subscript𝑃semiPoisson~𝑠P_{\operatorname{semi-Poisson}}(\widetilde{s})italic_P start_POSTSUBSCRIPT roman_semi - roman_Poisson end_POSTSUBSCRIPT ( over~ start_ARG italic_s end_ARG ) 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.

Refer to caption
Figure 5: Unfolded level spacing distribution of model (3) using the tunneling rate (4). The level spacings are collected from 106superscript10610^{6}10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT ensembles. The red line represents the semi-Poisson distribution.

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 T𝑇Titalic_T governing the entire systems by superposing classical orbits with different values of T𝑇Titalic_T (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 1superscriptPlanck-constant-over-2-pi1\hbar^{-1}roman_ℏ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, which can be expressed as lnγξ/proportional-to𝛾𝜉Planck-constant-over-2-pi\ln\gamma\propto-\xi/\hbarroman_ln italic_γ ∝ - italic_ξ / roman_ℏ where ξ𝜉\xiitalic_ξ is a parameter related to systems and resonant orbits (for a double-well system , ξ𝜉\xiitalic_ξ 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 1superscriptPlanck-constant-over-2-pi1\hbar^{-1}roman_ℏ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT [9, 10, 11, 12, 13]. We formulate this qualitatively as ξ=ξ(1)=ξ¯+iλiδ(1i1)𝜉𝜉superscriptPlanck-constant-over-2-pi1¯𝜉subscript𝑖subscript𝜆𝑖𝛿superscriptPlanck-constant-over-2-pi1superscriptsubscriptPlanck-constant-over-2-pi𝑖1\xi=\xi(\hbar^{-1})=\bar{\xi}+\sum_{i}\lambda_{i}\delta(\hbar^{-1}-\hbar_{i}^{% -1})italic_ξ = italic_ξ ( roman_ℏ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) = over¯ start_ARG italic_ξ end_ARG + ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_δ ( roman_ℏ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT - roman_ℏ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) where ξ¯¯𝜉\bar{\xi}over¯ start_ARG italic_ξ end_ARG is the average value over a shell of 1superscriptPlanck-constant-over-2-pi1\hbar^{-1}roman_ℏ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, i1superscriptsubscriptPlanck-constant-over-2-pi𝑖1\hbar_{i}^{-1}roman_ℏ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT is the enhanced (or suppressed) point and λisubscript𝜆𝑖\lambda_{i}italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the corresponding strength. In the limit of strong tunneling, the distribution of i1superscriptsubscriptPlanck-constant-over-2-pi𝑖1\hbar_{i}^{-1}roman_ℏ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT becomes dense on the shell, allowing us to treat ξ𝜉\xiitalic_ξ as a random variable centered around ξ¯¯𝜉\bar{\xi}over¯ start_ARG italic_ξ end_ARG. 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 θ𝜃\thetaitalic_θ can undergo three kinds of reflections by sides: θθ𝜃𝜃\theta\to-\thetaitalic_θ → - italic_θ (reflection at y=0𝑦0y=0italic_y = 0), θπθ𝜃𝜋𝜃\theta\to\pi-\thetaitalic_θ → italic_π - italic_θ (reflection at x=Lcosα𝑥𝐿𝛼x=L\cos\alphaitalic_x = italic_L roman_cos italic_α) and θ2αθ𝜃2𝛼𝜃\theta\to 2\alpha-\thetaitalic_θ → 2 italic_α - italic_θ (reflection at y=xtanα𝑦𝑥𝛼y=x\tan\alphaitalic_y = italic_x roman_tan italic_α). These can generate a dihedral group DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT (N=n1𝑁subscript𝑛1N=n_{1}italic_N = italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT) with dimension 2N2𝑁2N2 italic_N, yielding a conserved quantity T=cosNθ𝑇𝑁𝜃T=\cos N\thetaitalic_T = roman_cos italic_N italic_θ independent with the energy. T𝑇Titalic_T can be quantized after expanding to the polynomial of cosθ=Px/P𝜃subscript𝑃𝑥𝑃\cos\theta={P}_{x}/Proman_cos italic_θ = italic_P start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT / italic_P where Pxsubscript𝑃𝑥P_{x}italic_P start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT is the momentum along x𝑥xitalic_x-direction and P𝑃Pitalic_P is the norm of total momentum.

Refer to caption
Figure 6: Rational right triangle billiard. This billiard includes at least one vertex with angle α𝛼\alphaitalic_α being m1π/n1subscript𝑚1𝜋subscript𝑛1m_{1}\pi/n_{1}italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_π / italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT where n1subscript𝑛1n_{1}italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is even, m1subscript𝑚1m_{1}italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (m11subscript𝑚11m_{1}\neq 1italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≠ 1) and n1subscript𝑛1n_{1}italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT are coprime integers. L𝐿Litalic_L is the length of hypotenuse. The directions of billiard can be described by either of two alternative variables: the angle θ𝜃\thetaitalic_θ of inclination with respect to the x𝑥xitalic_x-axis, or the included angle θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG with respect to the outer normal vector of sides.

The conserved quantity T𝑇Titalic_T originates from the existence of quotient projection mapping each classical trajectory on the surface S1/DNsuperscript𝑆1subscript𝐷𝑁S^{1}/D_{N}italic_S start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT / italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT, where S1superscript𝑆1S^{1}italic_S start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT is the unit circle. It means that for a classical trajectory, the possible directions, denoted as ±θjplus-or-minussubscript𝜃𝑗\pm\theta_{j}± italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT, can be related to a reference θ0[0,π/N]subscript𝜃00𝜋𝑁\theta_{0}\in[0,\pi/N]italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ [ 0 , italic_π / italic_N ] so that

θj=θ0+2jπN,j=0,1,,N1.formulae-sequencesubscript𝜃𝑗subscript𝜃02𝑗𝜋𝑁𝑗01𝑁1\theta_{j}=\theta_{0}+\frac{2j\pi}{N},\quad j=0,1,\dots,N-1.italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG 2 italic_j italic_π end_ARG start_ARG italic_N end_ARG , italic_j = 0 , 1 , … , italic_N - 1 . (A5)

If the quantized operator T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG remains conserved in quantum mechanics, the eigenstate should also follow the orbits of group DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT, just as each classical trajectory does. In this manner, the eigenstate, denoted as |ψket𝜓|\psi\rangle| italic_ψ ⟩, can be expanded by plane waves with direction ±θjplus-or-minussubscript𝜃𝑗\pm\theta_{j}± italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT,

|ψ=j=0N1(cj+|kcosθj,ksinθj+cj|kcosθj,ksinθj)ket𝜓superscriptsubscript𝑗0𝑁1superscriptsubscript𝑐𝑗ket𝑘subscript𝜃𝑗𝑘subscript𝜃𝑗superscriptsubscript𝑐𝑗ket𝑘subscript𝜃𝑗𝑘subscript𝜃𝑗|\psi\rangle=\sum_{j=0}^{N-1}\left(c_{j}^{+}\left|k\cos\theta_{j},k\sin\theta_% {j}\right\rangle+c_{j}^{-}\left|k\cos\theta_{j},-k\sin\theta_{j}\right\rangle\right)| italic_ψ ⟩ = ∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT | italic_k roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_k roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ⟩ + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT | italic_k roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , - italic_k roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ⟩ ) (A6)

where cj±superscriptsubscript𝑐𝑗plus-or-minusc_{j}^{\pm}italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are the superposing parameters. Then, the Dirichlet boundary conditions at y=0𝑦0y=0italic_y = 0, y=xtanα𝑦𝑥𝛼y=x\tan\alphaitalic_y = italic_x roman_tan italic_α and x=Lcosα𝑥𝐿𝛼x=L\cos\alphaitalic_x = italic_L roman_cos italic_α can be evaluated sequentially:

j=0N1(cj++cj)eikxcosθj=0,superscriptsubscript𝑗0𝑁1superscriptsubscript𝑐𝑗superscriptsubscript𝑐𝑗superscript𝑒𝑖𝑘𝑥subscript𝜃𝑗0\displaystyle\sum_{j=0}^{N-1}\left(c_{j}^{+}+c_{j}^{-}\right)e^{ikx\cos\theta_% {j}}=0,∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_x roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = 0 , (A7)
j=0N1[cj+eikrcos(θjα)+cjeikrcos(θj+α)]=0,superscriptsubscript𝑗0𝑁1delimited-[]superscriptsubscript𝑐𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼superscriptsubscript𝑐𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼0\displaystyle\sum_{j=0}^{N-1}\left[c_{j}^{+}e^{ikr\cos\left(\theta_{j}-\alpha% \right)}+c_{j}^{-}e^{ikr\cos\left(\theta_{j}+\alpha\right)}\right]=0,∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT [ italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_α ) end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_α ) end_POSTSUPERSCRIPT ] = 0 , (A8)
j=0N1(cj+eikysinθj+cjeikysinθj)eikLcosαcosθj=0,superscriptsubscript𝑗0𝑁1superscriptsubscript𝑐𝑗superscript𝑒𝑖𝑘𝑦subscript𝜃𝑗superscriptsubscript𝑐𝑗superscript𝑒𝑖𝑘𝑦subscript𝜃𝑗superscript𝑒𝑖𝑘𝐿𝛼subscript𝜃𝑗0\displaystyle\sum_{j=0}^{N-1}\left(c_{j}^{+}e^{iky\sin\theta_{j}}+c_{j}^{-}e^{% -iky\sin\theta_{j}}\right)e^{ikL\cos\alpha\cos\theta_{j}}=0,∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_y roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_k italic_y roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_L roman_cos italic_α roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = 0 , (A9)

where r=x2+y2𝑟superscript𝑥2superscript𝑦2absentr=\sqrt{x^{2}+y^{2}}\initalic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∈ [0,L]0𝐿[0,L][ 0 , italic_L ]. Suppose the subscript of cj±superscriptsubscript𝑐𝑗plus-or-minusc_{j}^{\pm}italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT has a cyclic relation cj±:=cjmodN±assignsuperscriptsubscript𝑐𝑗plus-or-minussuperscriptsubscript𝑐modulo𝑗𝑁plus-or-minusc_{j}^{\pm}:=c_{j\bmod N}^{\pm}italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT := italic_c start_POSTSUBSCRIPT italic_j roman_mod italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT for j𝑗j\in\mathbb{Z}italic_j ∈ blackboard_Z, Eq. (A8) can be rewritten as

j=0N1(cj+m++cj)eikrcos(θj+α)=0.superscriptsubscript𝑗0𝑁1superscriptsubscript𝑐𝑗𝑚superscriptsubscript𝑐𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼0\sum_{j=0}^{N-1}\left(c_{j+m}^{+}+c_{j}^{-}\right)e^{ikr\cos\left(\theta_{j}+% \alpha\right)}=0.∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j + italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_α ) end_POSTSUPERSCRIPT = 0 . (A10)

We can prove that Eq. (A7) and Eq. (A10) are equivalent to

cj++cj=cj+m++cj=0,j=0,1,,N1.formulae-sequencesuperscriptsubscript𝑐𝑗superscriptsubscript𝑐𝑗superscriptsubscript𝑐𝑗𝑚superscriptsubscript𝑐𝑗0𝑗01𝑁1c_{j}^{+}+c_{j}^{-}=c_{j+m}^{+}+c_{j}^{-}=0,\quad j=0,1,\ldots,N-1.italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = italic_c start_POSTSUBSCRIPT italic_j + italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 , italic_j = 0 , 1 , … , italic_N - 1 . (A11)

We take the i𝑖iitalic_i-th (i=1,2,,N1)𝑖12𝑁1(i=1,2,\ldots,N-1)( italic_i = 1 , 2 , … , italic_N - 1 ) derivative of the both sides of Eq. (A7) and make x𝑥xitalic_x equal to zero. In this way, we obtain N𝑁Nitalic_N homogeneous linear equations by defining aj:=cj++cjassignsubscript𝑎𝑗superscriptsubscript𝑐𝑗superscriptsubscript𝑐𝑗a_{j}:=c_{j}^{+}+c_{j}^{-}italic_a start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT := italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT,

[111cosθ0cosθ1cosθN1cos2θ0cos2θ1cos2θN1cosN1θ0cosN1θ1cosN1θN1][a0a1aN1]=0.delimited-[]111missing-subexpressionsubscript𝜃0subscript𝜃1subscript𝜃𝑁1missing-subexpressionsuperscript2subscript𝜃0superscript2subscript𝜃1superscript2subscript𝜃𝑁1missing-subexpressionmissing-subexpressionsuperscript𝑁1subscript𝜃0superscript𝑁1subscript𝜃1superscript𝑁1subscript𝜃𝑁1missing-subexpressiondelimited-[]subscript𝑎0subscript𝑎1subscript𝑎𝑁10\left[\begin{array}[]{ccccc}1&1&\ldots&1\\ \cos\theta_{0}&\cos\theta_{1}&\ldots&\cos\theta_{N-1}\\ \cos^{2}\theta_{0}&\cos^{2}\theta_{1}&\ldots&\cos^{2}\theta_{N-1}\\ \vdots&\vdots&\ddots&\vdots\\ \cos^{N-1}\theta_{0}&\cos^{N-1}\theta_{1}&\ldots&\cos^{N-1}\theta_{N-1}\end{% array}\right]\left[\begin{array}[]{c}a_{0}\\ a_{1}\\ \vdots\\ a_{N-1}\end{array}\right]=0.[ start_ARRAY start_ROW start_CELL 1 end_CELL start_CELL 1 end_CELL start_CELL … end_CELL start_CELL 1 end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_cos italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL start_CELL roman_cos italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL … end_CELL start_CELL roman_cos italic_θ start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL start_CELL roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL … end_CELL start_CELL roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL ⋮ end_CELL start_CELL ⋮ end_CELL start_CELL ⋱ end_CELL start_CELL ⋮ end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_cos start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL start_CELL roman_cos start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL … end_CELL start_CELL roman_cos start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW end_ARRAY ] [ start_ARRAY start_ROW start_CELL italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL ⋮ end_CELL end_ROW start_ROW start_CELL italic_a start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL end_ROW end_ARRAY ] = 0 . (A12)

Its determinant is known as Vandermonde determinant, yielding

det[cosi1θj]N×N=0i<jN1(cosθjcosθi).\operatorname{det}\left[\cos^{i-1}\theta_{j}\right]_{N\times N}=\prod_{0\leq i% <j\leq N-1}\left(\cos\theta_{j}-\cos\theta_{i}\right).roman_det [ roman_cos start_POSTSUPERSCRIPT italic_i - 1 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT italic_N × italic_N end_POSTSUBSCRIPT = ∏ start_POSTSUBSCRIPT 0 ≤ italic_i < italic_j ≤ italic_N - 1 end_POSTSUBSCRIPT ( roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - roman_cos italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) . (A13)

If let cosθj=cosθisubscript𝜃𝑗subscript𝜃𝑖\cos\theta_{j}=\cos\theta_{i}roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = roman_cos italic_θ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, we obtain θ0=(Nij)π/Nsubscript𝜃0𝑁𝑖𝑗𝜋𝑁\theta_{0}=({N-i-j})\pi/{N}italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = ( italic_N - italic_i - italic_j ) italic_π / italic_N which cannot be confined in the domain (0,π/N)0𝜋𝑁\left(0,{\pi}/{N}\right)( 0 , italic_π / italic_N ). As a result, this determinant is nonzero for θ0(0,π/N)subscript𝜃00𝜋𝑁\theta_{0}\in\left(0,{\pi}/{N}\right)italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_π / italic_N ). In this way, the solution to Eq. (A12) is a null vector, i.e.,

aj=cj++cj=0.subscript𝑎𝑗superscriptsubscript𝑐𝑗superscriptsubscript𝑐𝑗0a_{j}=c_{j}^{+}+c_{j}^{-}=0.italic_a start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 . (A14)

When θ0=0,π/Nsubscript𝜃00𝜋𝑁\theta_{0}=0,\pi/Nitalic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0 , italic_π / italic_N, the billiard has N𝑁Nitalic_N directions corresponding to degenerate cases with θj=θNjsubscript𝜃𝑗subscript𝜃𝑁𝑗\theta_{j}=-\theta_{N-j}italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = - italic_θ start_POSTSUBSCRIPT italic_N - italic_j end_POSTSUBSCRIPT. We redefine the superposing parameter as cjsubscript𝑐𝑗c_{j}italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT, so that |ψ=j=0N1cj|kcosθj,ksinθjket𝜓superscriptsubscript𝑗0𝑁1subscript𝑐𝑗ket𝑘subscript𝜃𝑗𝑘subscript𝜃𝑗|\psi\rangle=\sum_{j=0}^{N-1}c_{j}\left|k\cos\theta_{j},k\sin\theta_{j}\right\rangle| italic_ψ ⟩ = ∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | italic_k roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_k roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ⟩. For θ0=0subscript𝜃00\theta_{0}=0italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0, the degenerate version of boundary conditions Eq. (A7) and Eq. (A8) become

c0eikx+cN/2eikx+j=1N/21(cj+cNj)eikxcosθj=0,subscript𝑐0superscript𝑒𝑖𝑘𝑥subscript𝑐𝑁2superscript𝑒𝑖𝑘𝑥superscriptsubscript𝑗1𝑁21subscript𝑐𝑗subscript𝑐𝑁𝑗superscript𝑒𝑖𝑘𝑥subscript𝜃𝑗0\displaystyle c_{0}e^{ikx}+c_{N/2}e^{-ikx}+\sum_{j=1}^{N/2-1}\left(c_{j}+c_{N-% j}\right)e^{ikx\cos\theta_{j}}=0,italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_x end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_N / 2 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_k italic_x end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N / 2 - 1 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_N - italic_j end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_x roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = 0 , (A15)
j=0(m1)/2(cj+cmj)eikrcos(θjα)+limit-fromsuperscriptsubscript𝑗0𝑚12subscript𝑐𝑗subscript𝑐𝑚𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼\displaystyle\sum_{j=0}^{(m-1)/2}(c_{j}+c_{m-j})e^{ikr\cos\left(\theta_{j}-% \alpha\right)}+∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m - 1 ) / 2 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_m - italic_j end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_α ) end_POSTSUPERSCRIPT +
j=m+1(m+N1)/2(cj+cm+Nj)eikrcos(θjα)=0.superscriptsubscript𝑗𝑚1𝑚𝑁12subscript𝑐𝑗subscript𝑐𝑚𝑁𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼0\displaystyle\sum_{j=m+1}^{(m+N-1)/2}(c_{j}+c_{m+N-j})e^{ikr\cos\left(\theta_{% j}-\alpha\right)}=0.∑ start_POSTSUBSCRIPT italic_j = italic_m + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m + italic_N - 1 ) / 2 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_m + italic_N - italic_j end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_α ) end_POSTSUPERSCRIPT = 0 . (A16)

These two relations are equivalent to cj+cNj=cj+cmj=0subscript𝑐𝑗subscript𝑐𝑁𝑗subscript𝑐𝑗subscript𝑐𝑚𝑗0c_{j}+c_{N-j}=c_{j}+c_{m-j}=0italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_N - italic_j end_POSTSUBSCRIPT = italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_m - italic_j end_POSTSUBSCRIPT = 0, yielding cmj=cjsubscript𝑐𝑚𝑗subscript𝑐𝑗c_{m-j}=c_{-j}italic_c start_POSTSUBSCRIPT italic_m - italic_j end_POSTSUBSCRIPT = italic_c start_POSTSUBSCRIPT - italic_j end_POSTSUBSCRIPT. Consequently, all cjsubscript𝑐𝑗c_{j}italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT must equal to c0=0subscript𝑐00c_{0}=0italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0. For θ0=π/Nsubscript𝜃0𝜋𝑁\theta_{0}=\pi/Nitalic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_π / italic_N, we have

j=0N/21(cj+cN1j)eikxcosθj=0,superscriptsubscript𝑗0𝑁21subscript𝑐𝑗subscript𝑐𝑁1𝑗superscript𝑒𝑖𝑘𝑥subscript𝜃𝑗0\displaystyle\sum_{j=0}^{N/2-1}\left(c_{j}+c_{N-1-j}\right)e^{ikx\cos\theta_{j% }}=0,∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N / 2 - 1 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_N - 1 - italic_j end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_x roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = 0 , (A17)
c(m1)/2eikr+c(m+N1)/2eikr+subscript𝑐𝑚12superscript𝑒𝑖𝑘𝑟limit-fromsubscript𝑐𝑚𝑁12superscript𝑒𝑖𝑘𝑟\displaystyle c_{(m-1)/2}e^{ikr}+c_{(m+N-1)/2}e^{-ikr}+italic_c start_POSTSUBSCRIPT ( italic_m - 1 ) / 2 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT ( italic_m + italic_N - 1 ) / 2 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_i italic_k italic_r end_POSTSUPERSCRIPT +
j=0(m3)/2(cj+cm1j)eikrcos(θjα)+limit-fromsuperscriptsubscript𝑗0𝑚32subscript𝑐𝑗subscript𝑐𝑚1𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼\displaystyle\sum_{j=0}^{(m-3)/2}\left(c_{j}+c_{m-1-j}\right)e^{ikr\cos\left(% \theta_{j}-\alpha\right)}+∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m - 3 ) / 2 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_m - 1 - italic_j end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_α ) end_POSTSUPERSCRIPT +
j=m(m+N3)/2(cj+cm+N1j)eikrcos(θjα)=0.superscriptsubscript𝑗𝑚𝑚𝑁32subscript𝑐𝑗subscript𝑐𝑚𝑁1𝑗superscript𝑒𝑖𝑘𝑟subscript𝜃𝑗𝛼0\displaystyle\sum_{j=m}^{(m+N-3)/2}\left(c_{j}+c_{m+N-1-j}\right)e^{ikr\cos% \left(\theta_{j}-\alpha\right)}=0.∑ start_POSTSUBSCRIPT italic_j = italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m + italic_N - 3 ) / 2 end_POSTSUPERSCRIPT ( italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_m + italic_N - 1 - italic_j end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_r roman_cos ( italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_α ) end_POSTSUPERSCRIPT = 0 . (A18)

They are equivalent to cj+cNj1=cj+cmj1=0subscript𝑐𝑗subscript𝑐𝑁𝑗1subscript𝑐𝑗subscript𝑐𝑚𝑗10c_{j}+c_{N-j-1}=c_{j}+c_{m-j-1}=0italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_N - italic_j - 1 end_POSTSUBSCRIPT = italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT italic_m - italic_j - 1 end_POSTSUBSCRIPT = 0, yielding cmj1=c1jsubscript𝑐𝑚𝑗1subscript𝑐1𝑗c_{m-j-1}=c_{1-j}italic_c start_POSTSUBSCRIPT italic_m - italic_j - 1 end_POSTSUBSCRIPT = italic_c start_POSTSUBSCRIPT 1 - italic_j end_POSTSUBSCRIPT. Consequently, all cjsubscript𝑐𝑗c_{j}italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT must equal to c(m1)/2=0subscript𝑐𝑚120c_{(m-1)/2}=0italic_c start_POSTSUBSCRIPT ( italic_m - 1 ) / 2 end_POSTSUBSCRIPT = 0. So far, we have proved an absence of common eigenstates for the degenerate cases of θ0=0,π/Nsubscript𝜃00𝜋𝑁\theta_{0}=0,\pi/Nitalic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0 , italic_π / italic_N. Using the same tricks on Eq. (A10), we can obtain cj+m++cj=0superscriptsubscript𝑐𝑗𝑚superscriptsubscript𝑐𝑗0c_{j+m}^{+}+c_{j}^{-}=0italic_c start_POSTSUBSCRIPT italic_j + italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0. Applying the relation Eq. (A11) on the eigenstates |ψket𝜓|\psi\rangle| italic_ψ ⟩, we get

ψ(x,y)=j=0N/21sin(kxcosθj)sin(kysinθj).𝜓𝑥𝑦superscriptsubscript𝑗0𝑁21𝑘𝑥subscript𝜃𝑗𝑘𝑦subscript𝜃𝑗\psi(x,y)=\sum_{j=0}^{N/2-1}\sin\left(kx\cos\theta_{j}\right)\sin\left(ky\sin% \theta_{j}\right).italic_ψ ( italic_x , italic_y ) = ∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N / 2 - 1 end_POSTSUPERSCRIPT roman_sin ( italic_k italic_x roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_sin ( italic_k italic_y roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) . (A19)

One can verify ψ(x,y)𝜓𝑥𝑦\psi(x,y)italic_ψ ( italic_x , italic_y ) is invariant under the action of dihedral group DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT, θj2βπ/N±θjsubscript𝜃𝑗plus-or-minus2𝛽𝜋𝑁subscript𝜃𝑗\theta_{j}\to 2\beta\pi/N\pm\theta_{j}italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT → 2 italic_β italic_π / italic_N ± italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT for β𝛽\beta\in\mathbb{Z}italic_β ∈ blackboard_Z. ψ(x,y)𝜓𝑥𝑦\psi(x,y)italic_ψ ( italic_x , italic_y ) is a irreducible representation of group DNsubscript𝐷𝑁D_{N}italic_D start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT. We apply the boundary condition Eq. (A9) to Eq. (A19), yielding

j=0N/21sin(kLcosαcosθj)sin(kysinθj)=0.superscriptsubscript𝑗0𝑁21𝑘𝐿𝛼subscript𝜃𝑗𝑘𝑦subscript𝜃𝑗0\sum_{j=0}^{N/2-1}\sin\left(kL\cos\alpha\cos\theta_{j}\right)\sin\left(ky\sin% \theta_{j}\right)=0.∑ start_POSTSUBSCRIPT italic_j = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N / 2 - 1 end_POSTSUPERSCRIPT roman_sin ( italic_k italic_L roman_cos italic_α roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_sin ( italic_k italic_y roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) = 0 . (A20)

We can prove that Eq. (A20) is equivalent to

sin(kLcosαcosθj)=0orkLcosαcosθj=Mjπformulae-sequence𝑘𝐿𝛼subscript𝜃𝑗0or𝑘𝐿𝛼subscript𝜃𝑗subscript𝑀𝑗𝜋\sin\left(kL\cos\alpha\cos\theta_{j}\right)=0\quad\text{or}\quad kL\cos\alpha% \cos\theta_{j}=M_{j}\piroman_sin ( italic_k italic_L roman_cos italic_α roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) = 0 or italic_k italic_L roman_cos italic_α roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = italic_M start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_π (A21)

where Mj,j=0,1,,N/21formulae-sequencesubscript𝑀𝑗for-all𝑗01𝑁21M_{j}\in\mathbb{Z},\forall j=0,1,\ldots,N/2-1italic_M start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ∈ blackboard_Z , ∀ italic_j = 0 , 1 , … , italic_N / 2 - 1. This proof is similar to the proof of Eq. (A11). We take the 2i12𝑖12i-12 italic_i - 1-th (i=1,2,,N/2)𝑖12𝑁2(i=1,2,\ldots,N/2)( italic_i = 1 , 2 , … , italic_N / 2 ) deviatives of both side of Eq. (A15) and make y𝑦yitalic_y equal to zero. We obtain N1𝑁1N-1italic_N - 1 homogeneous linear equations by defining bj:=sin(kLcosαcosθj)sinθjassignsubscript𝑏𝑗𝑘𝐿𝛼subscript𝜃𝑗subscript𝜃𝑗b_{j}:=\sin\left(kL\cos\alpha\cos\theta_{j}\right)\sin\theta_{j}italic_b start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT := roman_sin ( italic_k italic_L roman_cos italic_α roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_sin italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT,

[111sin2θ0sin2θ1sin2θN1sin4θ0sin4θ1sin4θN1sinN2θ0sinN2θ1sinN2θN1][b0b1bN1]=0.delimited-[]111missing-subexpressionsuperscript2subscript𝜃0superscript2subscript𝜃1superscript2subscript𝜃𝑁1missing-subexpressionsuperscript4subscript𝜃0superscript4subscript𝜃1superscript4subscript𝜃𝑁1missing-subexpressionmissing-subexpressionsuperscript𝑁2subscript𝜃0superscript𝑁2subscript𝜃1superscript𝑁2subscript𝜃𝑁1missing-subexpressiondelimited-[]subscript𝑏0subscript𝑏1subscript𝑏𝑁10\left[\begin{array}[]{ccccc}1&1&\ldots&1\\ \sin^{2}\theta_{0}&\sin^{2}\theta_{1}&\ldots&\sin^{2}\theta_{N-1}\\ \sin^{4}\theta_{0}&\sin^{4}\theta_{1}&\ldots&\sin^{4}\theta_{N-1}\\ \vdots&\vdots&\ddots&\vdots\\ \sin^{N-2}\theta_{0}&\sin^{N-2}\theta_{1}&\ldots&\sin^{N-2}\theta_{N-1}\end{% array}\right]\left[\begin{array}[]{c}b_{0}\\ b_{1}\\ \vdots\\ b_{N-1}\end{array}\right]=0.[ start_ARRAY start_ROW start_CELL 1 end_CELL start_CELL 1 end_CELL start_CELL … end_CELL start_CELL 1 end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL start_CELL roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL … end_CELL start_CELL roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_sin start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL start_CELL roman_sin start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL … end_CELL start_CELL roman_sin start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL ⋮ end_CELL start_CELL ⋮ end_CELL start_CELL ⋱ end_CELL start_CELL ⋮ end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_sin start_POSTSUPERSCRIPT italic_N - 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL start_CELL roman_sin start_POSTSUPERSCRIPT italic_N - 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL … end_CELL start_CELL roman_sin start_POSTSUPERSCRIPT italic_N - 2 end_POSTSUPERSCRIPT italic_θ start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL start_CELL end_CELL end_ROW end_ARRAY ] [ start_ARRAY start_ROW start_CELL italic_b start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL ⋮ end_CELL end_ROW start_ROW start_CELL italic_b start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT end_CELL end_ROW end_ARRAY ] = 0 . (A22)

Its determinant is nonzero, so the solution is bj=0subscript𝑏𝑗0b_{j}=0italic_b start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = 0, yielding sin(KLcosαcosθj)=0𝐾𝐿𝛼subscript𝜃𝑗0\sin\left(KL\cos\alpha\cos\theta_{j}\right)=0roman_sin ( italic_K italic_L roman_cos italic_α roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) = 0. The relation Eq. (A21) can be used as the quantization condition. We divide it by itself with a different subscript, and obtain

cosθjcosθj=MjMj,j,j=0,1,,N1formulae-sequencesubscript𝜃𝑗subscript𝜃superscript𝑗subscript𝑀𝑗subscript𝑀superscript𝑗𝑗superscript𝑗01𝑁1\frac{\cos\theta_{j}}{\cos\theta_{j^{\prime}}}=\frac{M_{j}}{M_{j^{\prime}}}\in% \mathbb{Q},\quad j,j^{\prime}=0,1,\ldots,N-1divide start_ARG roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG start_ARG roman_cos italic_θ start_POSTSUBSCRIPT italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT end_ARG = divide start_ARG italic_M start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG start_ARG italic_M start_POSTSUBSCRIPT italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT end_ARG ∈ blackboard_Q , italic_j , italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 0 , 1 , … , italic_N - 1 (A23)

where Mj=MNj1subscript𝑀𝑗subscript𝑀𝑁𝑗1M_{j}=-M_{N-j-1}italic_M start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = - italic_M start_POSTSUBSCRIPT italic_N - italic_j - 1 end_POSTSUBSCRIPT for j=N/2,N/2+1,,N1𝑗𝑁2𝑁21𝑁1j=N/2,N/2+1,\ldots,N-1italic_j = italic_N / 2 , italic_N / 2 + 1 , … , italic_N - 1. By using the identity

cosθj+cosθj=2cosθj+θj2cosθjθj2subscript𝜃𝑗subscript𝜃superscript𝑗2subscript𝜃𝑗subscript𝜃superscript𝑗2subscript𝜃𝑗subscript𝜃superscript𝑗2\cos\theta_{j}+\cos\theta_{j^{\prime}}=2\cos\frac{\theta_{j}+\theta_{j^{\prime% }}}{2}\cos\frac{\theta_{j}-\theta_{j^{\prime}}}{2}roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + roman_cos italic_θ start_POSTSUBSCRIPT italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 2 roman_cos divide start_ARG italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_θ start_POSTSUBSCRIPT italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG roman_cos divide start_ARG italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_θ start_POSTSUBSCRIPT italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG (A24)

and setting jsuperscript𝑗j^{\prime}italic_j start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT to be j+2𝑗2j+2italic_j + 2, we find that

2cos2πN=cosθjcosθj+1+cosθj+2cosθj+1,22𝜋𝑁subscript𝜃𝑗subscript𝜃𝑗1subscript𝜃𝑗2subscript𝜃𝑗12\cos\frac{2\pi}{N}=\frac{\cos\theta_{j}}{\cos\theta_{j+1}}+\frac{\cos\theta_{% j+2}}{\cos\theta_{j+1}}\in\mathbb{Q},2 roman_cos divide start_ARG 2 italic_π end_ARG start_ARG italic_N end_ARG = divide start_ARG roman_cos italic_θ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG start_ARG roman_cos italic_θ start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT end_ARG + divide start_ARG roman_cos italic_θ start_POSTSUBSCRIPT italic_j + 2 end_POSTSUBSCRIPT end_ARG start_ARG roman_cos italic_θ start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT end_ARG ∈ blackboard_Q , (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 ±1,±1/2plus-or-minus1plus-or-minus12\pm 1,\pm 1/2± 1 , ± 1 / 2 and 00 (A proof was provided in [41]). They correspond to the right triangle billiards with α=π/6𝛼𝜋6\alpha=\pi/6italic_α = italic_π / 6 and π/4𝜋4\pi/4italic_π / 4, which are the only two complete integrable cases. For others, Eq. (A25) are not satisfied. Therefore, we conclude that in rational right triangle billiards, T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG 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.,

un(s)𝐧^(s)ψn(x(s),y(s)),subscriptsuperscript𝑢𝑛𝑠^𝐧𝑠subscript𝜓𝑛𝑥𝑠𝑦𝑠u^{\prime}_{n}(s)\triangleq\hat{\mathbf{n}}(s)\cdot\nabla\psi_{n}(x(s),y(s)),italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_s ) ≜ over^ start_ARG bold_n end_ARG ( italic_s ) ⋅ ∇ italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_x ( italic_s ) , italic_y ( italic_s ) ) , (B26)

where (x(s),y(s))𝑥𝑠𝑦𝑠(x(s),y(s))( italic_x ( italic_s ) , italic_y ( italic_s ) ) is a point on the boundary parameterized by the arc-length s𝑠sitalic_s and 𝐧^(s)^𝐧𝑠\hat{\mathbf{n}}(s)over^ start_ARG bold_n end_ARG ( italic_s ) denotes the outer normal unit vector at (x(s),y(s))𝑥𝑠𝑦𝑠(x(s),y(s))( italic_x ( italic_s ) , italic_y ( italic_s ) ). Conversely, the wave function can be obtained by boundary integral

ψn(x,y)=B𝑑sun(s)G[x,y;x(s),y(s)],subscript𝜓𝑛𝑥𝑦subscriptcontour-integral𝐵differential-d𝑠subscriptsuperscript𝑢𝑛𝑠𝐺𝑥𝑦𝑥𝑠𝑦𝑠\psi_{n}(x,y)=-\oint_{\partial B}ds\;u^{\prime}_{n}(s)G[x,y;x(s),y(s)],italic_ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_x , italic_y ) = - ∮ start_POSTSUBSCRIPT ∂ italic_B end_POSTSUBSCRIPT italic_d italic_s italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_s ) italic_G [ italic_x , italic_y ; italic_x ( italic_s ) , italic_y ( italic_s ) ] , (B27)

where B𝐵\partial B∂ italic_B is the billiard boundary curve, (x,y)𝑥𝑦(x,y)( italic_x , italic_y ) is the point inside the billiard table and G[x,y;x(s),y(s)]𝐺𝑥𝑦𝑥𝑠𝑦𝑠G[x,y;x(s),y(s)]italic_G [ italic_x , italic_y ; italic_x ( italic_s ) , italic_y ( italic_s ) ] is the 2-dimensional free particle Green’s function. In phase space, the conjugate variable to s𝑠sitalic_s is sinθ¯¯𝜃\sin\bar{\theta}roman_sin over¯ start_ARG italic_θ end_ARG, where θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG 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 un(θ¯)subscript𝑢𝑛¯𝜃u_{n}(\bar{\theta})italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) is expressed as

un(θ¯)=B𝑑sexp(ikssinθ¯)un(s),subscript𝑢𝑛¯𝜃subscriptcontour-integral𝐵differential-d𝑠𝑖𝑘𝑠¯𝜃subscriptsuperscript𝑢𝑛𝑠u_{n}(\bar{\theta})=\oint_{\partial B}ds\;\exp\left(iks\sin\bar{\theta}\right)% u^{\prime}_{n}(s),italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) = ∮ start_POSTSUBSCRIPT ∂ italic_B end_POSTSUBSCRIPT italic_d italic_s roman_exp ( italic_i italic_k italic_s roman_sin over¯ start_ARG italic_θ end_ARG ) italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_s ) , (B28)

where k𝑘kitalic_k is the wavenumber. It should be emphasized that the relations between θ𝜃\thetaitalic_θ and θ¯¯𝜃\bar{\theta}over¯ start_ARG italic_θ end_ARG have three kind depending on the sides: θ¯=θ+π/2¯𝜃𝜃𝜋2\bar{\theta}=\theta+\pi/2over¯ start_ARG italic_θ end_ARG = italic_θ + italic_π / 2 (for y=0𝑦0y=0italic_y = 0), θ¯=θ¯𝜃𝜃\bar{\theta}=\thetaover¯ start_ARG italic_θ end_ARG = italic_θ (for x=Lcosα𝑥𝐿𝛼x=L\cos\alphaitalic_x = italic_L roman_cos italic_α), θ¯=π/2+αθ¯𝜃𝜋2𝛼𝜃\bar{\theta}=\pi/2+\alpha-\thetaover¯ start_ARG italic_θ end_ARG = italic_π / 2 + italic_α - italic_θ (for y=xtanα𝑦𝑥𝛼y=x\tan\alphaitalic_y = italic_x roman_tan italic_α). Hence, |T(cosθ¯)|𝑇¯𝜃|T(\cos\bar{\theta})|| italic_T ( roman_cos over¯ start_ARG italic_θ end_ARG ) | remains the classical conserved quantity. Without ambiguity, we also use the symbol T𝑇Titalic_T to denote |T(cosθ¯)|𝑇¯𝜃|T(\cos\bar{\theta})|| italic_T ( roman_cos over¯ start_ARG italic_θ end_ARG ) |. The values of T𝑇Titalic_T can be used to distinguish different trajectories. For example, in π/8𝜋8\pi/8italic_π / 8 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 T=1𝑇1T=1italic_T = 1, while UDOs have various values of T𝑇Titalic_T. Quantum mechanically, there exists “superscar” states, exhibiting a superposition of spatially parallel POs [22, 38, 39]. T𝑇Titalic_T cannot be used to distinguish different superscar states as they all equal to 1111. But, the complete details of trajectories can be obtained by projecting the eigenstate using Husimi functions on the phase space expanded by s𝑠sitalic_s and sinθ¯¯𝜃\sin\bar{\theta}roman_sin over¯ start_ARG italic_θ end_ARG. The Husimi function h(q,p)𝑞𝑝h(q,p)italic_h ( italic_q , italic_p ) with q=s𝑞𝑠q=sitalic_q = italic_s, p=sinθ¯𝑝¯𝜃p=\sin\bar{\theta}italic_p = roman_sin over¯ start_ARG italic_θ end_ARG is defined as

h(q,p)1A|Bc(q,p)k(l)un(l)𝑑l|2,𝑞𝑝1𝐴superscriptsubscriptcontour-integral𝐵superscriptsubscript𝑐𝑞𝑝𝑘𝑙subscriptsuperscript𝑢𝑛𝑙differential-d𝑙2h(q,p)\triangleq\frac{1}{A}\left|\oint_{\partial B}c_{(q,p)}^{k}(l)u^{\prime}_% {n}(l)dl\right|^{2},italic_h ( italic_q , italic_p ) ≜ divide start_ARG 1 end_ARG start_ARG italic_A end_ARG | ∮ start_POSTSUBSCRIPT ∂ italic_B end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT ( italic_q , italic_p ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_l ) italic_u start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_l ) italic_d italic_l | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (B29)

where A𝐴Aitalic_A is the normalization parameter. c(q,p)k(l)superscriptsubscript𝑐𝑞𝑝𝑘𝑙c_{(q,p)}^{k}(l)italic_c start_POSTSUBSCRIPT ( italic_q , italic_p ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_l ) is the standard coherent state, i.e.,

c(q,p)k(l)=superscriptsubscript𝑐𝑞𝑝𝑘𝑙absent\displaystyle c_{(q,p)}^{k}(l)=italic_c start_POSTSUBSCRIPT ( italic_q , italic_p ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_l ) = jexp[ikp(lq+j)]subscript𝑗𝑖𝑘𝑝𝑙𝑞𝑗\displaystyle\sum_{j\in\mathbb{Z}}\exp[ikp(l-q+j\mathcal{L})]∑ start_POSTSUBSCRIPT italic_j ∈ blackboard_Z end_POSTSUBSCRIPT roman_exp [ italic_i italic_k italic_p ( italic_l - italic_q + italic_j caligraphic_L ) ] (B30)
×exp[k2(lq+j)2]absent𝑘2superscript𝑙𝑞𝑗2\displaystyle\times\exp\left[-\frac{k}{2}(l-q+j\mathcal{L})^{2}\right]× roman_exp [ - divide start_ARG italic_k end_ARG start_ARG 2 end_ARG ( italic_l - italic_q + italic_j caligraphic_L ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ]

where \mathcal{L}caligraphic_L is the length of boundary and the sum of j𝑗jitalic_j ensures the coherent states being periodic with a period of \mathcal{L}caligraphic_L.

Refer to caption
Figure 7: Trajectories of typical eigenstates for π/8𝜋8\pi/8italic_π / 8 billiard. The left column are the distributions of Husimi functions h(s,sinθ¯)𝑠¯𝜃h(s,\sin\bar{\theta})italic_h ( italic_s , roman_sin over¯ start_ARG italic_θ end_ARG ) on the phase space spanned by s𝑠sitalic_s and sinθ¯¯𝜃\sin\bar{\theta}roman_sin over¯ start_ARG italic_θ end_ARG. c𝑐citalic_c, b𝑏bitalic_b and a𝑎aitalic_a respectively represent the hypotenuse, the adjacent side and the opposite side of π/8𝜋8\pi/8italic_π / 8 corner. The right column are the conjugate boundary function un(θ¯)subscript𝑢𝑛¯𝜃u_{n}(\bar{\theta})italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) with respect to the conserved quantity |T(cosθ¯)|𝑇¯𝜃\left|T(\cos\bar{\theta})\right|| italic_T ( roman_cos over¯ start_ARG italic_θ end_ARG ) |. Each row corresponds to the same eigenstate. (a) A superscar state with NWeyl10017.3subscript𝑁Weyl10017.3N_{\operatorname{Weyl}}\approx 10017.3italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 10017.3 and σ0.006𝜎0.006\sigma\approx 0.006italic_σ ≈ 0.006, typically exhibiting a superposition of parallel periodic orbits. (b) A spacial uniformly distributed orbit with NWeyl40090.0subscript𝑁Weyl40090.0N_{\operatorname{Weyl}}\approx 40090.0italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 40090.0 and σ0.002𝜎0.002\sigma\approx 0.002italic_σ ≈ 0.002. (c1) A superposition of a single superscar and a single UDO with NWeyl60144.6subscript𝑁Weyl60144.6N_{\operatorname{Weyl}}\approx 60144.6italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 60144.6 and σ0.563𝜎0.563\sigma\approx 0.563italic_σ ≈ 0.563. (c2) is the eigenstate next to (c1) with NWeyl60144.1subscript𝑁Weyl60144.1N_{\operatorname{Weyl}}\approx 60144.1italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 60144.1 and σ0.489𝜎0.489\sigma\approx 0.489italic_σ ≈ 0.489. (d1) A superposition of two UDOs with NWeyl61476.4subscript𝑁Weyl61476.4N_{\operatorname{Weyl}}\approx 61476.4italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 61476.4 and σ0.461𝜎0.461\sigma\approx 0.461italic_σ ≈ 0.461. (d2) is the eigenstate next to (d1) with NWeyl61476.6subscript𝑁Weyl61476.6N_{\operatorname{Weyl}}\approx 61476.6italic_N start_POSTSUBSCRIPT roman_Weyl end_POSTSUBSCRIPT ≈ 61476.6 and σ0.509𝜎0.509\sigma\approx 0.509italic_σ ≈ 0.509.

The results for π/8𝜋8\pi/8italic_π / 8 billiard are shown in Fig. 7. And the Poincaré-Husimi representations and the plot of the conjugate boundary function un(θ¯)subscript𝑢𝑛¯𝜃u_{n}(\bar{\theta})italic_u start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( over¯ start_ARG italic_θ end_ARG ) with respect to the conserved quantity T𝑇Titalic_T are compared to illustrate the trajectories of typical eigenstates with different σ𝜎\sigmaitalic_σ (the uncertainty of T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG). 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 T𝑇Titalic_T and thus possessing the conserved quantity T^^𝑇\hat{T}over^ start_ARG italic_T end_ARG. The eigenstate in Fig. 7.(c1) is a superposition of one superscar (including four parallel POs) and one UDO from the observation on h(s,sinθ¯)𝑠¯𝜃h(s,\sin\bar{\theta})italic_h ( italic_s , roman_sin over¯ start_ARG italic_θ end_ARG ). Consequently, it has two separated peaks on the distribution of T𝑇Titalic_T, causing an elevated σ𝜎\sigmaitalic_σ. 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 1/fα1superscript𝑓𝛼1/f^{\alpha}1 / italic_f start_POSTSUPERSCRIPT italic_α end_POSTSUPERSCRIPT 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).