Classical and quantum chaos in bean- and peanut-shaped billiards

The boundary of a billiard system plays a crucial role in shaping its dynamics, which may be integrable, mixed, or fully chaotic. When a boundary has varying curvature, it offers a unique setting to study the relation between classical chaos and quantum behaviour. In this study, we introduce two geometrically distinct billiards: a bean- and a peanut-shaped billiard. These systems incorporate both focusing and defocusing walls with no neutral segments. Our study reveals a strong correlation between classical and quantum dynamics. Analysis of billiard flow diagrams confirms sensitivity to initial conditions-a defining feature of chaos. Poincaré maps further show the phase space intricately woven with regions of chaotic motion and stability islands. We employ both statistical and dynamical measures to characterise quantum chaos. Statistical indicator includes nearest-neighbour spacing distribution, level spacing ratios, and spectral staircase function, while dynamical measures includes out-of-time-order correlators and spectral complexity. We also observe eigenfunction scarring in both the billiards.

💡 Research Summary

In this work the authors introduce two novel two‑dimensional billiard tables whose boundaries are defined by quartic curves with varying curvature: a “bean‑shaped” billiard based on the Bean curve and a “peanut‑shaped” billiard based on a Cassini oval. By choosing specific parameter sets (a₁=2, b₁=6 for the bean and a₂=1, b₂=√1.375 for the peanut) the authors obtain domains that contain both concave (focusing) and convex (defocusing) wall segments, but no flat portions. The bean billiard possesses a single mirror symmetry (x = 0) while the peanut billiard has two orthogonal mirror symmetries (x = 0 and y = 0).

The classical dynamics are explored through continuous billiard flow and the associated discrete Poincaré map (boundary coordinate s and momentum p = sin θ). Trajectory simulations demonstrate extreme sensitivity to initial conditions, and Lyapunov‑type exponential divergence is inferred from the flow diagrams. Phase‑space portraits reveal a mixed phase space: a chaotic sea interwoven with stability islands. The single‑symmetry bean table shows larger chaotic regions aligned with the symmetry axis, whereas the double‑symmetry peanut table exhibits a richer island structure due to the additional mirror line.

Quantum mechanically the stationary Schrödinger equation with Dirichlet boundary conditions is solved numerically for each geometry. The eigenvalue spectra are analysed with several statistical tools. Nearest‑Neighbour Spacing Distributions (NNSD) and Level‑Spacing Ratios (LSR) both display intermediate statistics between Poisson (integrable) and Wigner‑Dyson (GOE) forms, reflecting the underlying mixed classical dynamics. The spectral staircase function follows the Weyl law on average but shows pronounced oscillations that correspond to periodic‑orbit contributions.

Beyond traditional spectral statistics, the authors employ two recently proposed dynamical diagnostics of quantum chaos. The Out‑of‑Time‑Order Correlator (OTOC) is computed for a set of local operators; it grows exponentially at early times, with a growth rate that is larger for the peanut billiard, mirroring the larger classical Lyapunov exponent. Spectral Complexity (SC), defined as the entropy of the Fourier‑transformed level density, is higher for the peanut geometry, indicating a more intricate spectral structure. Both measures correlate well with the classical indicators, providing a unified picture of chaos across the classical‑quantum divide.

Eigenfunctions are examined for scarred patterns. In the bean billiard, a dominant scar aligns with the single mirror axis, while in the peanut billiard two distinct scars appear along each symmetry line. These scars are associated with unstable periodic orbits identified in the classical flow, confirming the semiclassical link between classical instability and quantum wavefunction localization.

Overall, the paper demonstrates that (i) mixed‑curvature boundaries generate simultaneous classical and quantum chaos; (ii) the presence and number of mirror symmetries strongly influence the distribution of stability islands, the strength of eigenfunction scarring, and the values of OTOC and SC; and (iii) modern quantum‑chaos diagnostics (OTOC, SC) are consistent with traditional spectral statistics. The results suggest that such asymmetric, curvature‑varying billiards can serve as versatile testbeds for studies of quantum chaos, wave‑guide design, and the development of quantum simulators that exploit controlled chaotic dynamics.