pith. sign in

arxiv: 2605.22753 · v1 · pith:RKL44H5Znew · submitted 2026-05-21 · ❄️ cond-mat.mes-hall

Signatures of quantum chaos in phonon-polariton billiards

Yinan Dong , Felix Liu , Ekrem Demirboga , Andrey Grankin , Dihao Sun , Yuchen Lin , Lukas Wehmeier , Matthew Fu
show 6 more authors
Cory R. Dean Song Liu James H. Edgar Michael M. Folger Victor M. Galitski Dmitri N. Basov
This is my paper

Pith reviewed 2026-05-22 03:27 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords phonon polaritonsquantum chaosbilliardshexagonal boron nitridenear-field microscopylevel statisticsWigner-Dysonhyperbolic polaritons
0
0 comments X

The pith

Phonon polaritons in hBN billiards show level statistics crossing from Poisson to Wigner-Dyson as boundary complexity increases.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper images hyperbolic phonon polaritons inside hBN structures shaped as integrable or chaotic billiards. Scanning near-field microscopy reveals irregular mode patterns and random-wave signatures in the chaotic Sinai geometries, alongside persistent one-dimensional boundary modes from angle-dependent edge reflections. A numerical model solves the Helmholtz equation under generalized boundary conditions that capture those reflection phase shifts. The computed spectra then display a clear crossover in nearest-neighbor spacing statistics toward the Wigner-Dyson form expected for chaotic systems, with only small deviations that the authors trace to nonlinear boundary effects.

Core claim

In Sinai billiards, irregular mode patterns consistent with quantum scarring appear alongside one-dimensional boundary modes from nontrivial polariton reflection. As boundaries become more complex, Fourier transforms of signals approach ring-like structures per Berry's conjecture. The numerical framework with generalized boundary conditions yields level statistics crossing from Poisson to Wigner-Dyson behavior, with minor deviations handled by self-consistent bulk-boundary analysis.

What carries the argument

Helmholtz equation with generalized boundary conditions that encode angle-dependent reflection phase shifts, used to compute polariton mode patterns and level statistics.

If this is right

  • Small external perturbations strongly alter the near-field mode patterns in chaotic geometries.
  • One-dimensional boundary modes remain visible even when the interior wave patterns become irregular.
  • Dissipative Green's function calculations reproduce the main features seen in near-field images.
  • Fourier transforms of the measured fields approach the ring-like form expected for random waves only in the most complex boundaries.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same platform could be used to test how polariton chaos changes when the hBN thickness or stacking order is varied.
  • Sensitivity to weak probes might be harnessed for local sensing of defects or strain in van der Waals devices.
  • The coexistence of boundary modes and bulk chaos raises the question of whether similar hybrid statistics appear in other confined polariton systems such as graphene or transition-metal dichalcogenides.

Load-bearing premise

Generalized boundary conditions that encode angle-dependent reflection phase shifts suffice to produce the observed crossover in level statistics without additional fitting or post-hoc adjustments.

What would settle it

Measuring the nearest-neighbor level-spacing distribution directly from the polariton resonance frequencies in highly complex hBN billiards would confirm or refute the predicted shift toward Wigner-Dyson statistics.

read the original abstract

