Dynamical billiards

A dynamical billiard is a dynamical system in which a particle alternates between free motion (typically as a straight line) and specular reflections from a boundary. When the particle hits the boundary it reflects from it without loss of speed (i.e. elastic collisions). Billiards are Hamiltonian idealizations of the game of billiards, but where the region contained by the boundary can have shapes other than rectangular and even be multidimensional. Dynamical billiards may also be studied on non-Euclidean geometries; indeed, the first studies of billiards established their ergodic motion on surfaces of constant negative curvature. The study of billiards which are kept out of a region, rather than being kept in a region, is known as outer billiard theory.

A particle moving inside the Bunimovich stadium, a well-known chaotic billiard. See the Software section for making such an animation.

The motion of the particle in the billiard is a straight line, with constant energy, between reflections with the boundary (a geodesic if the Riemannian metric of the billiard table is not flat). All reflections are specular: the angle of incidence just before the collision is equal to the angle of reflection just after the collision. The sequence of reflections is described by the billiard map that completely characterizes the motion of the particle.

Billiards capture all the complexity of Hamiltonian systems, from integrability to chaotic motion, without the difficulties of integrating the equations of motion to determine its Poincaré map. Birkhoff showed that a billiard system with an elliptic table is integrable.

Equations of motionEdit

The Hamiltonian for a particle of mass m moving freely without friction on a surface is:


where   is a potential designed to be zero inside the region   in which the particle can move, and infinity otherwise:


This form of the potential guarantees a specular reflection on the boundary. The kinetic term guarantees that the particle moves in a straight line, without any change in energy. If the particle is to move on a non-Euclidean manifold, then the Hamiltonian is replaced by:


where   is the metric tensor at point  . Because of the very simple structure of this Hamiltonian, the equations of motion for the particle, the Hamilton–Jacobi equations, are nothing other than the geodesic equations on the manifold: the particle moves along geodesics.

Notable billiards and billiard classesEdit

Hadamard's billiardsEdit

Hadamard's billiards concern the motion of a free point particle on a surface of constant negative curvature, in particular, the simplest compact Riemann surface with negative curvature, a surface of genus 2 (a two-holed donut). The model is exactly solvable, and is given by the geodesic flow on the surface. It is the earliest example of deterministic chaos ever studied, having been introduced by Jacques Hadamard in 1898.

Artin's billiardEdit

Artin's billiard considers the free motion of a point particle on a surface of constant negative curvature, in particular, the simplest non-compact Riemann surface, a surface with one cusp. It is notable for being exactly solvable, and yet not only ergodic but also strongly mixing. It is an example of an Anosov system. This system was first studied by Emil Artin in 1924.

Dispersing and semi-dispersing billiardsEdit

Let M be complete smooth Riemannian manifold without boundary, maximal sectional curvature of which is not greater than K and with the injectivity radius  . Consider a collection of n geodesically convex subsets (walls)  ,  , such that their boundaries are smooth submanifolds of codimension one. Let  , where   denotes the interior of the set  . The set   will be called the billiard table. Consider now a particle that moves inside the set B with unit speed along a geodesic until it reaches one of the sets Bi (such an event is called a collision) where it reflects according to the law “the angle of incidence is equal to the angle of reflection” (if it reaches one of the sets   ,  , the trajectory is not defined after that moment). Such dynamical system is called semi-dispersing billiard. If the walls are strictly convex, then the billiard is called dispersing. The naming is motivated by observation that a locally parallel beam of trajectories disperse after a collision with strictly convex part of a wall, but remain locally parallel after a collision with a flat section of a wall.

Dispersing boundary plays the same role for billiards as negative curvature does for geodesic flows causing the exponential instability of the dynamics. It is precisely this dispersing mechanism that gives dispersing billiards their strongest chaotic properties, as it was established by Yakov G. Sinai.[1] Namely, the billiards are ergodic, mixing, Bernoulli, having a positive Kolmogorov-Sinai entropy and an exponential decay of correlations.

