Level Spacing Distribution Research Articles

Distributions of Consecutive Level Spacings of Gaussian Unitary Ensemble and Their Ratio: ab initio Derivation

Abstract In recent studies of many-body localization in nonintegrable quantum systems, the distribution of the ratio of two consecutive energy level spacings, $r_n=(E_{n+1}-E_n)/(E_{n}-E_{n-1})$ or $\tilde{r}_n=\min (r_n,r_n^{-1})$, has been used as a measure to quantify the chaoticity, alternative to the more conventional distribution of the level spacings, $s_n=\bar{\rho }(E_n)(E_{n+1}-E_n)$, as the former makes unnecessary the unfolding required for the latter. Based on our previous work on the Tracy–Widom approach to the Jánossy densities, we present analytic expressions for the joint probability distribution of two consecutive eigenvalue spacings and the distribution of their ratio for the Gaussian unitary ensemble (GUE) of random Hermitian $N\times N$ matrices at $N\rightarrow \infty$, in terms of a system of differential equations. As a showcase of the efficacy of our results for characterizing an approach to quantum chaoticity, we contrast them to arguably the most ideal of all quantum-chaotic spectra: the zeroes of the Riemann $\zeta$ function on the critical line at increasing heights.

Anomalous quantum transport in fractal lattices

Fractal lattices are self-similar structures with repeated patterns on different scales. Quantum transport through such structures is subtle due to the possible co-existence of localized and extended states. Here, we study the dynamical properties of two fractal lattices, the Sierpiński gasket and the Sierpiński carpet. While the gasket exhibits sub-diffusive behavior, sub-ballistic transport occurs in the carpet. We show that the different dynamical behavior is in line with qualitative differences of the systems’ spectral properties. Specifically, in contrast to the Sierpiński carpet, the Sierpiński gasket exhibits an inverse power-law behavior of the level spacing distribution. As a possible technological application, we discuss a memory effect in the Sierpiński gasket which allows to read off the phase information of an initial state from the spatial distribution after long evolution times. We also show that interpolating between fractal and regular lattices allows for flexible tuning between different transport regimes.

Intermediate spectral statistics of rational triangular quantum billiards.

Triangular billiards whose angles are rational multiples of π are one of the simplest examples of pseudo-integrable models with intriguing classical and quantum properties. We perform an extensive numerical study of spectral statistics of eight quantized rational triangles, six belonging to the family of right-angled Veech triangles and two obtuse rational triangles. Large spectral samples of up to one million energy levels were calculated for each triangle, which permits one to determine their spectral statistics with great accuracy. It is demonstrated that they are of the intermediate type, sharing some features with chaotic systems, like level repulsion, and some with integrable systems, like exponential tails of the level spacing distributions. Another distinctive feature of intermediate spectral statistics is a finite value of the level compressibility. The short-range statistics such as the level spacing distributions, and long-range statistics such as the number variance and spectral form factors were analyzed in detail. An excellent agreement between the numerical data and the model of gamma distributions is revealed.

Two mass-imbalanced atoms in a hard-wall trap: Deep learning integrability of many-body systems.

The study of integrable systems has led to significant advancements in our understanding of many-body physics. We design a series of numerical experiments to analyze the integrability of a mass-imbalanced two-body system through energy-level statistics and deep learning of wave functions. The level spacing distributions are fitted by a Brody distribution and the fitting parameter ω is found to separate the integrable and nonintegrable mass ratios by a critical line ω=0. The convolutional neural network built from the probability density images could identify the transition points between integrable and nonintegrable systems with high accuracy, yet in a much shorter computation time. A brilliant example of the network's ability is to identify a new integrable mass ratio 1/3 by learning from the known integrable case of equal mass, with a remarkable network confidence of 98.22%. The robustness of our neural networks is further enhanced by adversarial learning, where samples are generated by standard and quantum perturbations mixed in the probability density images and the wave functions, respectively.

Spacing distribution for quantum Rabi models * *Dedicated to Mr. Hiromichi Shinohara, who was instrumental in establishing NTT-IFM, on his 70th birthday.

The asymmetric quantum Rabi model (AQRM) is a fundamental model in quantum optics describing the interaction of light and matter. Besides its immediate physical interest, the AQRM possesses an intriguing mathematical structure which is far from being completely understood. In this paper, we focus on the distribution of the level spacing, the difference between consecutive eigenvalues of the AQRM in the limit of high energies, i.e. large quantum numbers. In the symmetric case, that is the quantum Rabi model (QRM), the spacing distribution for each parity (given by the Z2 -symmetry) is fully clarified by an asymptotic expression derived by de Monvel and Zielinski, though some questions remain for the full spectrum spacing. However, in the general AQRM case, there is no parity decomposition for the eigenvalues. In connection with numerically exact studies and recent theoretical results, we describe the spacing distribution for the AQRM which is characterized by a new type of periodicity and symmetric behavior of the distribution with respect to the bias parameter, reflecting the hidden symmetry of the AQRM. In addition, we observe in the AQRM the excited state quantum phase transition for large values of the bias parameter, analogous to the QRM with large qubit energy, and an internal symmetry of the level spacing distribution for fixed bias. This novel symmetry is independent from the symmetry for half-integer bias and not explained by current theoretical knowledge.

