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

Physical Review Research

Many-body spectral statistics of relativistic quantum billiard systems

Xianzhang Chen, Zhen-Qi Chen, Liang Huang, Celso Grebogi, Ying‐Cheng Lai

Published Jan 26, 2023• 0 citations
Publisher

Abstract

In the field of quantum chaos, spectral statistics is one of the most extensively investigated characteristics. Despite a large body of existing literature, the effects of many-body interactions on the spectral statistics of relativistic quantum systems remain poorly understood. Treating electron-electron interactions through the one-orbital mean-field Hubbard model, we address this fundamental issue using graphene billiards with the geometric shape of a circular sector as prototypical systems. Our approach is to consider the two characteristically different cases where the statistics are Poisson and Gaussian orthogonal ensemble (GOE) so the corresponding classical dynamics are typically integrable and chaotic, respectively, and to systematically investigate how the statistics change as the Hubbard interaction strength increases from zero. We find that, for energies near the Dirac point, the Hubbard interactions have a significant effect on the spectral statistics. Regardless of the type of spectral statistics to begin with, increasing the Hubbard interaction strength up to a critical value causes the statistics to approach GOE, rendering more applicable the random matrix theory. As the interaction strength increases beyond the critical value, the statistics evolve toward Poisson, due to the emergence of an energy gap rendering the quasiparticles massive. We also find that the energy levels above and below the Dirac point can exhibit different statistics, and the many-body interactions have little effect on the statistics for levels far from the Dirac point. These results reveal the intriguing interplay between many-body interactions and the spectral statistics, which we develop a physical picture to understand.

DOI
10.1103/physrevresearch.5.013050

Related papers

Level statistics of the one-dimensional ionic Hubbard model

Jeannette De Marco, Luisa Tolle, Catalin-Mihai Halati +3 more

In this paper we analyze the spectral level statistics of the one-dimensional ionic Hubbard model, the Hubbard model with an alternating on-site potential. In particular, we focus on the statistics of the gap ratios between consecutive energy levels. This quantity is often used in order to signal whether a many-body system is integrable or chaotic. A chaotic system has typically the statistics of a Gaussian ensemble of random matrices while the spectral properties of the integrable system follow a Poisson statistics. We find that whereas the Hubbard model without alternating potential is known to be integrable and its spectral properties follow a Poissonian statistics, the presence of an alternating potential causes a drastic change in the spectral properties, which resemble the one of a Gaussian ensemble of random matrices. However, to uncover this behavior one has to separately consider the blocks of all symmetries of the ionic Hubbard model.

arXiv:2203.05379•Physical Review Research•Aug 11, 2022•13 citations

Energy-level statistics of interacting trapped bosons

Barnali Chakrabarti, Anindya Biswas, V. K. B. Kota +2 more

It is a well-established fact that statistical properties of energy-level spectra are the most efficient tool to characterize nonintegrable quantum systems. The statistical behavior of different systems such as complex atoms, atomic nuclei, two-dimensional Hamiltonians, quantum billiards, and noninteracting many bosons has been studied. The study of statistical properties and spectral fluctuations in interacting many-boson systems has developed interest in this direction. We are especially interested in weakly interacting trapped bosons in the context of Bose-Einstein condensation (BEC) as the energy spectrum shows a transition from a collective nature to a single-particle nature with an increase in the number of levels. However this has received less attention as it is believed that the system may exhibit Poisson-like fluctuations due to the existence of an external harmonic trap. Here we compute numerically the energy levels of the zero-temperature many-boson systems which are weakly interacting through the van der Waals potential and are confined in the three-dimensional harmonic potential. We study the nearest-neighbor spacing distribution and the spectral rigidity by unfolding the spectrum. It is found that an increase in the number of energy levels for repulsive BEC induces a transition from a Wigner-like form displaying level repulsion to the Poisson distribution for $P(s)$. It does not follow the Gaussian orthogonal ensemble prediction. For repulsive interaction, the lower levels are correlated and manifest level-repulsion. For intermediate levels $P(s)$ shows mixed statistics, which clearly signifies the existence of two energy scales: external trap and interatomic interaction, whereas for very high levels the trapping potential dominates, generating a Poisson distribution. Comparison with mean-field results for lower levels are also presented. For attractive BEC near the critical point we observe the Shnirelman-like peak near $s=0$, which signifies the presence of a large number of quasidegenerate states.

arXiv:1101.5469•Physical Review A•Jul 23, 2012•10 citations

Ensemble-averaged mean-field many-body level density: An indicator of integrable versus chaotic single-particle dynamics

Georg Maier, Carolyn Echter, Juan Diego Urbina +2 more

