We explore the connections between dissipative quantum phase transitions and non-Hermitian random matrix theory. For this, we work in the framework of the dissipative Dicke model which is archetypal of symmetry-breaking phase transitions in open quantum systems. We establish that the Liouvillian describing the quantum dynamics exhibits distinct spectral features of integrable and chaotic character on the two sides of the critical point. We follow the distribution of the spacings of the complex Liouvillian eigenvalues across the critical point. In the normal and superradiant phases, the distributions are two-dimensional Poisson and that of the Ginibre unitary random matrix ensemble, respectively. Our results are corroborated by computing a recently introduced complex-plane generalization of the consecutive level-spacing ratio distribution. Our approach can be readily adapted for classifying the nature of quantum dynamics across dissipative critical points in other open quantum systems.
We explore the connections between dissipative quantum phase transitions and non-Hermitian random matrix theory. For this, we work in the framework of the dissipative Dicke model which is archetypal of symmetry-breaking phase transitions in open quantum systems. We establish that the Liouvillian describing the quantum dynamics exhibits distinct spectral features of integrable and chaotic character on the two sides of the critical point. We follow the distribution of the spacings of the complex Liouvillian eigenvalues across the critical point. In the normal and superradiant phases, the distributions are two-dimensional Poisson and that of the Ginibre unitary random matrix ensemble, respectively. Our results are corroborated by computing a recently introduced complex-plane generalization of the consecutive level-spacing ratio distribution. Our approach can be readily adapted for classifying the nature of quantum dynamics across dissipative critical points in other open quantum systems.
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.
Motivated by the surge of interest in statistics of non-Hermitian random matrices as a framework for description of universal characteristics of dissipative chaotic quantum many-body systems, we address the problem of characterizing the parametric correlations of spectral densities. Considering parameter-dependent ensemble of complex Ginibre matrices we derive an explicit, closed-form expression for the parametric number covariance in the systems of symmetry class A for eigenvalues in a circular domain containing on average a finite number of eigenvalues in the spectral bulk. This behavior is expected to be universal, as further supported by numerical evidence for the real Ginibre ensemble, non-Hermitian Bernoulli Wigner matrices and bi-unitarily invariant ensembles. We also discuss a relation between parametric correlations of spectral densities and the distribution of the so-called eigenvector non-orthogonality factor, which attracted considerable interest in recent years.
A state of an open quantum system is described by a density matrix, whose dynamics is governed by a Liouvillian superoperator. Within a general framework, we explore fundamental properties of both first-order dissipative phase transitions and second-order dissipative phase transitions associated with a symmetry breaking. In the critical region, we determine the general form of the steady-state density matrix and of the Liouvillian eigenmatrix whose eigenvalue defines the Liouvillian spectral gap. We illustrate our exact results by studying some paradigmatic quantum optical models exhibiting critical behavior.
The complex spectral statistics of non-Hermitian random matrices provide a standard diagnostic of dissipative quantum chaos. Yet an analytical understanding of how symmetry enriches these universal correlations remains limited. We establish an exact duality between the spectral statistics of Hermitian and non-Hermitian random matrices in the limit of large matrix size. This duality resolves a longstanding analytical problem by yielding the previously unknown eigenvalue pair-correlation functions for complex symmetric (class AI†) and complex self-dual matrices (class AII†), completing the analytical description of bulk pair correlations in the non-Hermitian threefold way. It further gives closed-form spectral densities near the origin for seven additional symmetry classes. Numerical tests in interacting and dissipative quantum models support the universality of these predictions beyond Gaussian ensembles. The duality suggests a deeper connection between chaos in closed and open quantum systems.
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.
We study the transition between integrable and chaotic behavior in dissipative open quantum systems, exemplified by a boundary driven quantum spin chain. The repulsion between the complex eigenvalues of the corresponding Liouville operator in radial distance s is used as a universal measure. The corresponding level spacing distribution is well fitted by that of a static two-dimensional Coulomb gas with harmonic potential at inverse temperature β∈[0,2]. Here, β=0 yields the two-dimensional Poisson distribution, matching the integrable limit of the system, and β=2 equals the distribution obtained from the complex Ginibre ensemble, describing the fully chaotic limit. Our findings generalize the results of Grobe, Haake, and Sommers, who derived a universal cubic level repulsion for small spacings s. We collect mathematical evidence for the universality of the full level spacing distribution in the fully chaotic limit at β=2. It holds for all three Ginibre ensembles of random matrices with independent real, complex, or quaternion matrix elements.