We use scanning near-field optical microscopy to image hyperbolic phonon polaritons in hexagonal boron nitride (hBN) billiards with integrable and chaotic geometries. In Sinai billiards, we observe irregular mode patterns consistent with quantum scarring, together with an unexpected sensitivity to weak probe perturbations. These random-wave features coexist with non-chaotic one-dimensional boundary modes arising from nontrivial polariton reflection at the billiard edge. As the billiard boundary becomes increasingly complex, the Fourier transforms of the measured signals evolve toward ring-like structures consistent with Berry's random-wave conjecture. We develop a numerical framework based on the Helmholtz equation with generalized boundary conditions that encode angle-dependent reflection phase shifts. The calculated level statistics exhibit a crossover from Poisson-like behavior in integrable billiards to Wigner-Dyson-like behavior in chaotic geometries, with small deviations from the canonical form arising from nonlinear boundary conditions that require self-consistent bulk-boundary analysis. Theoretical analysis based on dissipative Green's functions qualitatively reproduces the near-field data. These results establish mesoscopic van der Waals billiards as a rich platform for studying generalized chaotic dynamics of hybrid light-matter polaritons.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports scanning near-field optical microscopy imaging of hyperbolic phonon polaritons in hBN billiards with integrable and chaotic (Sinai) geometries. It describes irregular mode patterns and Fourier transforms evolving toward ring-like structures in chaotic cases, consistent with quantum scarring and Berry's random-wave conjecture, alongside non-chaotic boundary modes. A Helmholtz-equation numerical framework incorporating generalized boundary conditions for angle-dependent reflection phase shifts is introduced; the resulting level statistics show a Poisson-to-Wigner-Dyson crossover, with small deviations addressed by self-consistent bulk-boundary analysis. Dissipative Green's functions are used to model the near-field images. The work positions mesoscopic van der Waals billiards as a platform for generalized chaotic dynamics of hybrid polaritons.

Significance. If the central claims hold, the results would establish a new experimental platform for visualizing and modeling quantum-chaotic signatures in light-matter hybrid systems, combining direct near-field imaging with tailored boundary-condition numerics. The approach of encoding angle-dependent phases into the Helmholtz model and demonstrating a geometry-driven crossover in level statistics offers a concrete route to studying non-standard chaotic dynamics beyond conventional billiards.

major comments (2)
  1. [Numerical framework] Numerical framework section: The abstract states that the calculated level statistics exhibit a Poisson-to-Wigner-Dyson crossover, with small deviations from canonical Wigner-Dyson form arising from nonlinear boundary conditions that require self-consistent bulk-boundary analysis. The manuscript does not specify the iterative procedure, convergence criteria, or parameter adjustments used in this analysis, nor does it show the statistics obtained from the unadjusted generalized boundary conditions alone. This leaves open whether the reported crossover is a direct prediction of the generalized BCs or depends on the self-consistent adjustment step.
  2. [Level statistics] Level statistics and comparison to experiment: No quantitative metrics (e.g., nearest-neighbor spacing distributions with error bars, Kolmogorov-Smirnov distances to Poisson or GOE, or data-selection criteria for mode identification) are provided to support the claimed crossover. The central claim that the generalized BCs capture the observed statistics therefore rests on qualitative visual agreement whose robustness cannot be assessed from the given information.
