Spectral Gaps and Midgap States in Random Quantum Master Equations

Dynamical Transitions from Slow to Fast Relaxation in Random Open Quantum Systems

We explore the effects of spatial locality on the dynamics of random quantum systems subject to a Markovian noise. To this end, we study a model in which the system Hamiltonian and its couplings to the noise are random matrices whose entries decay as power laws of distance, with distinct exponents α_{H}, α_{L}. The steady state is always featureless, but the rate at which it is approached exhibits three phases depending on α_{H} and α_{L}: a phase where the approach is asymptotically exponential as a result of a gap in the spectrum of the Lindblad superoperator that generates the dynamics, and two gapless phases with subexponential relaxation, distinguished by the manner in which the gap decreases with system size. Within perturbation theory, the phase boundaries in the (α_{H},α_{L}) plane differ for weak and strong decoherence, suggesting phase transitions as a function of noise strength. We identify nonperturbative effects that prevent such phase transitions in the thermodynamic limit.

Random sparse generators of Markovian evolution and their spectral properties

The evolution of a complex multistate system is often interpreted as a continuous-time Markovian process. To model the relaxation dynamics of such systems, we introduce an ensemble of random sparse matrices which can be used as generators of Markovian evolution. The sparsity is controlled by a parameter φ, which is the number of nonzero elements per row and column in the generator matrix. Thus, a member of the ensemble is characterized by the Laplacian of a directed regular graph with D vertices (number of system states) and 2φD edges with randomly distributed weights. We study the effects of sparsity on the spectrum of the generator. Sparsity is shown to close the large spectral gap that is characteristic of nonsparse random generators. We show that the first moment of the eigenvalue distribution scales as ∼φ, while its variance is ∼sqrt[φ]. By using extreme value theory, we demonstrate how the shape of the spectral edges is determined by the tails of the corresponding weight distributions and clarify the behavior of the spectral gap as a function of D. Finally, we analyze complex spacing ratio statistics of ultrasparse generators, φ=const, and find that starting already at φ⩾2, spectra of the generators exhibit universal properties typical of Ginibre's orthogonal ensemble.

Entanglement Entropy of Non-Hermitian Eigenstates and the Ginibre Ensemble

Entanglement entropy is a powerful tool in characterizing universal features in quantum many-body systems. In quantum chaotic Hermitian systems, typical eigenstates have near maximal entanglement with very small fluctuations. Here, we show that for Hamiltonians displaying non-Hermitian many-body quantum chaos, modeled by the Ginibre ensemble, the entanglement entropy of typical eigenstates is greatly suppressed. The entropy does not grow with the Hilbert space dimension for sufficiently large systems, and the fluctuations are of equal order. We derive the novel entanglement spectrum that has infinite support in the complex plane and strong energy dependence. We provide evidence of universality, and similar behavior is found in the non-Hermitian Sachdev-Ye-Kitaev model, indicating the general applicability of the Ginibre ensemble to dissipative many-body quantum chaos.

Spacing distribution in the two-dimensional Coulomb gas: Surmise and symmetry classes of non-Hermitian random matrices at noninteger β

A random matrix representation is proposed for the two-dimensional (2D) Coulomb gas at inverse temperature β. For 2×2 matrices with Gaussian distribution we analytically compute the nearest-neighbor spacing distribution of complex eigenvalues in radial distance. Because it does not provide such a good approximation as the Wigner surmise in 1D, we introduce an effective β_{eff}(β) in our analytic formula that describes the spacing obtained numerically from the 2D Coulomb gas well for small values of β. It reproduces the 2D Poisson distribution at β=0 exactly, that is valid for a large particle number. The surmise is used to fit data in two examples, from open quantum spin chains and ecology. The spacing distributions of complex symmetric and complex quaternion self-dual ensembles of non-Hermitian random matrices, that are only known numerically, are very well fitted by noninteger values β=1.4 and β=2.6 from a 2D Coulomb gas, respectively. These two ensembles have been suggested as the only two symmetry classes, where the 2D bulk statistics is different from the Ginibre ensemble.

Spectral Filtering Induced by Non-Hermitian Evolution with Balanced Gain and Loss: Enhancing Quantum Chaos

The dynamical signatures of quantum chaos in an isolated system are captured by the spectral form factor, which exhibits as a function of time a dip, a ramp, and a plateau, with the ramp being governed by the correlations in the level spacing distribution. While decoherence generally suppresses these dynamical signatures, the nonlinear non-Hermitian evolution with balanced gain and loss (BGL) in an energy-dephasing scenario can enhance manifestations of quantum chaos. In the Sachdev-Ye-Kitaev model and random matrix Hamiltonians, BGL increases the span of the ramp, lowering the dip as well as the value of the plateau, providing an experimentally realizable physical mechanism for spectral filtering. The chaos enhancement due to BGL is optimal over a family of filter functions that can be engineered with fluctuating Hamiltonians.

Spectral Statistics of Non-Hermitian Matrices and Dissipative Quantum Chaos