Chaotic properties of general semi-dispersing billiards are not understood that well, however, those of one important type of semi-dispersing billiards, hard ball gas were studied in some details since 1975 (see next section).

General results of Dmitri Burago and Serge Ferleger[2] on the uniform estimation on the number of collisions in non-degenerate semi-dispersing billiards allow to establish finiteness of its topological entropy and no more than exponential growth of periodic trajectories.[3] In contrast, degenerate semi-dispersing billiards may have infinite topological entropy.[4]

Lorentz gas, aka Sinai billiardEdit

A particle moving inside the Sinai billiard, also known as Lorentz gas.

The table of the Lorentz gas (also known as Sinai billiard) is a square with a disk removed from its center; the table is flat, having no curvature. The billiard arises from studying the behavior of two interacting disks bouncing inside a square, reflecting off the boundaries of the square and off each other. By eliminating the center of mass as a configuration variable, the dynamics of two interacting disks reduces to the dynamics in the Sinai billiard.

The billiard was introduced by Yakov G. Sinai as an example of an interacting Hamiltonian system that displays physical thermodynamic properties: almost all (up to a measure zero) of its possible trajectories are ergodic and it has a positive Lyapunov exponent.

Sinai's great achievement with this model was to show that the classical Boltzmann–Gibbs ensemble for an ideal gas is essentially the maximally chaotic Hadamard billiards.

Bunimovich stadiumEdit

The table called the Bunimovich stadium is a rectangle capped by semicircles, a shape called a stadium. Until it was introduced by Leonid Bunimovich, billiards with positive Lyapunov exponents were thought to need convex scatters, such as the disk in the Sinai billiard, to produce the exponential divergence of orbits. Bunimovich showed that by considering the orbits beyond the focusing point of a concave region it was possible to obtain exponential divergence.

Magnetic billiardsEdit

Movement of a charged particle inside a Sinai billiard with perpendicular magnetic field.

Magnetic billiards represent billiards where a charged particle is propagating under the presence of a perpendicular magnetic field. As a result, the particle trajectory changes from a straight line into an arc of a circle. The radius of this circle is inversely proportional to the magnetic field strength. Such billiards have been useful in real world applications of billiards, typically modelling nanodevices (see Applications).

Generalized billiardsEdit

Generalized billiards (GB) describe a motion of a mass point (a particle) inside a closed domain   with the piece-wise smooth boundary  . On the boundary   the velocity of point is transformed as the particle underwent the action of generalized billiard law. GB were introduced by Lev D. Pustyl'nikov in the general case,[5] and, in the case when   is a parallelepiped[6] in connection with the justification of the second law of thermodynamics. From the physical point of view, GB describe a gas consisting of finitely many particles moving in a vessel, while the walls of the vessel heat up or cool down. The essence of the generalization is the following. As the particle hits the boundary  , its velocity transforms with the help of a given function  , defined on the direct product   (where   is the real line,   is a point of the boundary and   is time), according to the following law. Suppose that the trajectory of the particle, which moves with the velocity  , intersects   at the point   at time  . Then at time   the particle acquires the velocity  , as if it underwent an elastic push from the infinitely-heavy plane  , which is tangent to   at the point  , and at time   moves along the normal to   at   with the velocity  . We emphasize that the position of the boundary itself is fixed, while its action upon the particle is defined through the function  .

We take the positive direction of motion of the plane   to be towards the interior of  . Thus if the derivative  , then the particle accelerates after the impact.

If the velocity  , acquired by the particle as the result of the above reflection law, is directed to the interior of the domain  , then the particle will leave the boundary and continue moving in   until the next collision with  . If the velocity   is directed towards the outside of  , then the particle remains on   at the point   until at some time   the interaction with the boundary will force the particle to leave it.

If the function   does not depend on time  ; i.e.,  , the generalized billiard coincides with the classical one.

This generalized reflection law is very natural. First, it reflects an obvious fact that the walls of the vessel with gas are motionless. Second the action of the wall on the particle is still the classical elastic push. In the essence, we consider infinitesimally moving boundaries with given velocities.

