Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2605.27693

Effects of stickiness in quantum chaotic billiards with nnn-fold symmetry

R. B. do Carmo, T. Araújo Lima

nlin.CDquant-phSubmitted May 26, 2026• Announced May 28, 2026
arXivHTMLPDF

Abstract

In this work, we study a family of fully chaotic billiards that exhibits only rotational symmetries, whose geometry is based on the $C_3$ symmetry system proposed by Leyvraz, Schmit, and Seligman~(LSS) in 1996. Quantum spectral analyses are performed on billiards with symmetry $C_n$~(the billiard repeats itself under rotations of $2π/n$), where $n$ is the symmetry parameter. In these systems, there are subspaces of singlets~(invariant under time reversal) and doublets~(not invariant under time reversal). For the LSS billiard, it has been established both numerically and experimentally that the corresponding subspectra follow the Gaussian Orthogonal Ensemble~(GOE) statistics for singlets and the Gaussian Unitary Ensemble~(GUE) statistics for doublets. From a classical perspective, the shapes of these billiards allow certain subregions of phase space to be visited more frequently by chaotic trajectories, a phenomenon known as stickiness. We investigate the relationship between the fraction of sticky regions in classical phase space and the deviations of the energy subspectra from GOE and GUE statistics. Our results suggest the existence of correlations between the energy distributions associated with different symmetry subspaces. In addition, we discuss aspects related to the superposition of the different energy subspectra.

Comments
8 pages, 5 figures

Related papers

Spectral statistics and localization properties of a $C_3$-symmetric billiard

Matic Orel, Marko Robnik

We revisit the spectral statistics of the C$_3$--symmetric billiard introduced by Dembowski [Phys. Rev. E, R4516 (2000)], which exhibits both GOE and GUE statistics depending on the symmetry block. Using high--precision Beyn's contour--integral method for the nonlinear Fredholm eigenvalue problem with built-in separation of irreducible subspaces, we compute 2.8x10$^5$ eigenvalues in each symmetry subspace, enabling statistically meaningful comparisons with random matrix theory. The improved spectra reveal clear GOE--GUE correspondence and resolve previously observed deviations in long--range spectral correlations. Furthermore, we analyze phase--space eigenstate localization through the distribution of entropy localization measures, which, for chaotic states follow a Beta distribution whose standard deviation decays as a power--law with energy, consistent with the onset of quantum ergodicity as described by Schnirelman's theorem.

quant-pharXiv:2603.03460•Mar 5, 2026

Expanded boundary integral method and chaotic time-reversal doublets in quantum billiards

G Veble, T Prosen, M Robnik

We present the expanded boundary integral method for solving the planar Helmholtz problem, which combines the ideas of the boundary integral method and the scaling method and is applicable to arbitrary shapes. We apply the method to a chaotic billiard with unidirectional transport, where we demonstrate the existence of doublets of chaotic eigenstates, which are quasi-degenerate due to time-reversal symmetry, and a very particular level spacing distribution that attains a chaotic Shnirelman peak at short energy ranges and exhibits Gaussian Unitary Ensemble (GUE) like statistics for large energy ranges. We show that, as a consequence of such particular level statistics or algebraic tunnelling between disjoint chaotic components connected by time-reversal operation, the system exhibits quantum current reversals.

arXiv:nlin/0612011•New Journal of Physics•Jan 25, 2007•39 citations

Test of a numerical approach to the quantization of billiards

Barbara Dietz, A. Heine, Vincent Heuveline +1 more

A method for computing large numbers of eigenvalues of self-adjoint elliptic operators [J. Comput. Phys. 184, 321 (2003)] is tested in numerical studies of the spectral properties of quantum billiards. To this extent, we study a time-reversal invariant quantum billiard of threefold symmetry, that undergoes a transformation in its symmetry properties from C(3v) to C3 . Thereby a transition from Gaussian orthogonal to Gaussian unitary ensemble statistics is observed, verifying earlier experimental indications and theoretical predictions. At the same time our numerical ansatz is shown to be applicable to arbitrary billiard shapes.

