⚛️ quantum physics

Analytic level statistics of the quantum Rabi model

This paper derives an analytic expression for the level spacing and ratio statistics of the quantum Rabi model by analyzing its high-energy asymptotic spectrum as a uniformly spaced ladder modulated by decaying oscillations, a theoretical prediction that is confirmed by exact diagonalization data.

Original authors: Zi-Min Li, Long-Jie Li, German Sierra, Murray T. Batchelor

Published 2026-11-02

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Zi-Min Li, Long-Jie Li, German Sierra, Murray T. Batchelor

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the quantum world, where particles behave more like waves than solid objects, the energy levels of a system are not random. They are like the rungs on a ladder, each one representing a specific state the system can occupy. For decades, physicists have studied how these rungs are spaced to understand the nature of the system itself. If the rungs are scattered without any pattern, the system is considered simple and predictable. If they push each other away in a very specific way, the system is chaotic and complex. Between these two extremes lies a vast landscape of quantum systems that do not fit neatly into either category. Understanding the precise spacing of these energy rungs is crucial because it reveals the hidden rules governing how light and matter interact at the smallest scales.

One such system is the quantum Rabi model, a fundamental description of a simple two-part system where a tiny particle with two possible states interacts with a single wave of light. This model is a cornerstone for understanding how atoms absorb and emit light, yet its energy levels have a peculiar structure that has long puzzled researchers. While the overall pattern of energy levels looks like a regular ladder, a closer look reveals a subtle, rhythmic wobble in the spacing between the rungs. For a long time, scientists could only describe this wobble through complex computer simulations, lacking a clear mathematical formula to explain exactly how the spacing was distributed.

A team of researchers has now bridged this gap by deriving a precise analytic expression for the distribution of these spacings. They focused on a specific high-energy range where the system behaves in a predictable way, allowing them to count the rungs and measure the gaps between them with mathematical certainty. Their work reveals that the wobble in the spacing is not random noise but a structured oscillation that splits the gaps into two distinct groups. By counting how often these gaps exceed a certain size, the researchers were able to map out the entire distribution of spacings. They found that the pattern forms a distinctive shape with two sharp peaks, resembling twin mountains, which are determined by the strength of the interaction between the particle and the light wave.

The researchers confirmed their mathematical predictions by comparing them against massive computer simulations of the model. The match was exact: the computer-generated data, which involved calculating millions of energy levels, aligned perfectly with the new formula. This agreement held true for the overall shape of the distribution, the range of possible gap sizes, and the specific positions of the twin peaks. The study also extended this analysis to the ratio of consecutive gaps, a measurement that does not require the complex adjustments usually needed to compare different energy scales. This ratio statistic also matched the computer simulations closely, showing a strong concentration near a specific value that shifts as the system moves to higher energies.

What makes this discovery significant is that it provides a clear, analytical window into a system that was previously understood only through numerical approximation. The researchers showed that as the energy levels get higher, the oscillation in the spacing becomes smaller, and the two groups of gaps move closer together, eventually merging into a single, uniform ladder. However, within any finite window of energy levels, the twin peaks remain visible, serving as a fingerprint of the underlying quantum interaction. The study demonstrates that even in a system that appears to be settling into a simple, regular pattern, there is a rich, decaying structure that can be captured and predicted with simple counting methods. This approach offers a new way to analyze not just the quantum Rabi model, but potentially other complex quantum systems where the energy levels follow a similar, slowly changing rhythm.

Technical Summary: Analytic Level Statistics of the Quantum Rabi Model

Problem Statement
Level statistics are fundamental for characterizing correlations in quantum spectra, typically distinguishing between Poisson statistics (integrable systems) and Wigner–Dyson statistics (random-matrix ensembles). While the Quantum Rabi Model (QRM)—a minimal model of light-matter interaction—has been studied numerically for its spectral correlations and as a diagnostic for quantum chaos, a rigorous analytic derivation of the spacing distribution for a finite high-energy window within a fixed parity sector has been lacking. Previous work established that the QRM spectrum approaches a delta measure at unit spacing in the asymptotic limit, but the detailed structure of the spacing density before this limit is reached, specifically the formation of twin peaks and the behavior of consecutive spacing ratios, remained unanalyzed.

