Eigenvalue Statistics as an Indicator of Integrability of Nonequilibrium Density Operators

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.