Level-spacing distribution of localized phases induced by quasiperiodic potentials

Level-spacing distribution of localized phases induced by quasiperiodic potentials

Mixed (regular and chaotic) dynamics under the channeling of the high energy positrons in [100] direction of the silicon crystal and Podolskiy-Narimanov distribution

The character of motion (regular or chaotic) on the quantum level manifests itself in the statistics of the set of the system's energy levels. The completely regular case generates the sequence of the levels with exponential (Poisson) level spacing distribution; the completely chaotic one — with Wigner distribution. The most interesting case is the co-existence of the regular and chaotic motion domains in the phase space of the system under consideration. The assumption of independent generation of two level sequencies by one chaotic domain and by all regular domains leads to Berry-Robnik distribution. However, the presence of the chaotic motion domain in the phase space affects the levels produced by the regular domains via the so-called chaos-assisted tunneling (CAT) that leads to Podolskiy-Narimanov distribution of the level spacings. This distribution needs the mean amplitude of the tunnel transition and the relative contribution of the regular domains to the mean level density as the parameters. Using their estimations for the transverse motion of the high energy positrons channeling near [100] direction in the silicon crystal we found that Podolskiy-Narimanov distribution demonstrates the best agreement with the level spacing distribution.

Classification of same-gate quantum circuits and their space-time symmetries with application to the level-spacing distribution

We study Floquet systems with translationally invariant nearest-neighbor two-site gates. Depending on the order in which the gates are applied on an N-site system with periodic boundary conditions, there are factorially many different circuit configurations. We prove that there are only N−1 different spectrally equivalent classes, which can be viewed either as a generalization of the brick wall or of the staircase configuration. Every class, characterized by two integers, has a nontrivial space-time symmetry with important implications for the level-spacing distribution—a standard indicator of quantum chaos. Namely, in order to study chaoticity one should not look at eigenphases of the Floquet propagator itself, but rather at the spectrum of an appropriate root of the propagator. Published by the American Physical Society 2024

Signatures of quantum chaos and fermionization in the incoherent transport of bosonic carriers in the Bose-Hubbard chain.

We analyze the stationary current of Bose particles across the Bose-Hubbard chain connected to a battery, focusing on the effect of interparticle interactions. It is shown that the current magnitude drastically decreases as the strength of interparticle interactions exceeds the critical value which marks the transition to quantum chaos in the Bose-Hubbard Hamiltonian. We found that this transition is well reflected in the nonequilibrium many-body density matrix of the system. Namely, the level-spacing distribution for eigenvalues of the density matrix changes from Poisson to Wigner-Dyson distributions. With the further increase of the interaction strength, the Wigner-Dyson spectrum statistics change back to the Poisson statistics which now marks fermionization of the Bose particles. With respect to the stationary current, this leads to the counter-intuitive dependence of the current magnitude on the particle number.

Quantum integrability and chaos in a periodic Toda lattice with balanced loss-gain.

We consider an equal-mass quantum Toda lattice with balanced loss-gain for two and three particles. The two-particle Toda lattice is integrable, and two integrals of motion that are in involution have been found. The bound-state energy and the corresponding eigenfunctions have been obtained numerically for a few low-lying states. The three-particle quantum Toda lattice with balanced loss-gain and velocity-mediated coupling admits mixed phases of integrability and chaos depending on the value of the loss-gain parameter. We have obtained analytic expressions for two integrals of motion that are in involution. Although an analytic expression for the third integral has not been found, the numerical investigation suggests integrability below a critical value of the loss-gain strength and chaos above this critical value. The level spacing distribution changes from the Wigner-Dyson to the Poisson distribution as the loss-gain parameter passes through this critical value and approaches zero. An identical behavior is seen in terms of the gap-ratio distribution of the energy levels. The existence of mixed phases of quantum integrability and chaos in the specified ranges of the loss-gain parameter has also been confirmed independently via the study of level repulsion and complexity in higher order excited states.

Quantum graphs and microwave networks as narrow-band filters for quantum and microwave devices.

