The study of dissipation and decoherence in generic open quantum systems recently led to the investigation of spectral and steady-state properties of random Lindbladian dynamics. A natural question is then how realistic and universal those properties are. Here, we address these issues by considering a different description of dissipative quantum systems, namely, the discrete-time Kraus map representation of completely positive quantum dynamics. Through random matrix theory (RMT) techniques and numerical exact diagonalization, we study random Kraus maps, allowing for a varying dissipation strength, and their local circuit counterpart. We find the spectrum of the random Kraus map to be either an annulus or a disk inside the unit circle in the complex plane, with a transition between the two cases taking place at a critical value of dissipation strength. The eigenvalue distribution and the spectral transition are well described by a simplified RMT model that we can solve exactly in the thermodynamic limit, by means of non-Hermitian RMT and quaternionic free probability. The steady state, on the contrary, is not affected by the spectral transition. It has, however, a perturbative crossover regime at small dissipation, inside which the steady state is characterized by uncorrelated eigenvalues. At large dissipation (or for any dissipation for a large-enough system), the steady state is well described by a random Wishart matrix. The steady-state properties thus coincide with those already observed for random Lindbladian dynamics, indicating their universality. Quite remarkably, the statistical properties of the local Kraus circuit are qualitatively the same as those of the nonlocal Kraus map, indicating that the latter, which is more tractable, already captures the realistic and universal physical properties of generic open quantum systems.
The study of dissipation and decoherence in generic open quantum systems recently led to the investigation of spectral and steady-state properties of random Lindbladian dynamics. A natural question is then how realistic and universal those properties are. Here, we address these issues by considering a different description of dissipative quantum systems, namely, the discrete-time Kraus map representation of completely positive quantum dynamics. Through random matrix theory (RMT) techniques and numerical exact diagonalization, we study random Kraus maps, allowing for a varying dissipation strength, and their local circuit counterpart. We find the spectrum of the random Kraus map to be either an annulus or a disk inside the unit circle in the complex plane, with a transition between the two cases taking place at a critical value of dissipation strength. The eigenvalue distribution and the spectral transition are well described by a simplified RMT model that we can solve exactly in the thermodynamic limit, by means of non-Hermitian RMT and quaternionic free probability. The steady state, on the contrary, is not affected by the spectral transition. It has, however, a perturbative crossover regime at small dissipation, inside which the steady state is characterized by uncorrelated eigenvalues. At large dissipation (or for any dissipation for a large-enough system), the steady state is well described by a random Wishart matrix. The steady-state properties thus coincide with those already observed for random Lindbladian dynamics, indicating their universality. Quite remarkably, the statistical properties of the local Kraus circuit are qualitatively the same as those of the nonlocal Kraus map, indicating that the latter, which is more tractable, already captures the realistic and universal physical properties of generic open quantum systems.
We study spectral and steady-state properties of generic Markovian dissipative systems described by quadratic fermionic Liouvillian operators of the Lindblad form. The Hamiltonian dynamics is modeled by a generic random quadratic operator, i.e., as a featureless superconductor of class D, whereas the Markovian dissipation is described by M M random linear jump operators. By varying the dissipation strength and the ratio of dissipative channels per fermion, m=M/(2N_F) m = M / ( 2 N F ) , we find two distinct phases where the support of the single-particle spectrum has one or two connected components. In the strongly dissipative regime, this transition occurs for m=1/2 m = 1 / 2 and is concomitant with a qualitative change in both the steady-state and the spectral gap that rules the large-time dynamics. Above this threshold, the spectral gap and the steady-state purity qualitatively agree with the fully generic (i.e., non-quadratic) case studied recently. Below m=1/2 m = 1 / 2 , the spectral gap closes in the thermodynamic limit and the steady-state decouples into an ergodic and a nonergodic sector yielding a non-monotonic steady-state purity as a function of the dissipation strength. Our results show that some of the universal features previously observed for fully random Liouvillians are generic for a sufficiently large number of jump operators. On the other hand, if the number of dissipation channels is decreased the system can exhibit nonergodic features, rendering it possible to suppress dissipation in protected subspaces even in the presence of strong system-environment coupling.
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.
Understanding the far-from-equilibrium dynamics of dissipative quantum systems, where dissipation and decoherence coexist with unitary dynamics, is an enormous challenge with immense rewards. Often, the only realistic approach is to forgo a detailed microscopic description and search for signatures of universal behavior shared by collections of many distinct, yet sufficiently similar, complex systems. Quantum chaos provides a powerful statistical framework for addressing this question, relying on symmetries to obtain information not accessible otherwise. This thesis examines how to reconcile chaos with dissipation, proceeding along two complementary lines. In Part I, we apply non-Hermitian random matrix theory to open quantum systems with Markovian dissipation and discuss the relaxation timescales and steady states of three representative examples of increasing physical relevance: single-particle Lindbladians and Kraus maps, open free fermions, and dissipative Sachdev-Ye-Kitaev (SYK) models. In Part II, we investigate the symmetries, correlations, and universality of many-body open quantum systems, classifying several models of dissipative quantum matter. From a theoretical viewpoint, this thesis lays out a generic framework for the study of the universal properties of realistic, chaotic, and dissipative quantum systems. From a practical viewpoint, it provides the concrete building blocks of dynamical dissipative evolution constrained by symmetry, with potential technological impact on the fabrication of complex quantum structures. (Full abstract in the thesis.)
We numerically study the spectral statistics of open quantum many-body systems (OQMBS) as signatures of quantum chaos (or the lack thereof), using the dissipative spectral form factor (DSFF), a generalization of the spectral form factor to complex spectra. We show that the DSFF of chaotic OQMBS generically displays the quadratic ramp-plateau behaviour of the Ginibre ensemble from random matrix theory, in contrast to the linear ramp-plateau behaviour of the Gaussian ensemble in closed quantum systems. Furthermore, in the presence of many-body interactions, such RMT behaviour emerges only after a time scale τdev, which generally increases with system size for sufficiently large system size, and can be identified as the non-Hermitian analogue of the many-body Thouless time. The universality of the random matrix theory behavior is demonstrated by surveying twelve models of OQMBS, including random Kraus circuits (quantum channels) and random Lindbladians (Liouvillians) in several symmetry classes, as well as Lindbladians of paradigmatic models such as the Sachdev-Ye-Kitaev (SYK), XXZ, and the transverse field Ising models. We devise an unfolding and filtering procedure to remove variations of the averaged density of states which would otherwise hide the universal RMT-like signatures in the DSFF for chaotic OQMBS. Beyond chaotic OQMBS, we study the spectral statistics of non-chaotic OQMBS, specifically the integrable XX model and a system in the many-body localized (MBL) regime in the presence of dissipation, which exhibit DSFF behaviours distinct from the ramp-plateau behaviour of random matrix theory. Lastly, we study the DSFF of Lindbladians with the Hamiltonian term set to zero, i.e. only the jump operators are present, and demonstrate that the results of RMT universality and scaling of many-body Thouless time survive even without coherent evolution.
Generators of continuous-time quantum Markovian evolution, known as Lindblad operators, play a key role in the theory of open quantum systems. The question of how to define and characterize random Lindblad operators has attracted considerable attention in recent years. Purely dissipative random Lindbladians have been extensively studied using toolboxes of free probability and non-Hermitian random matrix theories, which helped to uncover a characteristic lemon-like spectral support. In this work, we focus on random Lindbladians with a Hamiltonian component. We analyze three ensembles of increasing complexity, designed to illustrate the joint effects of dissipative and Hamiltonian dynamics. While the first two allow for an exact analytic treatment with standard probabilistic and scalar-valued free-probability methods, the most general case requires the more general framework of operator-valued free probability. We explain how the spectral support changes with the relative Hamiltonian strength and the rank of the dissipative component and characterize the effect of Hamiltonian-induced spectral squeezing.
Understanding the dynamics of open quantum systems is intrinsically important to control and utilize quantum hardware in various applications. The generic behavior of purely dissipative open quantum many-body systems with k-local, normal dissipation processes can be investigated using random matrix theory, revealing different decay timescales of observables determined by their complexity, as shown by Wang et al. [K. Wang et al., Phys. Rev. Lett. 124, 100604 (2020)]. The time evolution of an open quantum system is entirely described by the eigenvalue spectrum of the Lindbladian, in which the hierarchy of decay timescales is reflected in the formation of well-separated eigenvalue clusters. Here, we extend the analysis of this spectrum to the presence of unitary dynamics using a random matrix model. In the case of strong dissipation, the unitary dynamics can be treated perturbatively and the k locality of the Hamiltonian determines how susceptible the spectrum is to such a perturbation. For the physically most relevant case of (dissipative) two-body interactions, the correction in the first order of the perturbation vanishes, leading to the relative robustness of the spectral features, which in turn implies the same robustness of the corresponding separation of decay timescales. We find the perturbative approach for the spectrum to yield an excellent agreement with exact-diagonalization results. For weak dissipation, the spectrum flows into clusters except for outlier eigenmodes, which we identify to be the local symmetries of the Hamiltonian. We support our analytic findings by numerical simulations and show that our generic model applies to a nonrandom Heisenberg chain.