It is considered the reflection from the boundary   both in the framework of classical mechanics (Newtonian case) and the theory of relativity (relativistic case).

Main results: in the Newtonian case the energy of particle is bounded, the Gibbs entropy is a constant,[6][7][8] (in Notes) and in relativistic case the energy of particle, the Gibbs entropy, the entropy with respect to the phase volume grow to infinity,[6][8] (in Notes), references to generalized billiards.

Quantum chaosEdit

The quantum version of the billiards is readily studied in several ways. The classical Hamiltonian for the billiards, given above, is replaced by the stationary-state Schrödinger equation   or, more precisely,


where   is the Laplacian. The potential that is infinite outside the region   but zero inside it translates to the Dirichlet boundary conditions:


As usual, the wavefunctions are taken to be orthonormal:


Curiously, the free-field Schrödinger equation is the same as the Helmholtz equation,




This implies that two and three-dimensional quantum billiards can be modelled by the classical resonance modes of a radar cavity of a given shape, thus opening a door to experimental verification. (The study of radar cavity modes must be limited to the transverse magnetic (TM) modes, as these are the ones obeying the Dirichlet boundary conditions).

The semi-classical limit corresponds to   which can be seen to be equivalent to  , the mass increasing so that it behaves classically.

As a general statement, one may say that whenever the classical equations of motion are integrable (e.g. rectangular or circular billiard tables), then the quantum-mechanical version of the billiards is completely solvable. When the classical system is chaotic, then the quantum system is generally not exactly solvable, and presents numerous difficulties in its quantization and evaluation. The general study of chaotic quantum systems is known as quantum chaos.

A particularly striking example of scarring on an elliptical table is given by the observation of the so-called quantum mirage.


Billiards, both quantum and classical, have been applied in several areas of physics to model quite diverse real world systems. Examples include ray-optics,[9] lasers,[10][11] acoustics,[12] optical fibers (e.g. double-clad fibers [13][14]), or quantum-classical correspondence.[15] One of their most frequent application is to model particles moving inside nanodevices, for example quantum dots,[16][17] pn-junctions,[18] antidot superlattices,[19][20] among others. The reason for this broadly spread effectiveness of billiards as physical models resides on the fact that in situations with small amount of disorder or noise, the movement of e.g. particles like electrons, or light rays, is very much similar to the movement of the point-particles in billiards. In addition, the energy conserving nature of the particle collisions is a direct reflection of the energy conservation of Hamiltonian mechanics.


Open source software to simulate billiards exist for various programming languages. From most recent to oldest, existing software are: DynamicalBilliards.jl (Julia), Bill2D (C++) and Billiard Simulator (Matlab). The animations present on this page were done with DynamicalBilliards.jl.