We investigate properties of the transmission amplitude of quantum graphs and microwave networks composed of regular polygons such as triangles and squares. We show that for the graphs composed of regular polygons, with the edges of the length l, the transmission amplitude displays a band of transmission suppression with some narrow peaks of full transmission. The peaks are distributed symmetrically with respect to the symmetry axis kl=π, where k is the wave vector. For microwave networks the transmission peak amplitudes are reduced and their symmetry is broken due to the influence of internal absorption. We demonstrate that for the graphs composed of the same polygons but separated by the edges of length l^{'}<l, the transmission spectrum is generally not symmetric according to the axis kl^{'}=π. We also show that graphs composed of regular polygons of different size with the edges being irrational numbers are not fully chaotic and their level spacing distribution and the spectral rigidity are well described by the Berry-Robnik distributions. Moreover, the transmission spectrum of such a graph displays peaks which are very close to one. Furthermore, the microwave networks are investigated in the time-domain using short Gaussian pulses. In this case the delay-time distributions, though very sensitive to the internal structure of the networks, show the sequences of transmitted peaks with the amplitudes much smaller than the input one. The analyzed properties of the graphs and networks suggest that they can be effectively used to manipulate quantum and wave transport.

Quantum chaos in interacting Bose-Bose mixtures

The appearance of chaotic quantum dynamics significantly depends on the symmetry properties of a system, and in cold atomic systems many of these can be experimentally controlled. In this work, we systematically study the emergence of quantum chaos in a minimal system describing one-dimensional harmonically trapped Bose-Bose mixtures by tuning the particle-particle interactions. Using an improved exact diagonalization scheme, we examine the transition from integrability to chaos when the inter-component interaction changes from weak to strong. Our study is based on the analysis of the level spacing distribution and the distribution of the matrix elements of observables in terms of the eigenstate thermalization hypothesis and their dynamics. We show that one can obtain strong signatures of chaos by increasing the inter-component interaction strength and breaking the symmetry of intra-component interactions.

Each state in a one-dimensional disordered system has two localization lengths when the Hilbert space is constrained

In disordered systems, the amplitudes of the localized states will decrease exponentially away from their centers and the localization lengths characterize such decrease. In this paper, we find a model in which each eigenstate is decreasing at two distinct rates. The model is a one-dimensional disordered system with a constrained Hilbert space: all eigenstates should be orthogonal to a state , where is a given exponentially localized state. Although the dimension of the Hilbert space is only reduced by 1, the amplitude of each state will decrease at one rate near its center and at another rate in the rest of the region. Depending on , it is also possible that all states are changed from localized states to extended states. In such a case, the level spacing distribution is different from that of the three well-known ensembles of the random matrices. This indicates that a new ensemble of random matrices exists in this model. Finally we discuss the physics behind such phenomena and propose an experiment to observe them.

Many-Body Quantum Chaos and Emergence of Ginibre Ensemble

We show that non-Hermitian Ginibre random matrix behaviors emerge in spatially extended many-body quantum chaotic systems in the space direction, just as Hermitian random matrix behaviors emerge in chaotic systems in the time direction. Starting with translational invariant models, which can be associated with dual transfer matrices with complex-valued spectra, we show that the linear ramp of the spectral form factor necessitates that the dual spectra have nontrivial correlations, which in fact fall under the universality class of the Ginibre ensemble, demonstrated by computing the level spacing distribution and the dissipative spectral form factor. As a result of this connection, the exact spectral form factor for the Ginibre ensemble can be used to universally describe the spectral form factor for translational invariant many-body quantum chaotic systems in the scaling limit where t and L are large, while the ratio between L and L_{Th}, the many-body Thouless length is fixed. With appropriate variations of Ginibre models, we analytically demonstrate that our claim generalizes to models without translational invariance as well. The emergence of the Ginibre ensemble is a genuine consequence of the strongly interacting and spatially extended nature of the quantum chaotic systems we consider, unlike the traditional emergence of Hermitian random matrix ensembles.

Energy-level statistics in planar fractal tight-binding models

In this study, we examine the statistics of level spectra in a noninteracting electron gas confined to a Sierpi\ifmmode \acute{n}\else \'{n}\fi{}ski carpet lattice. These lattices are constructed using two types of the $\mathit{self}$ and $gene$ patterns, and they are categorized by the area-perimeter scaling law. The singularly continuous spectra, along with the nearest level spacing distribution and gap-ratio distribution, reveal a critical phase different from both extended and localized phases. This critical phase differs from the behavior near the Anderson model's metal-insulator transition. The Wigner-like conjecture is confirmed for both lattice classes, indicating Gaussian orthogonal symmetry. A similar observation was made in a quasiperiodic lattice [Phys. Rev. Lett. 80, 3996 (1998)]. The self-similar nature of fractals also leads to level clustering behavior.