Physical Review E•Feb 24, 2005•16 citations

Effects of stickiness on the quantum states of strongly chaotic open systems

Miguel A. Prado Reynoso, Edson M. Signor, Sandra D. Prado +1 more

Focusing on the standard map of the kicked rotor, we investigate the effects of classical stickiness (orbits temporarily confined to a region of a chaotic phase space) on an open quantum system. The regions of stickiness that we identify emerge in the strongly chaotic map. By scanning the system's phase space with a leak, we show that stickiness reduces the degree of delocalization of the quantum states and the decay rate of generic initial states. We find an excellent correspondence between the classical dwell time and finite-time Lyapunov exponents with the quantum dwell time and Wehrl entropies of the quantum states. Knowledge of the structure of the classically chaotic trajectories can thus be used to determine the best placement of leaks to enhance or decrease the degree of delocalization of quantum states and to modify the quantum dynamics of open systems.

arXiv:2407.08780•Physical review. E•Dec 27, 2024•3 citations

Geometry of quantum states in random matrix ensembles and the chaos-integrability transition

Anonymous, Keun-young Kim A, Kunal Pal +1 more

We consider the geometry of quantum states associated with random matrix Hamiltonians belonging to ensembles that exhibit an integrable-to-chaotic transition in terms of the nearest-neighbour energy level spacing distribution, focusing on the $β$-Gaussian ensembles with generic Dyson index. For the tridiagonal Gaussian $β$-ensemble, which shows chaos-to-integrability transition with varying Dyson index, we first calculate an analytical expression for the ensemble-averaged fidelity susceptibility for a two-by-two matrix representation of the Hamiltonian for generic values of the Dyson index, and show that it diverges as the total Hamiltonian, which is the sum of a diagonal matrix with independent elements and a tridiagonal $β$-matrix, goes over to the integrable phase. Next, for large-dimensional matrices, we numerically compute the fidelity susceptibility and the quantum metric tensor, respectively, for a one-parameter and a two-parameter class of Hamiltonians that we construct from the $β$-ensemble, and find scaling relations of these quantities for chaotic and integrable phases. We also consider variations of the $β$-ensemble, such as the one that preserves the rotational invariance of the ensemble, and compute the relevant ensemble-averaged fidelity susceptibility to show that it has similar features as the tridiagonal ensemble in the integrable as well as the chaotic phase, thereby establishing the universality of these properties. Finally, the presence of non-vanishing non-diagonal elements, which arises due to the rotational non-invariance of the $β$-ensembles, is a specific feature of the corresponding metric tensor, and we use this to illuminate the difference between the quantum state space geometry of these ensembles and the rotationally invariant ones, such as the classical Gaussian ensembles.

arXiv:2507.13067•Physical review. D/Physical review. D.•Aug 6, 2026•0 citations

Quantum signatures of chaos in relativistic quantum billiards with shapes of circle- and ellipse-sectors*

Pei Yu, Weihua Zhang, Barbara Dietz +1 more

Abstract According to the Berry–Tabor conjecture, the spectral properties of typical nonrelativistic quantum systems with an integrable classical counterpart agree with those of Poissonian random numbers. We investigate to what extend it applies to relativistic neutrino billiards (NBs) consisting of a spin-1/2 particle confined to a bounded planar domain by imposing suitable boundary conditions (BCs). In distinction to nonrelativistic quantum billiards (QBs), NBs do not have a well-defined classical counterpart. However, the peaks in the length spectra, that is, the modulus of the Fourier transform of the spectral density from wave number to length, of NBs are just like for QBs at the lengths of periodic orbits of the classical billiard (CB). This implies that there must be a connection between NBs and the dynamic of the CB. We demonstrate that NBs with shapes of circle- and ellipse-sectors with an integrable classical dynamic, obtained by cutting the circle and ellipse NB along symmetry lines, have no common eigenstates with the latter and that, indeed, their spectral properties can be similar to those of classically chaotic QBs. These features orginate from the intermingling of symmetries of the spinor components and the discontinuity in the BCs leading to contradictory conditional equations at corners connecting curved and straight boundary parts. To corroborate the necessity of the curved boundary part in order to generate GOE-like behavior, we furthermore consider the right-angled triangle NB constructed by halving the equilateral-triangle NB along a symmetry axis. For an understanding of these findings in terms of purely classical quantities we use the semiclassical approach recently developed for massive NBs, and Poincaré–Husimi distributions of the eigenstates in classical phase space. The results indicate, that in the ultrarelativistic limit these NBs do not show the behavior expected for classically chaotic QBs.