According to the quantum chaos paradigm, the nature of a system's classical dynamics, whether integrable or chaotic, is universally reflected in the fluctuations of its quantum spectrum. However, since many-body spectra in the mean-field limit are composed of independent single-particle energy levels, their spectral fluctuations always display Poissonian behavior and hence cannot be used to distinguish underlying chaotic from integrable single-particle dynamics. We demonstrate that this distinction can, instead, be revealed from the mean many-body level density (at fixed energy) and its variance after averaging over ensembles representing different types of single-particle dynamics. This is in strong contrast to the energy-averaged mean level density (of a given system) that is assumed not to carry such information and is routinely removed to focus on universal signatures. To support our claim we systematically analyze the role of single-particle level correlations, that enter through Poisson and random matrix statistics (of various symmetry classes) into the ensemble-averaged density of states and its variance, contrasting bosonic and fermionic many-body systems. Our analytical study, together with extensive numerical simulations for systems with N≥5 particles consistently reveal significant differences (up to an order of magnitude for fermions and even larger for bosons) in the mean many-body level densities, depending on the nature of the underlying dynamics. Notably, in the fermionic case Poisson-type single-particle level fluctuations precisely cancel contributions from indistinguishability, such that the average many-body spectral density equals the (Thomas-Fermi) volume term. We further highlight the difference between the mean level density and its variance as functions of the total energy E and the excitation energy Q.

arXiv:2501.02918•Physical review. E•May 5, 2025•2 citations

Transitions in eigenvalue and wavefunction structure in (1+2) -body random matrix ensembles with spin

Manan Vyas, V. K. B. Kota, N. D. Chavda

Finite interacting Fermi systems with a mean-field and a chaos generating two-body interaction are modeled by one plus two-body embedded Gaussian orthogonal ensemble of random matrices with spin degree of freedom [called EGOE(1+2)-s]. Numerical calculations are used to demonstrate that, as lambda , the strength of the interaction (measured in the units of the average spacing of the single-particle levels defining the mean-field), increases, generically there is Poisson to GOE transition in level fluctuations, Breit-Wigner to Gaussian transition in strength functions (also called local density of states) and also a duality region where information entropy will be the same in both the mean-field and interaction defined basis. Spin dependence of the transition points lambda_{c} , lambdaF, and lambdad , respectively, is described using the propagator for the spectral variances and the formula for the propagator is derived. We further establish that the duality region corresponds to a region of thermalization. For this purpose we compared the single-particle entropy defined by the occupancies of the single-particle orbitals with thermodynamic entropy and information entropy for various lambda values and they are very close to each other at lambda=lambdad.

Physical Review E•Mar 12, 2010•25 citations

Manifestations of chaos in relativistic quantum systems - A study based on out-of-time-order correlator

Chen-Di Han, Hong-Ya Xu, Liang Huang +1 more

Previous work in the field of relativistic quantum chaos has revealed initial evidence that the manifestations of classical chaos in relativistic quantum systems tend to be weakened as compared with those in nonrelativistic quantum systems. To place this finding on a firmer ground, we investigate the relativistic quantum fingerprints of classical chaos using the out-of-time-order correlator (OTOC). OTOC has recently been applied to a number of fields in physics and it holds special promises for the field of quantum chaos, but existing work focused exclusively on nonrelativistic quantum systems. Calculating and analyzing OTOC for relativistic quantum billiard systems described by the massless Dirac equation, we find that the signatures of classical chaos are indeed characteristically less pronounced in relativistic than in nonrelativistic quantum systems. The finding is substantiated by studying four different aspects of OTOC. Firstly, in the energy eigenspace, in the short time regime there is a complete lack of any signature of chaos associated with the evolution of OTOC in the Dirac billiard due to the relativistic quantum phenomenon of Zitterbewegung, in contrast to nonrelativistic quantum billiard systems where chaos leaves behind a quite distinct signature. Secondly, weakening of the relativistic quantum manifestations of classical chaos occurs in the long time regime as well. Thirdly, evolution of the wave packet based OTOC also reveals that the fingerprints of chaos are much less pronounced in the Dirac billiard systems as compared with the Schrödinger billiards. Fourthly, the level spacing statistics of the OTOC operators in systems with classically integrable and chaotic dynamics, which are characteristically distinct in the Schrödinger billiards, bear a strong similarity in the Dirac billiards, indicating again the dwindling effects of chaos. The OTOC based results suggest that the impacts of classical chaos are generally weaker in relativistic than in nonrelativistic quantum systems, a fundamental issue that warrants further investigation.

Physics Open•Sep 4, 2019•4 citations

Circular and Elliptical Neutrino Billiards: A Semiclassical Approach

Barbara Dietz