See alsoEdit


  1. ^ "Archived copy" (PDF). Archived from the original (PDF) on 2013-12-31. Retrieved 2014-06-06.{{cite web}}: CS1 maint: archived copy as title (link)
  2. ^ Burago, D.; Ferleger, S.; Kononenko, A. (1 January 1998). "Uniform Estimates on the Number of Collisions in Semi-Dispersing Billiards". Annals of Mathematics. 147 (3): 695–708. doi:10.2307/120962. JSTOR 120962.
  3. ^ Burago, D.; Ferleger, S. (26 May 1997). "Topological Entropy Of Semi-Dispersing Billiards". Ergodic Theory and Dynamical Systems. 18 (4): 791. doi:10.1017/S0143385798108246. S2CID 122549772.
  4. ^ Burago, D. (1 February 2006). "Semi-dispersing billiards of infinite topological entropy". Ergodic Theory and Dynamical Systems. 26 (1): 45–52. doi:10.1017/S0143385704001002. S2CID 121644309.
  5. ^ Pustyl'nikov, L. D. (1999). "The law of entropy increase and generalized billiards". Russian Mathematical Surveys. 54 (3): 650–651. Bibcode:1999RuMaS..54..650P. doi:10.1070/rm1999v054n03abeh000168. S2CID 250902640.
  6. ^ a b c Pustyl'nikov, L. D. (1995). "Poincaré models, rigorous justification of the second law of thermodynamics from mechanics, and the Fermi acceleration mechanism". Russian Mathematical Surveys. 50 (1): 145–189. Bibcode:1995RuMaS..50..145P. doi:10.1070/rm1995v050n01abeh001663. S2CID 250875392.
  7. ^ Pustyl'nikov, L. D. (2005). "Generalized Newtonian periodic billiards in a ball". Russian Mathematical Surveys. 60 (2): 365–366. Bibcode:2005RuMaS..60..365P. doi:10.1070/RM2005v060n02ABEH000839. S2CID 250856558.
  8. ^ a b Deryabin, Mikhail V.; Pustyl'nikov, Lev D. (2007). "Nonequilibrium Gas and Generalized Billiards". Journal of Statistical Physics. 126 (1): 117–132. Bibcode:2007JSP...126..117D. doi:10.1007/s10955-006-9250-4. S2CID 55957240.
  9. ^ Kouznetsov, Dmitrii; Moloney, Jerome V. (September 2004). "Boundary behaviour of modes of a Dirichlet Laplacian". Journal of Modern Optics. 51 (13): 1955–1962. Bibcode:2004JMOp...51.1955K. doi:10.1080/09500340408232504. ISSN 0950-0340. S2CID 30880255.
  10. ^ Stone, A. Douglas (June 2010). "Chaotic billiard lasers". Nature. 465 (7299): 696–697. doi:10.1038/465696a. ISSN 1476-4687. PMID 20535191.
  11. ^ Gmachl, C. (1998-06-05). "High-Power Directional Emission from Microlasers with Chaotic Resonators". Science. 280 (5369): 1556–1564. arXiv:cond-mat/9806183. Bibcode:1998Sci...280.1556G. doi:10.1126/science.280.5369.1556. PMID 9616111. S2CID 502055.
  12. ^ Koyanagi, Sin’ichiro; Nakano, Takeru; Kawabe, Tetsuji (2008-08-01). "Application of Hamiltonian of ray motion to room acoustics". The Journal of the Acoustical Society of America. 124 (2): 719–722. Bibcode:2008ASAJ..124..719K. doi:10.1121/1.2946714. ISSN 0001-4966. PMID 18681564.
  13. ^ Leproux, P.; S. Fevrier; V. Doya; P. Roy; D. Pagnoux (2003). "Modeling and optimization of double-clad fiber amplifiers using chaotic propagation of pump". Optical Fiber Technology. 7 (4): 324–339. Bibcode:2001OptFT...7..324L. doi:10.1006/ofte.2001.0361.
  14. ^ a b B. D. Lubachevsky and F. H. Stillinger, Geometric properties of random disk packings, J. Statistical Physics 60 (1990), 561-583 http://www.princeton.edu/~fhs/geodisk/geodisk.pdf
  15. ^ Stöckmann, H.-J.; Stein, J. (1990-05-07). "Quantum chaos in billiards studied by microwave absorption". Physical Review Letters. 64 (19): 2215–2218. Bibcode:1990PhRvL..64.2215S. doi:10.1103/PhysRevLett.64.2215. ISSN 0031-9007. PMID 10041617.
  16. ^ Ponomarenko, L. A.; Schedin, F.; Katsnelson, M. I.; Yang, R.; Hill, E. W.; Novoselov, K. S.; Geim, A. K. (2008-04-18). "Chaotic Dirac Billiard in Graphene Quantum Dots". Science. 320 (5874): 356–358. arXiv:0801.0160. Bibcode:2008Sci...320..356P. doi:10.1126/science.1154663. ISSN 0036-8075. PMID 18420930. S2CID 206511356.
  17. ^ Bird, Jonathan P., ed. (2003). Electron Transport in Quantum Dots. doi:10.1007/978-1-4615-0437-5. ISBN 978-1-4020-7459-2.
  18. ^ Chen, Shaowen; Han, Zheng; Elahi, Mirza M.; Habib, K. M. Masum; Wang, Lei; Wen, Bo; Gao, Yuanda; Taniguchi, Takashi; Watanabe, Kenji; Hone, James; Ghosh, Avik W. (2016-09-30). "Electron optics with p-n junctions in ballistic graphene". Science. 353 (6307): 1522–1525. arXiv:1602.08182. Bibcode:2016Sci...353.1522C. doi:10.1126/science.aaf5481. ISSN 0036-8075. PMID 27708099. S2CID 118443999.
  19. ^ Weiss, D.; Roukes, M. L.; Menschig, A.; Grambow, P.; von Klitzing, K.; Weimann, G. (1991-05-27). "Electron pinball and commensurate orbits in a periodic array of scatterers" (PDF). Physical Review Letters. 66 (21): 2790–2793. Bibcode:1991PhRvL..66.2790W. doi:10.1103/PhysRevLett.66.2790. ISSN 0031-9007. PMID 10043617.
  20. ^ Datseris, George; Geisel, Theo; Fleischmann, Ragnar (2019-04-30). "Robustness of ballistic transport in antidot superlattices". New Journal of Physics. 21 (4): 043051. Bibcode:2019NJPh...21d3051D. doi:10.1088/1367-2630/ab19cc. ISSN 1367-2630.


