Among the fundamental symmetry classes of quantum chaotic systems in Dyson's threefold way, the symplectic class is rarely observed in nature. Characterized by the strongest possible level repulsion in the energy spectrum, the symplectic symmetry class also implies a double (Kramers) degeneracy of levels. Studying the spectral statistics of three quantum particles (identical bosons or mass-imbalanced fermions) in a harmonic trap, we find numerical evidence for strong level repulsion in the regime of weak contact interactions. While the statistical indicators are consistent with quantum chaos in systems with symplectic symmetry, the absence of Kramers degeneracy rules out this symmetry. In the strongly-interacting unitary limit either Poissonian or stick statistics are observed (depending on commensurability of the mass ratio) indicating regular dynamics.
Among the fundamental symmetry classes of quantum chaotic systems in Dyson's threefold way, the symplectic class is rarely observed in nature. Characterized by the strongest possible level repulsion in the energy spectrum, the symplectic symmetry class also implies a double (Kramers) degeneracy of levels. Studying the spectral statistics of three quantum particles (identical bosons or mass-imbalanced fermions) in a harmonic trap, we find numerical evidence for strong level repulsion in the regime of weak contact interactions. While the statistical indicators are consistent with quantum chaos in systems with symplectic symmetry, the absence of Kramers degeneracy rules out this symmetry. In the strongly-interacting unitary limit either Poissonian or stick statistics are observed (depending on commensurability of the mass ratio) indicating regular dynamics.
We consider identical quantum bosons with weak contact interactions in a two-dimensional isotropic harmonic trap. When the interactions are turned off, the energy levels are equidistant and highly degenerate. At linear order in the coupling parameter, these degenerate levels split, and we study the patterns of this splitting. It turns out that the problem is mathematically identical to diagonalizing the quantum resonant system of the two-dimensional Gross-Pitaevskii equation, whose classical counterpart has been studied previously in the mathematical literature on turbulence. Our purpose is to explore the implications of the symmetries and energy bounds of this resonant system, previously studied for the classical case, for the quantum level splitting. Simplifications in computing the splitting spectrum numerically result from exploiting the symmetries. The highest energy state emanating from each unperturbed level is explicitly described by our analytics. We furthermore discuss the energy-level spacing distributions in the spirit of quantum chaos theory. After separating the eigenvalues into blocks with respect to the known conservation laws, we observe the Wigner-Dyson statistics within specific large blocks, which leaves little room for further integrable structures in the problem beyond the symmetries that are already explicitly known.
This is the first in a pair of articles that classify the configuration space and kinematic symmetry groups for N identical particles in one-dimensional traps experiencing Galilean-invariant two-body interactions. These symmetries explain degeneracies in the few-body spectrum and demonstrate how tuning the trap shape and the particle interactions can manipulate these degeneracies. The additional symmetries that emerge in the non-interacting limit and in the unitary limit of an infinitely strong contact interaction are sufficient to algebraically solve for the spectrum and degeneracy in terms of the one-particle observables. Symmetry also determines the degree to which the algebraic expressions for energy level shifts by weak interactions or nearly-unitary interactions are universal, i.e. independent of trap shape and details of the interaction. Identical fermions and bosons with and without spin are considered. This article sequentially analyzes the symmetries of one, two and three particles in asymmetric, symmetric, and harmonic traps; the sequel article treats the N particle case.
We investigate the quantum-classical correspondence in open quantum many-body systems using the SU(3) Bose-Hubbard trimer as a minimal model. Combining exact diagonalization with semiclassical Langevin dynamics, we establish a direct connection between classical trajectories characterized by fixed-point attractors, limit cycles, or chaos and the spectral and structural properties of the quantum steady state. We show that classical dynamical behavior, as quantified by the sign of the Lyapunov exponent, governs the level statistics of the steady-state density matrix: non-positive exponents associated with regular dynamics yield Poissonian statistics, while positive exponents arising from chaotic dynamics lead to Wigner-Dyson statistics. Strong symmetries constrain the system to lower-dimensional manifolds, suppressing chaos and enforcing localization, while weak symmetries preserve the global structure of the phase space and allow chaotic behavior to persist. To characterize phase-space localization, we introduce the phase-space inverse participation ratio IPR, which defines an effective dimension D of the Husimi distribution's support. We find that the entropy scales as S∝lnND, consistently capturing the classical nature of the underlying dynamics. This semiclassical framework, based on stochastic mixtures of coherent states, successfully reproduces not only observable averages but also finer features such as spectral correlations and localization properties. Our results demonstrate that dissipative quantum chaos is imprinted in the steady-state density matrix, much like in closed systems, and that the interplay between dynamical regimes and symmetry constraints can be systematically probed using spectral and phase-space diagnostics. These tools offer a robust foundation for studying ergodicity, localization, and non-equilibrium phases of open quantum systems.
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.
We present an exact diagonalization study of the spectral properties of bosons harmonically confined in a quasi-2D plane and interacting via repulsive Gaussian potential. We consider the lowest 100 energy levels for systems of N=12,16, and 20 bosons in two distinct regimes: (a) when the interaction energy is small compared to the trap energy (moderate interaction) and (b) when the interaction energy is comparable to the trap energy (strong interaction), for the nonrotating (L_{z}=0) as well as the rotating single-vortex state (L_{z}=N). For higher angular momenta, L_{z}=2N and L_{z}=3N, only the strong interaction regime is considered. While the nearest-neighbor spacing distribution (NNSD) P(s) and the ratios of consecutive level spacings distribution P(r) are used to study the short-range correlations, the Dyson-Mehta Δ_{3} statistic and the level number variance Σ^{2}(L) are used to examine the long-range correlations. In the moderate interaction regime, the nonrotating system exhibits Poisson distribution, a characteristic of the regular energy spectra. In the strong interaction regime, the nonrotating system exhibits chaotic behavior signified by GOE distribution. Furthermore, in the rotating case for the single-vortex state (L_{z}=N) in the moderate interaction regime, the system exhibits signatures of weak chaos with some degree of regularity in the energy-level spectra. However, in the strong interaction regime for the rotating case with L_{z}=N, 2N and 3N, the system exhibits strong chaotic behavior. The rotation is found to contribute to enhancement of chaotic behavior in the system for both the moderate and the strong interaction regimes. Our results of NNSD analysis are supported by the analysis of the ratios of consecutive level spacings distribution P(r), which does not involve unfolding.
Ergodicity is a fundamental principle of statistical mechanics underlying the behavior of generic quantum many-body systems.However, how this universal many-body quantum chaotic regime emerges due to interactions remains largely a puzzle.This paper demonstrates using both heuristic arguments and a microscopic calculation that a dephasing mechanism, similar to Altshuler-Aronov-Khmelnitskii dephasing in the theory of localization, underlies this transition to chaos.We focus on the behavior of the spectral form factor (SFF) as a function of "time" t, which characterizes level correlations in the many-body spectrum.The SFF can be expressed as a sum over periodic classical orbits and its behavior hinges on the interference of trajectories related to each other by a time translation.In the absence of interactions, time-translation symmetry is present for each individual particle, which leads to a fast exponential growth of the SFF and correspondingly loss of correlations between many-body levels.Interactions lead to dephasing, which disrupts interference, and breaks the massive time-translation symmetry down to a global time-translation/energy conservation.This in turn gives rise to the hallmark linear-in-t ramp in the SFF reflecting Wigner-Dyson level repulsion.This general picture is supported by a microscopic analysis of an interacting many-body model.Specifically, we study the complex SYK 2 + SYK 2 2 model, which allows to tune between an integrable and chaotic regime.It is shown that the dephasing mass vanishes in the former case, which maps to the noninteracting complex SYK 2 model via a time reparameterization.In contrast, the chaotic regime gives rise to dephasing, which suppresses the exponential ramp of the noninteracting theory and induces correlations between many-body levels.