We study a model of a quantum collective spin weakly coupled to a spin-polarized Markovian environment and find that the spectrum is divided into two regions that we name normal and exceptional Liouvillian spectral phases. In the thermodynamic limit, the exceptional spectral phase displays the unique property of being made up exclusively of second-order exceptional points. As a consequence, the evolution of any initial density matrix populating this region is slowed down and cannot be described by a linear combination of exponential decays. This phase is separated from the normal one by a critical line in which the density of Liouvillian eigenvalues diverges, a phenomenon analogous to that of excited-state quantum phase transitions observed in some closed quantum systems. In the limit of no bath polarization, this criticality is transferred onto the steady state, implying a dissipative quantum phase transition and the formation of a boundary time crystal.
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.
We consider a collective quantum spin s in contact with Markovian spin-polarized baths. Using a conserved superoperator charge, a differential representation of the Liouvillian is constructed to find its exact spectrum and eigenmodes. We study the spectral properties of the model in the large-s limit using a semiclassical quantization condition and show that the spectral density may diverge along certain curves in the complex plane. We exploit our exact solution to characterize steady-state properties, in particular at the discontinuous phase transition that arises for unpolarized environments, and to determine the decay rates of coherences and populations. Our approach provides a systematic way of finding integrable Liouvillian operators with nontrivial steady states as well as a way to study their spectral properties and eigenmodes.
Non-Hermitian degeneracies of Lindblad generators (Liouvillian exceptional points) can induce non-exponential relaxation and higher-order poles in dynamical response functions. A collective spin coupled to a polarized Markovian bath exhibits an \emph{exceptional spectral phase} in which defective Liouvillian modes imprint super-Lorentzian features in frequency-resolved spectra. We compute the emission spectrum via the Liouvillian resolvent, identify symmetry-sector selection rules, and demonstrate that exceptional-point signatures are strongly state-dependent: they are suppressed in steady-state fluorescence yet become unambiguous for generic (infinite-temperature or random) initial states. Our results provide an experimentally accessible spectroscopic diagnostic of many-body Liouvillian exceptional phases and clarify when steady-state emission can (and cannot) reveal them.
In specific open systems with collective dissipation the Liouvillian can be mapped to a non-Hermitian Hamiltonian. We here consider such a system where the Liouvillian is mapped to an XXZ Richardson-Gaudin integrable model and detail its exact Bethe ansatz solution. While no longer Hermitian, the Hamiltonian is
pseudo-Hermitian/PT symmetric, and as the strength of the coupling to the environment is increased, the spectrum in a fixed symmetry sector changes from a broken pseudo-Hermitian phase with complex conjugate eigenvalues to a pseudo-Hermitian phase with real eigenvalues, passing through a series of exceptional points and associated dissipative quantum phase transitions. The homogeneous limit supports a nontrivial steady state, and away from this limit this state gives rise to a slow logarithmic growth of the decay rate (spectral gap) with system size. Using the exact solution, it is furthermore shown how at large coupling strengths the ratio of the imaginary to the real part of the eigenvalues becomes approximately quantized in the remaining symmetry sectors.
The out of equilibrium dynamics of the Sachdev-Ye-Kitaev model (SYK), comprising
N Majoranas with random all-to-all four-body interactions, minimally coupled to a Markovian bath modeled by the Lindblad formalism, displays intriguing nontrivial features. In particular, the decay rate towards the steady state is a non-monotonic function of the bath coupling
μ, and an analogue of the Loschmidt echo for dissipative quantum systems undergoes a first order dynamical phase transitions that eventually becomes a crossover for sufficiently large
μ. We provide evidence that these features have their origin in the presence of exceptional points in the purely real eigenvalues of the SYK Liouvillian closest to the zero eigenvalue associated with the steady state. An analytic calculation at small
N, supported by numerical results for larger
N, reveals that the value of
μ∼0.1 at which the exceptional point corresponding to the longest living modes occurs is close to a local maximum of the decay rate. This value marks the start of a region of anomalous equilibration where the relaxation rate diminishes as the coupling to the bath becomes stronger. Moreover, the mentioned change from transition to crossover in the Loschmidt echo occurs at a larger
μ∼0.3 corresponding with a proliferation of exceptional points in the low energy limit of the Liouvillian spectrum. We expect these features to be generic in the approach to equilibrium in quantum strongly interacting many-body Liouvillians.
We demonstrate the existence of a dynamical quantum phase transition (DQPT) in a dissipative collective-spin model that exhibits the boundary time crystal (BTC) phase. We initialize the system in the ground state of the Hamiltonian in either the BTC or the non-BTC phase, and drive it across the BTC transition. The driving is done by an abrupt quench or by a finite-time linear ramp of a Hamiltonian control parameter under Markovian Lindblad dynamics. We diagnose DQPTs through zeros of the fidelity-based Loschmidt echo between the initial state and the evolving mixed state, which induce nonanalytic cusp-like features in the associated rate function. For quenches into the BTC phase, the Loschmidt echo exhibits repeated zeros due to the emergent time-periodic steady state, whereas for quenches into the non-BTC phase, the overlap vanishes and remains zero once the dynamics relaxes to a stationary state. We further show that the DQPT persists under the ramp protocol followed by unitary evolution with the final Hamiltonian. Finally, we analyze the finite-size scaling of the first critical time and find convergence to a constant in the thermodynamic limit, with distinct power-law approaches for the quench and the ramp protocols.