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:nlin/0612011 · New Journal of Physics

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

G Veble, T Prosen, M Robnik

Published Jan 25, 2007• 39 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1088/1367-2630/9/1/015

Related papers

Constraint operator solution to quantum billiard problems

David A. McGrew, W. Bauer

We introduce an additional method to solve Schr\"odinger's equation for a free particle in an infinite well of arbitrary shape (the Helmholtz equation with Dirichlet boundary conditions), a problem of interest in the area of quantum chaos. We expand the wave function in a basis of products of sine functions, then use the constraint operator to contain the wave function to a region within the domain of the basis functions. In this manner, a quantum billiard problem of arbitrary shape can be solved. Several methods exist to solve problems of this sort, but as recent work reviewing these methods has shown, all have shortcomings. Our work represents a different direction in the solution of these problems. Our method is different in that it provides a means of computing an eigenbasis. It is also interesting from a physical standpoint in that it can represent the Hamiltonian of a classically chaotic system in the basis of a classically regular system. \textcopyright{} 1996 The American Physical Society.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Nov 1, 1996•4 citations

Boundary Integral Method for Stationary States of Two-Dimensional Quantum Systems

Ioan Kosztin, Klaus J. Schulten

The boundary integral method for calculating the stationary states of a quantum particle in nano-devices and quantum billiards is presented in detail at an elementary level. According to the method, wave functions inside the domain of the device or billiard are expressed in terms of line integrals of the wavefunction and its normal derivative along the domain's boundary; the respective energy eigenvalues are obtained as the roots of Fredholm determinants. Numerical implementations of the method are described and applied to determine the energy level statistics of billiards with circular and stadium shapes and demonstrate the quantum mechanical characteristics of chaotic motion. The treatment of other examples as well as the advantages and limitations of the boundary integral method are discussed.

arXiv:physics/9702022•International Journal of Modern Physics C•Apr 1, 1997•42 citations

Expansion method for stationary states of quantum billiards

David L. Kaufman, Ioan Kosztin, Klaus J. Schulten

A simple expansion method for numerically calculating the energy levels and the corresponding wave functions of a quantum particle in a two-dimensional infinite potential well with arbitrary shape (quantum billiard) is presented. The method permits the study of quantum billiards in an introductory quantum mechanics course. According to the method, wave functions inside the billiard are expressed in terms of an expansion of a complete set of orthonormal functions defined in a surrounding rectangle for which the Dirichlet boundary conditions apply, while approximating the billiard boundary by a potential energy step of a sufficiently large size. Numerical implementations of the method are described and applied to determine the energies and wave functions for quarter-circle, circle, and triangle billiards. Finally, the expansion method is applied to investigate the quantum signatures of chaos in a classically chaotic generic-triangle billiard.

American Journal of Physics•Feb 1, 1999•56 citations

Scarring of Dirac fermions in chaotic billiards

Xuan Ni, Liang Huang, Ying‐Cheng Lai +1 more

Scarring in quantum systems with classical chaotic dynamics is one of the most remarkable phenomena in modern physics. Previous works were concerned mostly with nonrelativistic quantum systems described by the Schrödinger equation. The question remains outstanding of whether truly relativistic quantum particles that obey the Dirac equation can scar. A significant challenge is the lack of a general method for solving the Dirac equation in closed domains of arbitrary shape. In this paper, we develop a numerical framework for obtaining complete eigensolutions of massless fermions in general two-dimensional confining geometries. The key ingredients of our method are the proper handling of the boundary conditions and an efficient discretization scheme that casts the original equation in a matrix representation. The method is validated by (1) comparing the numerical solutions to analytic results for a geometrically simple confinement and (2) verifying that the calculated energy level-spacing statistics of integrable and chaotic geometries agree with the known results. Solutions of the Dirac equation in a number of representative chaotic geometries establish firmly the existence of scarring of Dirac fermions.

Physical Review E•Jul 11, 2012•25 citations

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

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

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.

nlin.CDCross-listedarXiv:2605.27693•May 28, 2026

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
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.
We introduce an additional method to solve Schr\"odinger's equation for a free particle in an infinite well of arbitrary shape (the Helmholtz equation with Dirichlet boundary conditions), a problem of interest in the area of quantum chaos. We expand the wave function in a basis of products of sine functions, then use the constraint operator to contain the wave function to a region within the domain of the basis functions. In this manner, a quantum billiard problem of arbitrary shape can be solved. Several methods exist to solve problems of this sort, but as recent work reviewing these methods has shown, all have shortcomings. Our work represents a different direction in the solution of these problems. Our method is different in that it provides a means of computing an eigenbasis. It is also interesting from a physical standpoint in that it can represent the Hamiltonian of a classically chaotic system in the basis of a classically regular system. \textcopyright{} 1996 The American Physical Society.
The boundary integral method for calculating the stationary states of a quantum particle in nano-devices and quantum billiards is presented in detail at an elementary level. According to the method, wave functions inside the domain of the device or billiard are expressed in terms of line integrals of the wavefunction and its normal derivative along the domain's boundary; the respective energy eigenvalues are obtained as the roots of Fredholm determinants. Numerical implementations of the method are described and applied to determine the energy level statistics of billiards with circular and stadium shapes and demonstrate the quantum mechanical characteristics of chaotic motion. The treatment of other examples as well as the advantages and limitations of the boundary integral method are discussed.
A simple expansion method for numerically calculating the energy levels and the corresponding wave functions of a quantum particle in a two-dimensional infinite potential well with arbitrary shape (quantum billiard) is presented. The method permits the study of quantum billiards in an introductory quantum mechanics course. According to the method, wave functions inside the billiard are expressed in terms of an expansion of a complete set of orthonormal functions defined in a surrounding rectangle for which the Dirichlet boundary conditions apply, while approximating the billiard boundary by a potential energy step of a sufficiently large size. Numerical implementations of the method are described and applied to determine the energies and wave functions for quarter-circle, circle, and triangle billiards. Finally, the expansion method is applied to investigate the quantum signatures of chaos in a classically chaotic generic-triangle billiard.
Scarring in quantum systems with classical chaotic dynamics is one of the most remarkable phenomena in modern physics. Previous works were concerned mostly with nonrelativistic quantum systems described by the Schrödinger equation. The question remains outstanding of whether truly relativistic quantum particles that obey the Dirac equation can scar. A significant challenge is the lack of a general method for solving the Dirac equation in closed domains of arbitrary shape. In this paper, we develop a numerical framework for obtaining complete eigensolutions of massless fermions in general two-dimensional confining geometries. The key ingredients of our method are the proper handling of the boundary conditions and an efficient discretization scheme that casts the original equation in a matrix representation. The method is validated by (1) comparing the numerical solutions to analytic results for a geometrically simple confinement and (2) verifying that the calculated energy level-spacing statistics of integrable and chaotic geometries agree with the known results. Solutions of the Dirac equation in a number of representative chaotic geometries establish firmly the existence of scarring of Dirac fermions.
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.
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.