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

American Journal of Physics

Expansion method for stationary states of quantum billiards

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

Published Feb 1, 1999• 56 citations
Publisher

Abstract

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.

DOI
10.1119/1.19208

Related papers

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

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

Wave packet construction in two-dimensional quantum billiards: Blueprints for the square, equilateral triangle, and circular cases

M. A. Doncheski, S. Heppelmann, R. W. Robinett +1 more

We present quasianalytical and numerical calculations of Gaussian wave packet solutions of the Schrödinger equation for two-dimensional infinite well and quantum billiard problems with equilateral triangle, square, and circular footprints. These cases correspond to N=3, N=4, and N→∞ regular polygonal billiards and infinite wells, respectively. In each case the energy eigenvalues and wave functions are given in terms of familiar special functions. For the first two systems, we obtain closed form expressions for the expansion coefficients for localized Gaussian wave packets in terms of the eigenstates of the particular geometry. For the circular case, we discuss numerical approaches. We use these results to discuss the short-time, quasiclassical evolution in these geometries and the structure of wave packet revivals. We also show how related half-well problems can be easily solved in each of the three cases.

arXiv:quant-ph/0307070•American Journal of Physics•May 12, 2003•33 citations

Wave packet construction in two-dimensional quantum billiards: Blueprints for the square, equilateral triangle, and circular cases

M. A. Doncheski, S. Heppelmann, Richard W. Robinett +1 more

We present quasi-analytical and numerical calculations of Gaussian wave packet solutions of the Schr\\"odinger equation for two-dimensional infinite well and quantum billiard problems with equilateral triangle, square, and circular footprints. These cases correspond to N=3, N=4, and $N \\to \\infty$ regular polygonal billiards and infinite wells, respectively. In each case the energy eigenvalues and wavefunctions are given in terms of familiar special functions. For the first two systems, we obtain closed form expressions for the expansion coefficients for localized Gaussian wavepackets in terms of the eigenstates of the particular geometry. For the circular case, we discuss numerical approaches. We use these results to discuss the short-time, quasi-classical evolution in these geometries and the structure of wave packet revivals. We also show how related half-well problems can be easily solved in each of the three cases.

•Jul 9, 2003•20 citations

Modeling Quantum Billiards with the Finite Element Method: Searching for Quantum Scarring Candidates

Daniel Pierce, Renuka Rajapakse

An electron in quantum confinement takes on a discrete energy spectrum which is defined based on the solution to the Schrodinger Equation for a given potential. Well defined closed-form energy spectra are known for the particle in a box, circular potential, quarter circle potential, and an equilateral triangle. A closed-form solution for more complex shapes may not be known, but numerical methods can be used to find an approximate solution. In this research, an application of the Finite Element Method (FEM) in Wolfram Mathematica is presented and applied to Quantum Billiards with a variety of geometries. To assess the accuracy of the method, the computed energy states are analyzed in the limit of a polygon with an increasing number of sides, the numerical results are validated against analytical solutions for geometries with known exact forms, and a standard convergence test is conducted. The FEM results closely match analytical solutions for known potentials, demonstrating its high accuracy. For high energy index n, quantum scarring may emerge for certain geometries. The nature of quantum scarring and its presence in the computed models is also investigated qualitatively.

quant-pharXiv:2603.24864•Mar 27, 2026

Behaviour of boundary functions for quantum billiards

A B cker, S F rstberger, R Schubert +1 more

We study the behaviour of the normal derivative of eigenfunctions of the Helmholtz equation inside billiards with Dirichlet boundary condition. These boundary functions are of particular importance because they uniquely determine the eigenfunctions inside the billiard and also other physical quantities of interest. Therefore, they form a reduced representation of the quantum system, analogous to the Poincaré section of the classical system. For the normal derivatives we introduce an equivalent to the standard Green function and derive an integral equation on the boundary. Based on this integral equation we compute the first two terms of the mean asymptotic behaviour of the boundary functions for large energies. The first term is universal and independent of the shape of the billiard. The second one is proportional to the curvature of the boundary. The asymptotic behaviour is compared with numerical results for the stadium billiard, different limaçon billiards and the circle billiard, and good agreement is found. Furthermore, we derive an asymptotic completeness relation for the boundary functions.

arXiv:nlin/0207034•Journal of Physics A Mathematical and General•Nov 19, 2002•21 citations
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.
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.
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.
We present quasianalytical and numerical calculations of Gaussian wave packet solutions of the Schrödinger equation for two-dimensional infinite well and quantum billiard problems with equilateral triangle, square, and circular footprints. These cases correspond to N=3, N=4, and N→∞ regular polygonal billiards and infinite wells, respectively. In each case the energy eigenvalues and wave functions are given in terms of familiar special functions. For the first two systems, we obtain closed form expressions for the expansion coefficients for localized Gaussian wave packets in terms of the eigenstates of the particular geometry. For the circular case, we discuss numerical approaches. We use these results to discuss the short-time, quasiclassical evolution in these geometries and the structure of wave packet revivals. We also show how related half-well problems can be easily solved in each of the three cases.
We present quasi-analytical and numerical calculations of Gaussian wave packet solutions of the Schr\\"odinger equation for two-dimensional infinite well and quantum billiard problems with equilateral triangle, square, and circular footprints. These cases correspond to N=3, N=4, and NtoinftyN \\to \\inftyNtoinfty regular polygonal billiards and infinite wells, respectively. In each case the energy eigenvalues and wavefunctions are given in terms of familiar special functions. For the first two systems, we obtain closed form expressions for the expansion coefficients for localized Gaussian wavepackets in terms of the eigenstates of the particular geometry. For the circular case, we discuss numerical approaches. We use these results to discuss the short-time, quasi-classical evolution in these geometries and the structure of wave packet revivals. We also show how related half-well problems can be easily solved in each of the three cases.
An electron in quantum confinement takes on a discrete energy spectrum which is defined based on the solution to the Schrodinger Equation for a given potential. Well defined closed-form energy spectra are known for the particle in a box, circular potential, quarter circle potential, and an equilateral triangle. A closed-form solution for more complex shapes may not be known, but numerical methods can be used to find an approximate solution. In this research, an application of the Finite Element Method (FEM) in Wolfram Mathematica is presented and applied to Quantum Billiards with a variety of geometries. To assess the accuracy of the method, the computed energy states are analyzed in the limit of a polygon with an increasing number of sides, the numerical results are validated against analytical solutions for geometries with known exact forms, and a standard convergence test is conducted. The FEM results closely match analytical solutions for known potentials, demonstrating its high accuracy. For high energy index n, quantum scarring may emerge for certain geometries. The nature of quantum scarring and its presence in the computed models is also investigated qualitatively.
We study the behaviour of the normal derivative of eigenfunctions of the Helmholtz equation inside billiards with Dirichlet boundary condition. These boundary functions are of particular importance because they uniquely determine the eigenfunctions inside the billiard and also other physical quantities of interest. Therefore, they form a reduced representation of the quantum system, analogous to the Poincaré section of the classical system. For the normal derivatives we introduce an equivalent to the standard Green function and derive an integral equation on the boundary. Based on this integral equation we compute the first two terms of the mean asymptotic behaviour of the boundary functions for large energies. The first term is universal and independent of the shape of the billiard. The second one is proportional to the curvature of the boundary. The asymptotic behaviour is compared with numerical results for the stadium billiard, different limaçon billiards and the circle billiard, and good agreement is found. Furthermore, we derive an asymptotic completeness relation for the boundary functions.