We propose to quantify the complexity of nonequilibrium steady state density operators, as well as of long-lived Liouvillian decay modes, in terms of the level spacing distribution of their spectra. Based on extensive numerical studies in a variety of models, some solvable and some unsolved, we conjecture that the integrability of density operators (e.g., the existence of an algebraic procedure for their construction in finitely many steps) is signaled by a Poissonian level statistics, whereas in the generic nonintegrable cases one finds level statistics of a Gaussian unitary ensemble of random matrices. Eigenvalue statistics can therefore be used as an efficient tool to identify integrable quantum nonequilibrium systems.
We propose to quantify the complexity of nonequilibrium steady state density operators, as well as of long-lived Liouvillian decay modes, in terms of the level spacing distribution of their spectra. Based on extensive numerical studies in a variety of models, some solvable and some unsolved, we conjecture that the integrability of density operators (e.g., the existence of an algebraic procedure for their construction in finitely many steps) is signaled by a Poissonian level statistics, whereas in the generic nonintegrable cases one finds level statistics of a Gaussian unitary ensemble of random matrices. Eigenvalue statistics can therefore be used as an efficient tool to identify integrable quantum nonequilibrium systems.
The Lindblad description of an open quantum system gives rise to two types of integrability, since the nonequilibrium steady state can be integrable independently of the Liouvillian. Taking boundary-driven and dephasing spin chains as a representative example, we discriminate Liouvillian and steady-state chaos by combining level spacing statistics and an extension of the eigenstate thermalization hypothesis to open quantum systems. Moreover, we analyze the structure of the steady states by expanding it in the basis of Pauli strings and comparing the weight of strings of different lengths. We show that the natural expectation that integrable steady states are "simple" (i.e., built from few-body local operators) does not hold: the steady states of both chaotic and integrable models have relevant contributions coming from Pauli strings of all possible lengths, including long-range and many-body interactions. Nevertheless, we show that one can effectively use the operator-size distribution to distinguish chaotic and integrable steady states.
We study generic open quantum systems with Markovian dissipation, focusing on a class of stochastic Liouvillian operators of Lindblad form with independent random dissipation channels (jump operators) and a random Hamiltonian. We perform a thorough numerical study, focusing on global spectral features, the spectral gap, and the steady-state purity and statistics. We establish that all properties follow three different regimes as a function of the dissipation strength, whose boundaries depend on the particular observable. Within each regime, we determine the scaling exponents with the dissipation strength and system size. On the analytical side, we compute the average spectral gap at arbitrary dissipation and provide simple closed-form expressions for the asymptotic values at strong and weak dissipation strength. We also consider spectral correlations in generic complex spectra, such as that of the Liouvillian. We introduce and study, both analytically and numerically, the complex ratio of the nearest-neighbor spacing by next-to-nearest-neighbor spacing. Besides the usual level repulsion, we find that complex spacing ratios have a nontrivial angular dependence which offers a distinct signature of quantum chaos. It distinguishes integrable systems, which support the flat ratio distribution characteristic of Poisson statistics, and chaotic dissipative quantum systems, which conform to Random Matrix Theory statistics. We apply this new signature of quantum chaos to random Liouvillians and boundary-driven dissipative spin-chains. The results of this thesis can help understand the long-time dynamics, steady-state properties, and spectral correlations of generic dissipative systems.
While the concepts of quantum many-body integrability and chaos are of fundamental importance for the understanding of quantum matter, their precise definition has so far remained an open question. In this Letter, we introduce an alternative indicator for quantum many-body integrability and chaos, which is based on the statistics of eigenstates by means of nearest-neighbor subsystem trace distances. We show that this provides us with a faithful classification through extensive numerical simulations for a large variety of paradigmatic model systems including random matrix theories, free fermions, Bethe-ansatz solvable systems, and models of many-body localization. While existing indicators, such as those obtained from level-spacing statistics, have already been utilized with great success, they also face limitations. This concerns, for instance, the quantum many-body kicked top, which is exactly solvable but classified as chaotic in certain regimes based on the level-spacing statistics, while our introduced indicator signals the expected quantum many-body integrability. We discuss the universal behaviors we observe for the nearest-neighbor trace distances and point out that our indicator might be useful also in other contexts such as for the many-body localization transition.
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.
The statistics of gaps between quantum energy levels is a hallmark criterion in quantum chaos and quantum integrability studies. The relevant distributions corresponding to exactly integrable vs. fully chaotic systems are universal and described by the Poisson vs. Wigner-Dyson curves. In the transitional regime between integrability and chaos, the distributions are much less universal and have not been understood quantitatively until now. We point out that the relevant statistics that controls these distributions is that of the matrix elements of the nonintegrable perturbation Hamiltonian in the energy eigenbasis of the unperturbed integrable system. With this insight, we formulate a simple random matrix ensemble that correctly reproduces the level spacing distributions in a variety of test systems. For the distribution of matrix elements appearing in our construction, we furthermore discover surprising universal features: across a variety of physical systems with diverse degrees of freedom, these distributions are dominated by simple power laws.
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.