Quantum chaos and circuit parameter optimization

We consider quantum chaos diagnostics of the variational circuit states at random parameters and explore their connection to the circuit expressibility and optimizability. By measuring the operator spreading coefficient and the eigenvalue spectrum of the modular Hamiltonian of the reduced density matrix, we identify the emergence of universal random matrix ensembles in high-depth circuit states. The diagnostics that use the eigenvalue spectrum, e.g. operator spreading and entanglement entropy, turn out to be more accurate measures of the variational quantum algorithm optimization efficiency than those that use the level spacing distribution of the entanglement spectrum, such as r-statistics or spectral form factors.

Implications of Spectral Interlacing for Quantum Graphs

Quantum graphs are ideally suited to studying the spectral statistics of chaotic systems. Depending on the boundary conditions at the vertices, there are Neumann and Dirichlet graphs. The latter ones correspond to totally disassembled graphs with a spectrum being the superposition of the spectra of the individual bonds. According to the interlacing theorem, Neumann and Dirichlet eigenvalues on average alternate as a function of the wave number, with the consequence that the Neumann spectral statistics deviate from random matrix predictions. There is, e.g., a strict upper bound for the spacing of neighboring Neumann eigenvalues given by the number of bonds (in units of the mean level spacing). Here, we present analytic expressions for level spacing distribution and number variance for ensemble averaged spectra of Dirichlet graphs in dependence of the bond number, and compare them with numerical results. For a number of small Neumann graphs, numerical results for the same quantities are shown, and their deviations from random matrix predictions are discussed.

Quantum Entanglement and Spectral Form Factor

We replace a Hamiltonian with a modular Hamiltonian in the spectral form factor and the level spacing distribution function. This study establishes a connection between quantities within Quantum Entanglement and Quantum Chaos. To have a universal study for Quantum Entanglement, we consider the Gaussian random 2-qubit model. The maximum violation of Bell’s inequality demonstrates a positive correlation with the entanglement entropy. Thus, the violation plays an equivalent role as Quantum Entanglement. We first provide an analytical estimation of the relation between quantum entanglement quantities and the dip when a subregion only has one qubit. The time of the first dip is monotone for entanglement entropy. The dynamics in a subregion is independent of the initial state at a late time. It is one of the signaling conditions for classical chaos. We also extend our analysis to the Gaussian random 3-qubit state, and it indicates a similar result. The simulation shows that the level spacing distribution function approaches GUE at a late time. In the end, we develop a technique within QFT to the spectral form factor for its relation to an n-sheet manifold. We apply the technology to a single interval in CFT2 and the spherical entangling surface in \(\mathcal {N}=4\) super Yang-Mills theory. The result is one for both cases, but the Rényi entropy can depend on the Rényi index. For the case of CFT2, it indicates the difference between the continuum and discrete spectrum.

Patterned random matrices: deviations from universality

We investigate the level spacing distribution for three ensembles of real symmetric matrices having additional structural constraint to reduce the number of independent entries to only in contrast to the for a real symmetric matrix of size n × n. We derive all the results analytically exactly for the matrices and show that spacing distribution display a range of behaviour based on the structural constraint. The spacing distribution of the ensemble of reverse circulant matrices with additional zeros is found to fall slower than exponential for larger spacing while that of symmetric circulant matrices has poisson spacings. The palindromic symmetric Toeplitz matrices on the other hand show level repulsion but the distribution is significantly different from Wigner. The behaviour of spacings for all the three ensembles clearly show the departure from universal result of Wigner distribution for real symmetric matrices. The deviation from universality continues in large n cases as well, which we study numerically.

Wave localization in number-theoretic landscapes

We investigate the localization of waves in aperiodic structures that manifest the characteristic multiscale complexity of certain arithmetic functions with a central role in number theory. In particular, we study the eigenspectra and wave localization properties of tight-binding Schr\"{o}dinger equation models with on-site potentials distributed according to the Liouville function $\lambda(n)$, the M\"{o}bius function $\mu(n)$, and the Legendre sequence of quadratic residues modulo a prime (QRs). We employ Multifractal Detrended Fluctuation Analysis (MDFA) and establish the multifractal scaling properties of the energy spectra in these systems. Moreover, by systematically analyzing the spatial eigenmodes and their level spacing distributions, we show the absence of level repulsion with broadband localization across the entire energy spectra. Our study introduces deterministic aperiodic systems whose eigenmodes are all strongly localized in realistic finite one-dimensional systems and provides opportunities for novel quantum and classical devices of particular importance to cold-atom experiments in engineered speckle potentials and enhanced light-matter interactions.