Methodology
The authors derive an analytic expression for the level spacing density, $P(s)$ , and the ratio of consecutive spacings, $\rho_n$ , by leveraging the high-energy asymptotic spectrum of the QRM within a fixed parity sector.

Asymptotic Spectrum: The analysis utilizes the large- $n$ asymptotic expansion of the eigenvalues $\lambda_p(n)$ in a parity sector $p = \pm 1$ . The spectrum is described as a uniformly spaced ladder ( $n - g^2$ ) modulated by a decaying oscillation with amplitude proportional to $n^{-1/4}$ and an alternating sign determined by the parity and the index $n$ .
Branch Separation: The alternating sign of the oscillation separates consecutive spacings into two distinct branches, $s_+(n)$ and $s_-(n)$ , which are reflections of each other about $s=1$ .
Counting Approach: Instead of standard spectral unfolding, the authors employ a "counting" method. They define a finite index window $[m_-, m_+)$ containing many oscillations. By counting the number of intervals where the spacing branches exceed a given threshold $s$ , they derive a cumulative count function $N(s)$ .
Continuum Approximation: The discrete count is approximated by an integral over the continuous variable $n$ , assuming the oscillation amplitude and phase vary slowly relative to the index step. This yields an analytic form for the spacing density $P(s)$ by differentiating the cumulative count.
Ratio Statistics: A similar leading-order analysis is applied to the ratio of consecutive spacings, $\rho_n = \min(s_{n-1}, s_n) / \max(s_{n-1}, s_n)$ , which inherently cancels local energy scale variations and requires no unfolding.

Key Contributions and Results

Analytic Spacing Density $P(s)$ : The paper derives an explicit formula for the spacing density (Eq. 23) valid for $s > 1$ , with the full distribution obtained via reflection symmetry $P(s) = P(2-s)$ .
- Support: The support of the distribution is determined by the oscillation amplitude at the lower end of the window ( $m_-$ ), setting the range $[1 - r m_-^{-1/4}, 1 + r m_-^{-1/4}]$ .
- Twin Peaks: The density exhibits two square-root cusps (twin peaks) located at $s = 1 \pm r m_+^{-1/4}$ , determined by the oscillation amplitude at the upper end of the window ( $m_+$ ).
- Validation: Histograms generated from exact diagonalization of the QRM Hamiltonian show excellent agreement with the analytic density, accurately reproducing the support boundaries and the positions of the twin peaks.
Ratio Statistics: The authors derive an analytic expression for the distribution of the ratio of consecutive spacings.
- Characteristic Ratio: A characteristic ratio $\rho^*$ is identified, determined by the oscillation amplitude at the upper window boundary ( $m_+$ ). The ratio density shows an enhancement near this value.
- Comparison: The resulting ratio statistics agree closely with exact diagonalization data and differ significantly from the smooth Poisson and Gaussian Orthogonal Ensemble (GOE) references, reflecting the correlated alternating modulation of the QRM spectrum.
High-Energy Behavior: As the index window shifts to higher energies (increasing $m_-$ and $m_+$ ) while maintaining the condition of sampling many oscillations, the decaying oscillation amplitude causes the twin peaks of $P(s)$ and the characteristic ratio $\rho^*$ to converge toward unity. This reflects the spectrum's approach to a uniformly spaced ladder.

Significance and Claims
The paper claims to provide the first analytic derivation of the spacing density for the QRM in a finite high-energy window within a fixed parity sector. The significance lies in:

Connecting Correlations to Origin: The analytic expression directly links the observed spectral correlations (twin peaks, specific support) to the physical origin of the high-energy oscillation (the decaying Laguerre-modulated tunneling).
Methodological Extension: The counting approach is presented as a tool applicable to other systems possessing controlled spectral asymptotics, slowly varying spacing branches, and windows containing many oscillations.
Refinement of Limits: The work bridges the gap between the known asymptotic limit (delta measure at unit spacing) and the finite-energy regime, characterizing the "fine structure" of the spectrum before the limiting law is reached.

The authors explicitly state that the analysis is a leading-order continuum approximation and that low-energy statistics and finer structures beyond this approximation remain topics for future study. They do not propose new experimental setups but suggest that the derived asymptotics for related models (asymmetric and two-photon QRM) offer concrete settings for applying this counting method.

More like this