We propose a measure, which we call the dissipative spectral form factor (DSFF), to characterize the spectral statistics of non-Hermitian (and nonunitary) matrices. We show that DSFF successfully diagnoses dissipative quantum chaos and reveals correlations between real and imaginary parts of the complex eigenvalues up to arbitrary energy scale (and timescale). Specifically, we provide the exact solution of DSFF for the complex Ginibre ensemble (GinUE) and for a Poissonian random spectrum (Poisson) as minimal models of dissipative quantum chaotic and integrable systems, respectively. For dissipative quantum chaotic systems, we show that the DSFF exhibits an exact rotational symmetry in its complex time argument τ. Analogous to the spectral form factor (SFF) behavior for Gaussian unitary ensemble, the DSFF for GinUE shows a "dip-ramp-plateau" behavior in |τ|: the DSFF initially decreases, increases at intermediate timescales, and saturates after a generalized Heisenberg time, which scales as the inverse mean level spacing. Remarkably, for large matrix size, the "ramp" of the DSFF for GinUE increases quadratically in |τ|, in contrast to the linear ramp in the SFF for Hermitian ensembles. For dissipative quantum integrable systems, we show that the DSFF takes a constant value, except for a region in complex time whose size and behavior depend on the eigenvalue density. Numerically, we verify the above claims and additionally show that the DSFF for real and quaternion real Ginibre ensembles coincides with the GinUE behavior, except for a region in the complex time plane of measure zero in the limit of large matrix size. As a physical example, we consider the quantum kicked top model with dissipation and show that it falls under the Ginibre universality class and Poisson as the "kick" is switched on or off. Lastly, we study spectral statistics of ensembles of random classical stochastic matrices or Markov chains and show that these models again fall under the Ginibre universality class.

Random generators of Markovian evolution: A quantum-classical transition by superdecoherence

Continuous-time Markovian evolution appears to be manifestly different in classical and quantum worlds. We consider ensembles of random generators of N-dimensional Markovian evolution, quantum and classical ones, and evaluate their universal spectral properties. We then show how the two types of generators can be related by superdecoherence. In analogy with the mechanism of decoherence, which transforms a quantum state into a classical one, superdecoherence can be used to transform a Lindblad operator (generator of quantum evolution) into a Kolmogorov operator (generator of classical evolution). We inspect spectra of random Lindblad operators undergoing superdecoherence and demonstrate that, in the limit of complete superdecoherence, the resulting operators exhibit spectral density typical to random Kolmogorov operators. By gradually increasing strength of superdecoherence, we observe a sharp quantum-to-classical transition. Furthermore, we define an inverse procedure of supercoherification that is a generalization of the scheme used to construct a quantum state out of a classical one. Finally, we study microscopic correlation between neighboring eigenvalues through the complex spacing ratios and observe the horseshoe distribution, emblematic of the Ginibre universality class, for both types of random generators. Remarkably, it survives both superdecoherence and supercoherification.

Full Spectrum of the Liouvillian of Open Dissipative Quantum Systems in the Zeno Limit

We consider an open quantum system with dissipation, described by a Lindblad Master equation (LME). For dissipation locally acting and sufficiently strong, a separation of the relaxation timescales occurs, which, in terms of the eigenvalues of the Liouvillian, implies a grouping of the latter in distinct vertical stripes in the complex plane at positions determined by the eigenvalues of the dissipator. We derive effective LME equations describing the modes within each stripe separately, and solve them perturbatively, obtaining for the full set of eigenvalues and eigenstates of the Liouvillian explicit expressions correct at order 1/Γ included, where Γ is the strength of the dissipation. As an example, we apply our general results to quantum XYZ spin chains coupled, at one boundary, to a dissipative bath of polarization.

Hierarchy of Relaxation Timescales in Local Random Liouvillians

To characterize the generic behavior of open quantum systems, we consider random, purely dissipative Liouvillians with a notion of locality. We find that the positivity of the map implies a sharp separation of the relaxation timescales according to the locality of observables. Specifically, we analyze a spin-1/2 system of size ℓ with up to n-body Lindblad operators, which are n local in the complexity-theory sense. Without locality (n=ℓ), the complex Liouvillian spectrum densely covers a "lemon"-shaped support, in agreement with recent findings [S. Denisov et al., Phys. Rev. Lett. 123, 140403 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.140403]. However, for local Liouvillians (n<ℓ), we find that the spectrum is composed of several dense clusters with random matrix spacing statistics, each featuring a lemon-shaped support wherein all eigenvectors correspond to n-body decay modes. This implies a hierarchy of relaxation timescales of n-body observables, which we verify to be robust in the thermodynamic limit. Our findings for n locality generalize immediately to the case of spatial locality, introducing further splitting of timescales due to the additional structure.

Tenfold Way for Quadratic Lindbladians

We uncover a topological classification applicable to open fermionic systems governed by a general class of Lindblad master equations. These "quadratic Lindbladians" can be captured by a non-Hermitian single-particle matrix which describes internal dynamics as well as system-environment coupling. We show that this matrix must belong to one of ten non-Hermitian Bernard-LeClair symmetry classes which reduce to the Altland-Zirnbauer classes in the closed limit. The Lindblad spectrum admits a topological classification, which we show results in gapless edge excitations with finite lifetimes. Unlike previous studies of purely Hamiltonian or purely dissipative evolution, these topological edge modes are unconnected to the form of the steady state. We provide one-dimensional examples where the addition of dissipators can either preserve or destroy the closed classification of a model, highlighting the sensitivity of topological properties to details of the system-environment coupling.