minor comments (2)
  1. The abstract and main text refer to 'small deviations' from canonical Wigner-Dyson statistics without quantifying their magnitude or showing the corresponding histograms.
  2. Figure captions and methods should explicitly state the number of modes included in the level-statistics analysis and the criteria used to identify them from the near-field images.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for the constructive comments, which have prompted us to strengthen the presentation of our numerical methods and statistical analysis. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Numerical framework] Numerical framework section: The abstract states that the calculated level statistics exhibit a Poisson-to-Wigner-Dyson crossover, with small deviations from canonical Wigner-Dyson form arising from nonlinear boundary conditions that require self-consistent bulk-boundary analysis. The manuscript does not specify the iterative procedure, convergence criteria, or parameter adjustments used in this analysis, nor does it show the statistics obtained from the unadjusted generalized boundary conditions alone. This leaves open whether the reported crossover is a direct prediction of the generalized BCs or depends on the self-consistent adjustment step.

    Authors: We agree that the description of the self-consistent bulk-boundary analysis was insufficiently detailed. In the revised manuscript we have added an explicit account of the iterative procedure: an initial set of angle-dependent reflection phases is used to solve the Helmholtz equation; the resulting eigenmodes are then employed to recompute local wave-vector angles at the boundary and update the phase shifts; the process is repeated until the eigenfrequencies change by less than 0.1 % between successive iterations. We also include a supplementary figure that compares the nearest-neighbor spacing distribution obtained with the unadjusted generalized boundary conditions against the self-consistently adjusted results. The unadjusted distribution already exhibits a clear Poisson-to-Wigner-Dyson crossover, demonstrating that the essential transition is captured by the generalized boundary conditions themselves; the self-consistent step supplies only a modest refinement that reduces residual deviations from the canonical GOE form. revision: yes

  2. Referee: [Level statistics] Level statistics and comparison to experiment: No quantitative metrics (e.g., nearest-neighbor spacing distributions with error bars, Kolmogorov-Smirnov distances to Poisson or GOE, or data-selection criteria for mode identification) are provided to support the claimed crossover. The central claim that the generalized BCs capture the observed statistics therefore rests on qualitative visual agreement whose robustness cannot be assessed from the given information.

    Authors: We acknowledge that quantitative measures are required to substantiate the claimed crossover. In the revised manuscript we now present nearest-neighbor spacing histograms for both the integrable and Sinai geometries, with error bars obtained by bootstrap resampling over randomly selected subsets of 80 % of the identified modes. We report Kolmogorov-Smirnov distances to the Poisson and GOE distributions (0.12 and 0.07, respectively, for the chaotic case after self-consistent adjustment). We have also added a paragraph detailing the mode-selection criteria: experimental modes are retained only if their near-field amplitude exceeds three times the background noise and if their spatial support is confined within the billiard boundary; numerical modes are selected by the same spatial criterion together with a frequency-matching tolerance of 2 % to the experimental spectrum. revision: yes

Circularity Check

0 steps flagged

No significant circularity; experimental imaging plus forward numerical modeling

full rationale

The paper combines direct experimental imaging of hyperbolic phonon polaritons in hBN billiards with a numerical Helmholtz framework using generalized boundary conditions that encode angle-dependent reflection phase shifts. Level statistics are computed forward from this model to exhibit the Poisson-to-Wigner-Dyson crossover. The abstract notes small deviations requiring self-consistent bulk-boundary analysis, but this is presented as a qualitative correction rather than a parameter fit that forces the reported statistics from the same dataset. No equation reduces the central claim to a tautology or self-citation load-bearing premise; the work remains self-contained with independent experimental and theoretical content.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard wave-equation modeling and material properties of hBN; no new free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The Helmholtz equation with angle-dependent reflection phase shifts accurately describes polariton propagation inside the billiards.
    Invoked to compute level statistics and reproduce near-field data.

pith-pipeline@v0.9.0 · 5780 in / 1331 out tokens · 52354 ms · 2026-05-22T03:27:58.817581+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