Sinai's billiardsEdit

  • Sinai, Ya. G. (1963). "[On the foundations of the ergodic hypothesis for a dynamical system of statistical mechanics]". Doklady Akademii Nauk SSSR (in Russian). 153 (6): 1261–1264. (in English, Sov. Math Dokl. 4 (1963) pp. 1818–1822).
  • Ya. G. Sinai, "Dynamical Systems with Elastic Reflections", Russian Mathematical Surveys, 25, (1970) pp. 137–191.
  • V. I. Arnold and A. Avez, Théorie ergodique des systèms dynamiques, (1967), Gauthier-Villars, Paris. (English edition: Benjamin-Cummings, Reading, Mass. 1968). (Provides discussion and references for Sinai's billiards.)
  • D. Heitmann, J.P. Kotthaus, "The Spectroscopy of Quantum Dot Arrays", Physics Today (1993) pp. 56–63. (Provides a review of experimental tests of quantum versions of Sinai's billiards realized as nano-scale (mesoscopic) structures on silicon wafers.)
  • S. Sridhar and W. T. Lu, "Sinai Billiards, Ruelle Zeta-functions and Ruelle Resonances: Microwave Experiments", (2002) Journal of Statistical Physics, Vol. 108 Nos. 5/6, pp. 755–766.
  • Linas Vepstas, Sinai's Billiards, (2001). (Provides ray-traced images of Sinai's billiards in three-dimensional space. These images provide a graphic, intuitive demonstration of the strong ergodicity of the system.)
  • N. Chernov and R. Markarian, "Chaotic Billiards", 2006, Mathematical survey and monographs nº 127, AMS.

Strange billiardsEdit

  • T. Schürmann and I. Hoffmann, The entropy of strange billiards inside n-simplexes. J. Phys. A28, page 5033ff, 1995. PDF-Document

Bunimovich stadiumEdit

Generalized billiardsEdit

  • M. V. Deryabin and L. D. Pustyl'nikov, "Generalized relativistic billiards", Reg. and Chaotic Dyn. 8(3), pp. 283–296 (2003).
  • M. V. Deryabin and L. D. Pustyl'nikov, "On Generalized Relativistic Billiards in External Force Fields", Letters in Mathematical Physics, 63(3), pp. 195–207 (2003).
  • M. V. Deryabin and L. D. Pustyl'nikov, "Exponential attractors in generalized relativistic billiards", Comm. Math. Phys. 248(3), pp. 527–552 (2004).

External linksEdit