Journal of Physics A Mathematical and Theoretical•Apr 19, 2022•10 citations
In this work, we study a family of fully chaotic billiards that exhibits only rotational symmetries, whose geometry is based on the C3C_3C3​ symmetry system proposed by Leyvraz, Schmit, and Seligman~(LSS) in 1996. Quantum spectral analyses are performed on billiards with symmetry CnC_nCn​(the billiard repeats itself under rotations of 2π/n2π/n2π/n), where nnn is the symmetry parameter. In these systems, there are subspaces of singlets(invariant under time reversal) and doublets~(not invariant under time reversal). For the LSS billiard, it has been established both numerically and experimentally that the corresponding subspectra follow the Gaussian Orthogonal Ensemble~(GOE) statistics for singlets and the Gaussian Unitary Ensemble~(GUE) statistics for doublets. From a classical perspective, the shapes of these billiards allow certain subregions of phase space to be visited more frequently by chaotic trajectories, a phenomenon known as stickiness. We investigate the relationship between the fraction of sticky regions in classical phase space and the deviations of the energy subspectra from GOE and GUE statistics. Our results suggest the existence of correlations between the energy distributions associated with different symmetry subspaces. In addition, we discuss aspects related to the superposition of the different energy subspectra.
Spectral statistics and localization properties of a C3C_3C3​-symmetric billiard
We revisit the spectral statistics of the C3_33​--symmetric billiard introduced by Dembowski [Phys. Rev. E, R4516 (2000)], which exhibits both GOE and GUE statistics depending on the symmetry block. Using high--precision Beyn's contour--integral method for the nonlinear Fredholm eigenvalue problem with built-in separation of irreducible subspaces, we compute 2.8x105^55 eigenvalues in each symmetry subspace, enabling statistically meaningful comparisons with random matrix theory. The improved spectra reveal clear GOE--GUE correspondence and resolve previously observed deviations in long--range spectral correlations. Furthermore, we analyze phase--space eigenstate localization through the distribution of entropy localization measures, which, for chaotic states follow a Beta distribution whose standard deviation decays as a power--law with energy, consistent with the onset of quantum ergodicity as described by Schnirelman's theorem.
We present the expanded boundary integral method for solving the planar Helmholtz problem, which combines the ideas of the boundary integral method and the scaling method and is applicable to arbitrary shapes. We apply the method to a chaotic billiard with unidirectional transport, where we demonstrate the existence of doublets of chaotic eigenstates, which are quasi-degenerate due to time-reversal symmetry, and a very particular level spacing distribution that attains a chaotic Shnirelman peak at short energy ranges and exhibits Gaussian Unitary Ensemble (GUE) like statistics for large energy ranges. We show that, as a consequence of such particular level statistics or algebraic tunnelling between disjoint chaotic components connected by time-reversal operation, the system exhibits quantum current reversals.
A method for computing large numbers of eigenvalues of self-adjoint elliptic operators [J. Comput. Phys. 184, 321 (2003)] is tested in numerical studies of the spectral properties of quantum billiards. To this extent, we study a time-reversal invariant quantum billiard of threefold symmetry, that undergoes a transformation in its symmetry properties from C(3v) to C3 . Thereby a transition from Gaussian orthogonal to Gaussian unitary ensemble statistics is observed, verifying earlier experimental indications and theoretical predictions. At the same time our numerical ansatz is shown to be applicable to arbitrary billiard shapes.
Focusing on the standard map of the kicked rotor, we investigate the effects of classical stickiness (orbits temporarily confined to a region of a chaotic phase space) on an open quantum system. The regions of stickiness that we identify emerge in the strongly chaotic map. By scanning the system's phase space with a leak, we show that stickiness reduces the degree of delocalization of the quantum states and the decay rate of generic initial states. We find an excellent correspondence between the classical dwell time and finite-time Lyapunov exponents with the quantum dwell time and Wehrl entropies of the quantum states. Knowledge of the structure of the classically chaotic trajectories can thus be used to determine the best placement of leaks to enhance or decrease the degree of delocalization of quantum states and to modify the quantum dynamics of open systems.
We consider the geometry of quantum states associated with random matrix Hamiltonians belonging to ensembles that exhibit an integrable-to-chaotic transition in terms of the nearest-neighbour energy level spacing distribution, focusing on the βββ-Gaussian ensembles with generic Dyson index. For the tridiagonal Gaussian βββ-ensemble, which shows chaos-to-integrability transition with varying Dyson index, we first calculate an analytical expression for the ensemble-averaged fidelity susceptibility for a two-by-two matrix representation of the Hamiltonian for generic values of the Dyson index, and show that it diverges as the total Hamiltonian, which is the sum of a diagonal matrix with independent elements and a tridiagonal βββ-matrix, goes over to the integrable phase. Next, for large-dimensional matrices, we numerically compute the fidelity susceptibility and the quantum metric tensor, respectively, for a one-parameter and a two-parameter class of Hamiltonians that we construct from the βββ-ensemble, and find scaling relations of these quantities for chaotic and integrable phases. We also consider variations of the βββ-ensemble, such as the one that preserves the rotational invariance of the ensemble, and compute the relevant ensemble-averaged fidelity susceptibility to show that it has similar features as the tridiagonal ensemble in the integrable as well as the chaotic phase, thereby establishing the universality of these properties. Finally, the presence of non-vanishing non-diagonal elements, which arises due to the rotational non-invariance of the βββ-ensembles, is a specific feature of the corresponding metric tensor, and we use this to illuminate the difference between the quantum state space geometry of these ensembles and the rotationally invariant ones, such as the classical Gaussian ensembles.
Abstract According to the Berry–Tabor conjecture, the spectral properties of typical nonrelativistic quantum systems with an integrable classical counterpart agree with those of Poissonian random numbers. We investigate to what extend it applies to relativistic neutrino billiards (NBs) consisting of a spin-1/2 particle confined to a bounded planar domain by imposing suitable boundary conditions (BCs). In distinction to nonrelativistic quantum billiards (QBs), NBs do not have a well-defined classical counterpart. However, the peaks in the length spectra, that is, the modulus of the Fourier transform of the spectral density from wave number to length, of NBs are just like for QBs at the lengths of periodic orbits of the classical billiard (CB). This implies that there must be a connection between NBs and the dynamic of the CB. We demonstrate that NBs with shapes of circle- and ellipse-sectors with an integrable classical dynamic, obtained by cutting the circle and ellipse NB along symmetry lines, have no common eigenstates with the latter and that, indeed, their spectral properties can be similar to those of classically chaotic QBs. These features orginate from the intermingling of symmetries of the spinor components and the discontinuity in the BCs leading to contradictory conditional equations at corners connecting curved and straight boundary parts. To corroborate the necessity of the curved boundary part in order to generate GOE-like behavior, we furthermore consider the right-angled triangle NB constructed by halving the equilateral-triangle NB along a symmetry axis. For an understanding of these findings in terms of purely classical quantities we use the semiclassical approach recently developed for massive NBs, and Poincaré–Husimi distributions of the eigenstates in classical phase space. The results indicate, that in the ultrarelativistic limit these NBs do not show the behavior expected for classically chaotic QBs.