Quantum signatures of classical chaos in the spectral properties of non-relativistic quantum billiards are well understood by now based on the semiclassical and the random matrix theory aproach.The experimental detection of relativistic phenomena featured by graphene implicated an increasing interest in the properties of graphene billiards and relativistic quantum billiards consisting of massless spin-1/2 particles governed by the Dirac equation and confined to a planar domain.Confinement is attained by imposing appropriate boundary conditions on the wave function components.We chose those proposed for neutrino billiards.By now numerical studies have been performed on the spectral properties of neutrino billiards of various shapes.Our main focus is the question to what extent the semiclassical approach in terms of trace formulae applies to such systems.We derive trace formulae for circular and elliptical neutrino billiards.The associated eigenvalue equations have been solved analytically for the circular billiard.We present analytical solutions for elliptical neutrino billiards.

Acta Physica Polonica A•Nov 1, 2019•7 citations
In the field of quantum chaos, spectral statistics is one of the most extensively investigated characteristics. Despite a large body of existing literature, the effects of many-body interactions on the spectral statistics of relativistic quantum systems remain poorly understood. Treating electron-electron interactions through the one-orbital mean-field Hubbard model, we address this fundamental issue using graphene billiards with the geometric shape of a circular sector as prototypical systems. Our approach is to consider the two characteristically different cases where the statistics are Poisson and Gaussian orthogonal ensemble (GOE) so the corresponding classical dynamics are typically integrable and chaotic, respectively, and to systematically investigate how the statistics change as the Hubbard interaction strength increases from zero. We find that, for energies near the Dirac point, the Hubbard interactions have a significant effect on the spectral statistics. Regardless of the type of spectral statistics to begin with, increasing the Hubbard interaction strength up to a critical value causes the statistics to approach GOE, rendering more applicable the random matrix theory. As the interaction strength increases beyond the critical value, the statistics evolve toward Poisson, due to the emergence of an energy gap rendering the quasiparticles massive. We also find that the energy levels above and below the Dirac point can exhibit different statistics, and the many-body interactions have little effect on the statistics for levels far from the Dirac point. These results reveal the intriguing interplay between many-body interactions and the spectral statistics, which we develop a physical picture to understand.
In this paper we analyze the spectral level statistics of the one-dimensional ionic Hubbard model, the Hubbard model with an alternating on-site potential. In particular, we focus on the statistics of the gap ratios between consecutive energy levels. This quantity is often used in order to signal whether a many-body system is integrable or chaotic. A chaotic system has typically the statistics of a Gaussian ensemble of random matrices while the spectral properties of the integrable system follow a Poisson statistics. We find that whereas the Hubbard model without alternating potential is known to be integrable and its spectral properties follow a Poissonian statistics, the presence of an alternating potential causes a drastic change in the spectral properties, which resemble the one of a Gaussian ensemble of random matrices. However, to uncover this behavior one has to separately consider the blocks of all symmetries of the ionic Hubbard model.
It is a well-established fact that statistical properties of energy-level spectra are the most efficient tool to characterize nonintegrable quantum systems. The statistical behavior of different systems such as complex atoms, atomic nuclei, two-dimensional Hamiltonians, quantum billiards, and noninteracting many bosons has been studied. The study of statistical properties and spectral fluctuations in interacting many-boson systems has developed interest in this direction. We are especially interested in weakly interacting trapped bosons in the context of Bose-Einstein condensation (BEC) as the energy spectrum shows a transition from a collective nature to a single-particle nature with an increase in the number of levels. However this has received less attention as it is believed that the system may exhibit Poisson-like fluctuations due to the existence of an external harmonic trap. Here we compute numerically the energy levels of the zero-temperature many-boson systems which are weakly interacting through the van der Waals potential and are confined in the three-dimensional harmonic potential. We study the nearest-neighbor spacing distribution and the spectral rigidity by unfolding the spectrum. It is found that an increase in the number of energy levels for repulsive BEC induces a transition from a Wigner-like form displaying level repulsion to the Poisson distribution for P(s)P(s)P(s). It does not follow the Gaussian orthogonal ensemble prediction. For repulsive interaction, the lower levels are correlated and manifest level-repulsion. For intermediate levels P(s)P(s)P(s) shows mixed statistics, which clearly signifies the existence of two energy scales: external trap and interatomic interaction, whereas for very high levels the trapping potential dominates, generating a Poisson distribution. Comparison with mean-field results for lower levels are also presented. For attractive BEC near the critical point we observe the Shnirelman-like peak near s=0s=0s=0, which signifies the presence of a large number of quasidegenerate states.
According to the quantum chaos paradigm, the nature of a system's classical dynamics, whether integrable or chaotic, is universally reflected in the fluctuations of its quantum spectrum. However, since many-body spectra in the mean-field limit are composed of independent single-particle energy levels, their spectral fluctuations always display Poissonian behavior and hence cannot be used to distinguish underlying chaotic from integrable single-particle dynamics. We demonstrate that this distinction can, instead, be revealed from the mean many-body level density (at fixed energy) and its variance after averaging over ensembles representing different types of single-particle dynamics. This is in strong contrast to the energy-averaged mean level density (of a given system) that is assumed not to carry such information and is routinely removed to focus on universal signatures. To support our claim we systematically analyze the role of single-particle level correlations, that enter through Poisson and random matrix statistics (of various symmetry classes) into the ensemble-averaged density of states and its variance, contrasting bosonic and fermionic many-body systems. Our analytical study, together with extensive numerical simulations for systems with N≥5 particles consistently reveal significant differences (up to an order of magnitude for fermions and even larger for bosons) in the mean many-body level densities, depending on the nature of the underlying dynamics. Notably, in the fermionic case Poisson-type single-particle level fluctuations precisely cancel contributions from indistinguishability, such that the average many-body spectral density equals the (Thomas-Fermi) volume term. We further highlight the difference between the mean level density and its variance as functions of the total energy E and the excitation energy Q.
Finite interacting Fermi systems with a mean-field and a chaos generating two-body interaction are modeled by one plus two-body embedded Gaussian orthogonal ensemble of random matrices with spin degree of freedom [called EGOE(1+2)-s]. Numerical calculations are used to demonstrate that, as lambda , the strength of the interaction (measured in the units of the average spacing of the single-particle levels defining the mean-field), increases, generically there is Poisson to GOE transition in level fluctuations, Breit-Wigner to Gaussian transition in strength functions (also called local density of states) and also a duality region where information entropy will be the same in both the mean-field and interaction defined basis. Spin dependence of the transition points lambda_{c} , lambdaF, and lambdad , respectively, is described using the propagator for the spectral variances and the formula for the propagator is derived. We further establish that the duality region corresponds to a region of thermalization. For this purpose we compared the single-particle entropy defined by the occupancies of the single-particle orbitals with thermodynamic entropy and information entropy for various lambda values and they are very close to each other at lambda=lambdad.
Previous work in the field of relativistic quantum chaos has revealed initial evidence that the manifestations of classical chaos in relativistic quantum systems tend to be weakened as compared with those in nonrelativistic quantum systems. To place this finding on a firmer ground, we investigate the relativistic quantum fingerprints of classical chaos using the out-of-time-order correlator (OTOC). OTOC has recently been applied to a number of fields in physics and it holds special promises for the field of quantum chaos, but existing work focused exclusively on nonrelativistic quantum systems. Calculating and analyzing OTOC for relativistic quantum billiard systems described by the massless Dirac equation, we find that the signatures of classical chaos are indeed characteristically less pronounced in relativistic than in nonrelativistic quantum systems. The finding is substantiated by studying four different aspects of OTOC. Firstly, in the energy eigenspace, in the short time regime there is a complete lack of any signature of chaos associated with the evolution of OTOC in the Dirac billiard due to the relativistic quantum phenomenon of Zitterbewegung, in contrast to nonrelativistic quantum billiard systems where chaos leaves behind a quite distinct signature. Secondly, weakening of the relativistic quantum manifestations of classical chaos occurs in the long time regime as well. Thirdly, evolution of the wave packet based OTOC also reveals that the fingerprints of chaos are much less pronounced in the Dirac billiard systems as compared with the Schrödinger billiards. Fourthly, the level spacing statistics of the OTOC operators in systems with classically integrable and chaotic dynamics, which are characteristically distinct in the Schrödinger billiards, bear a strong similarity in the Dirac billiards, indicating again the dwindling effects of chaos. The OTOC based results suggest that the impacts of classical chaos are generally weaker in relativistic than in nonrelativistic quantum systems, a fundamental issue that warrants further investigation.
Quantum signatures of classical chaos in the spectral properties of non-relativistic quantum billiards are well understood by now based on the semiclassical and the random matrix theory aproach.The experimental detection of relativistic phenomena featured by graphene implicated an increasing interest in the properties of graphene billiards and relativistic quantum billiards consisting of massless spin-1/2 particles governed by the Dirac equation and confined to a planar domain.Confinement is attained by imposing appropriate boundary conditions on the wave function components.We chose those proposed for neutrino billiards.By now numerical studies have been performed on the spectral properties of neutrino billiards of various shapes.Our main focus is the question to what extent the semiclassical approach in terms of trace formulae applies to such systems.We derive trace formulae for circular and elliptical neutrino billiards.The associated eigenvalue equations have been solved analytically for the circular billiard.We present analytical solutions for elliptical neutrino billiards.