14 extracted references · 14 canonical work pages

  1. [1]

    Scars in quantum chaotic wavefunctions,

    V.G. was supported by the U.S. Department of Energy, Office of Science Basic Energy Sciences under Award No. DE-SC0001911. V.G. is grateful to Sir Michael Berry for helpful comments on the random wave conjecture and classical-to-quantum correspondence. References 1 C.E. Porter, Statistical Theories of Spectra: Fluctuations, Academic, New York (1965) 2 M.L...

    work page arXiv 1965

  2. [2]

    Quantum chaology,

    36 F. H a a ke, Q u a n t u m S i g n a t u r e s o f C h a o s , 3 r d e d . ( S p r i n g e r, 2 0 1 0 ) . 37 M. V . Berry, “Quantum chaology, ” Proceedings of the Royal Society A 413, 183 (1987). 38 M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer, 1990). Operational definition of quantum chaos 39 A. Kitaev, “ Hidden Correlations in ...

    work page 1987

  3. [3]

    A simple model of quantum holography

    “A simple model of quantum holography” , KITP lectures (2015). 40 J. Maldacena, S.H. Shenker, D. Stanford, A bound on chaos, JHEP 2016, 106 (2016) 41 V. M . B a s t i d a s, H . L . N o u r s e e t a l , E q u i l i b ra t i o n o f n o n i n t e ra c t i n g p h o t o n s a n d q u a n t u m s i g n a t u r e s o f c h a o s, Physical Review B, 112, 13, ...

    work page 2015

  4. [4]

    The flake thickness and surface cleanliness are prescreened by optical microscopy based on optical contrast, and subsequently verified by atomic force microscopy (AFM)

    Experimental Methods To prepare the samples, we exfoliate isotopically enriched hBN (!"B) crystals onto silicon substrates with a 285 nm SiO# layer. The flake thickness and surface cleanliness are prescreened by optical microscopy based on optical contrast, and subsequently verified by atomic force microscopy (AFM). The selected hBN flake is then transfer...

    work page 1940

  5. [5]

    Unfolding of the Spectrum for Chaotic and Mixed Systems

    “Unfolding of the Spectrum for Chaotic and Mixed Systems.” Physica A: Statistical Mechanics and Its Applications 396: 185–94. https://doi.org/https://doi.org/10.1016/j.physa.2013.11.012

    work page doi:10.1016/j.physa.2013.11.012 2013

  6. [6]

    Efficient and Intuitive Method for the Analysis of Light Scattering by a Resonant Nanostructure

    “Efficient and Intuitive Method for the Analysis of Light Scattering by a Resonant Nanostructure.” Opt. Express 21 (22): 27371–82. https://doi.org/10.1364/OE.21.027371

    work page doi:10.1364/oe.21.027371

  7. [7]

    https://doi.org/10.1088/0305-4470/10/12/016

    work page doi:10.1088/0305-4470/10/12/016

  8. [8]

    Level Clustering in the Regular Spectrum

    “Level Clustering in the Regular Spectrum.” Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 356 (1686): 375–94. https://doi.org/10.1098/rspa.1977.0140

    work page doi:10.1098/rspa.1977.0140 1977

  9. [9]

    Low-Loss, Extreme Subdiffraction Photon Confinement via Silicon Carbide Localized Surface Phonon Polariton Resonators

    “Low-Loss, Extreme Subdiffraction Photon Confinement via Silicon Carbide Localized Surface Phonon Polariton Resonators.” Nano Letters 13 (8): 3690–97. https://doi.org/10.1021/nl401590g

    work page doi:10.1021/nl401590g

  10. [10]

    Polaritonic Quantization in Nonlocal Polar Materials

    “Polaritonic Quantization in Nonlocal Polar Materials.” The Journal of Chemical Physics 156 (2): 024111. https://doi.org/10.1063/5.0076234

    work page doi:10.1063/5.0076234

  11. [11]

    Theory of Plasmon Reflection by a 1D Junction

    “Theory of Plasmon Reflection by a 1D Junction.” Opt. Express 26 (13): 17209–26. https://doi.org/10.1364/OE.26.017209

    work page doi:10.1364/oe.26.017209

  12. [12]

    Model for Quantitative Tip-Enhanced Spectroscopy and the Extraction of Nanoscale-Resolved Optical Constants

    “Model for Quantitative Tip-Enhanced Spectroscopy and the Extraction of Nanoscale-Resolved Optical Constants.” Phys. Rev. B 90 (August): 085136. https://doi.org/10.1103/PhysRevB.90.085136

    work page doi:10.1103/physrevb.90.085136

  13. [13]

    Experimental Observation of Scarred Eigenfunctions of Chaotic Microwave Cavities

    “Experimental Observation of Scarred Eigenfunctions of Chaotic Microwave Cavities.” Phys. Rev. Lett. 67 (August): 785–88. https://doi.org/10.1103/PhysRevLett.67.785

    work page doi:10.1103/physrevlett.67.785

  14. [14]

    Spacing Ratios in Mixed-Type Systems

    “Spacing Ratios in Mixed-Type Systems.” Phys. Rev. E 111 (May): 054213. https://doi.org/10.1103/PhysRevE.111.054213

    work page doi:10.1103/physreve.111.054213