Dissipative Phase Transition Research Articles

Dissipative Floquet Dynamics: from Steady State to Measurement Induced Criticality in Trapped-ion Chains

Quantum systems evolving unitarily and subject to quantum measurements exhibit various types of non-equilibrium phase transitions, arising from the competition between unitary evolution and measurements. Dissipative phase transitions in steady states of time-independent Liouvillians and measurement induced phase transitions at the level of quantum trajectories are two primary examples of such transitions. Investigating a many-body spin system subject to periodic resetting measurements, we argue that many-body dissipative Floquet dynamics provides a natural framework to analyze both types of transitions. We show that a dissipative phase transition between a ferromagnetic ordered phase and a paramagnetic disordered phase emerges for long-range systems as a function of measurement probabilities. A measurement induced transition of the entanglement entropy between volume law scaling and sub-volume law scaling is also present, and is distinct from the ordering transition. The two phases correspond to an error-correcting and a quantum-Zeno regimes, respectively. The ferromagnetic phase is lost for short range interactions, while the volume law phase of the entanglement is enhanced. An analysis of multifractal properties of wave function in Hilbert space provides a common perspective on both types of transitions in the system. Our findings are immediately relevant to trapped ion experiments, for which we detail a blueprint proposal based on currently available platforms.

Unveiling quantum entanglement and correlation of sub-Ohmic and Ohmic baths for quantum phase transitions in dissipative systems

By employing the spin-boson model in a dense limit of environmental modes, quantum entanglement and correlation of sub-Ohmic and Ohmic baths for dissipative quantum phase transitions are numerically investigated based on the variational principle. With several measures borrowed from quantum information theory, three different types of singularities are found for the first-order, second-order, and Kosterlitz-Thouless phase transitions, respectively, and the values of transition points and critical exponents are accurately determined. Besides, the scaling form of the quantum discord in the Ohmic case is identified, quite different from that in the sub-Ohmic regime. In a two-spin model, two distinct behaviors of the quantum discord are uncovered: one is related to the quantum entanglement between two spins and the other is decided by the correlation function in the position space rather than the entanglement.

Probing a Dissipative Phase Transition with a Trapped Ion through Reservoir Engineering

Dissipation is often considered as a detrimental effect in quantum systems for unitary quantum operations. However, it has been shown that suitable dissipation can be useful resources in both quantum information and quantum simulation. Here, we propose and experimentally simulate a dissipative phase transition (DPT) model using a single trapped ion with an engineered reservoir. We show that the ion’s spatial oscillation mode reaches a steady state after the alternating application of unitary evolution under a quantum Rabi model Hamiltonian and sideband cooling of the oscillator. The average phonon number of the oscillation mode is used as the order parameter to provide evidence for the DPT. Our work highlights the suitability of trapped ions for simulating open quantum systems and shall facilitate further investigations of DPT with various dissipation terms.

Dissipative dynamics in open XXZ Richardson-Gaudin models

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 $\text{pseudo-Hermitian}/\mathcal{P}\mathcal{T}$ 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.

Liouvillian spectral collapse in the Scully-Lamb laser model

Phase transitions of thermal systems and the laser threshold were first connected more than forty years ago. Despite the nonequilibrium nature of the laser, the Landau theory of thermal phase transitions, applied directly to the Scully-Lamb laser model (SLLM), indicates that the laser threshold is a second-order phase transition, associated with a $U(1)$ spontaneous symmetry breaking (SSB). To capture the genuine nonequilibrium phase transition of the SLLM (i.e., a single-mode laser without a saturable absorber), here we employ a quantum theory of dissipative phase transitions. Our results confirm that the $U(1)$ SSB can occur at the lasing threshold but, in contrast to the Landau theory and semiclassical approximation, they signal that the SLLM "fundamental" transition is a different phenomenon, which we call Liouvillian spectral collapse; that is, the emergence of diabolic points of infinite degeneracy. By considering a generalized SLLM with additional dephasing, we witness a second-order phase transition, with a Liouvillian spectral collapse, but in the absence of symmetry breaking. Most surprisingly, the phase transition corresponds to the emergence of dynamical multistability even without SSB. Normally, bistability is suppressed by quantum fluctuations, while in this case, the very presence of quantum fluctuations enables bistability. This rather anomalous bistability, characterizing the truly dissipative and quantum origin of lasing, can be an experimental signature of our predictions, and we show that it is associated with an emergent dynamical hysteresis.

Dissipative Topological Phase Transition with Strong System-Environment Coupling.

A primary motivation for studying topological matter regards the protection of topological order from its environment. In this work, we study a topological emitter array coupled to an electromagnetic environment. The photon-emitter coupling produces nonlocal interactions between emitters. Using periodic boundary conditions for all ranges of environment-induced interactions, the chiral symmetry inherent to the emitter array is preserved. This chiral symmetry protects the Hamiltonian and induces parity in the Lindblad operator. A topological phase transition occurs at a critical photon-emitter coupling related to the energy spectrum width of the emitter array. Interestingly, the critical point nontrivially changes the dissipation rates of edge states, yielding a dissipative topological phase transition. In the protected topological phase, edge states suffer from environment-induced dissipation for weak photon-emitter coupling. However, strong coupling leads to robust dissipationless edge states with a window at the emitter spacing. Our work shows the potential to manipulate topological quantum matter with electromagnetic environments.

Continuous dissipative phase transitions with or without symmetry breaking

The paradigm of second-order phase transitions (PTs) induced by spontaneous symmetry breaking (SSB) in thermal and quantum systems is a pillar of modern physics that has been fruitfully applied to out-of-equilibrium open quantum systems. Dissipative phase transitions (DPTs) of second order are often connected with SSB, in close analogy with well-known thermal second-order PTs in closed quantum and classical systems. That is, a second-order DPT should disappear by preventing the occurrence of SSB. Here, we prove this statement to be wrong, showing that, surprisingly, SSB is not a necessary condition for the occurrence of second-order DPTs in out-of-equilibrium open quantum systems. We analytically prove this result using the Liouvillian theory of DPTs, and demonstrate this anomalous transition in a paradigmatic laser model, where we can arbitrarily remove SSB while retaining criticality, and on a Z 2-symmetric model of a two-photon Kerr resonator. This new type of PT cannot be interpreted as a ‘semiclassical’ bifurcation, because, after the DPT, the system steady state remains unique.

Determination of the critical exponents in dissipative phase transitions: Coherent anomaly approach

We propose a generalization of the coherent anomaly method to extract the critical exponents of a phase transition occurring in the steady-state of an open quantum many-body system. The method, originally developed by Suzuki [J. Phys. Soc. Jpn. {\bf 55}, 4205 (1986)] for equilibrium systems, is based on the scaling properties of the singularity in the response functions determined through cluster mean-field calculations. We apply this method to the dissipative transverse-field Ising model and the dissipative XYZ model in two dimensions obtaining convergent results already with small clusters.

First-order dissipative phase transition in an exciton-polariton condensate

We investigate the phase diagram of a two-dimensional driven-dissipative system of polaritons coupled to the excitonic reservoir. We find that two critical points exists. The first corresponds to the quasi-condensation and the second to a first-order phase transition from the non-uniform state with spatially modulated density to a uniform state. The latter is related to the modulational instability of a homogeneous state due to the repulsive interactions with the non-condensed reservoir. The first-order character of the transition is evidenced by a discontinuity in the density and the correlation length as well as the phase coexistence and metastability. Moreover, we show that a signature of a Berezinskii-Kosterlitz-Thouless-like transition can be observed in the non-uniform phase.

Metastable quantum entrainment

Quantum van der Pol oscillators are driven-dissipative systems displaying quantum synchronization phenomena. When forced by a squeezed drive, the frequency adjusts to half of the forcing displaying multiple preferred phases. Here we analyze the physical origin of this entrained response, establishing a connection with metastability in open quantum systems. We report a dynamical regime characterized by a huge separation of time scales, in which a dynamical mode displays a lifetime that can be orders of magnitude larger than the rest. In this regime, the long-time dynamics is captured by an incoherent process between two metastable states, which correspond to the preferred phases of the synchronized oscillator. In fact, we show that quantum entrainment is here characterized by fluctuations driving an incoherent process between two metastable phases, which ultimately limits its temporal coherence when moving into the quantum regime. Finally, we discuss connections with the phenomena of dissipative phase transitions and transient synchronization in open quantum systems.

Subradiant-to-Subradiant Phase Transition in the Bad Cavity Laser.

We show that the onset of steady-state superradiance in a bad cavity laser is preceded by a dissipative phase transition between two distinct phases of steady-state subradiance. The transition is marked by a nonanalytic behavior of the cavity output power and the mean atomic inversion, as well as a discontinuity in the variance of the collective atomic inversion. In particular, for repump rates below a critical value, the cavity output power is strongly suppressed and does not increase with the atom number, while it scales linearly with atom number above this value. Remarkably, we find that the atoms are in a macroscopically entangled steady state near the critical region with a vanishing fraction of unentangled atoms in the large atom number limit.

Theory of classical metastability in open quantum systems

We present a general theory of classical metastability in open quantum systems. Metastability is a consequence of a large separation in timescales in the dynamics, leading to the existence of a regime when states of the system appear stationary, before eventual relaxation toward a true stationary state at much larger times. In this work, we focus on the emergence of classical metastability, i.e., when metastable states of an open quantum system with separation of timescales can be approximated as probabilistic mixtures of a finite number of states. We find that a number of classical features follow from this approximation, for the manifold of metastable states, long-time dynamics between them, and symmetries of the dynamics. Namely, those states are approximately disjoint and thus play the role of metastable phases, the relaxation toward the stationary state is approximated by a classical stochastic dynamics between them, and weak symmetries correspond to their permutations. Importantly, the classical dynamics is observed not only on average but also at the level of individual quantum trajectories: We show that time coarse-grained continuous measurement records can be viewed as noisy classical trajectories, while their statistics can be approximated by that of the classical dynamics. Among others, this explains how first-order dynamical phase transitions arise from metastability. Finally, to verify the presence of classical metastability in a given open quantum system, we develop an efficient numerical approach that delivers the set of metastable phases together with the effective classical dynamics. Since the proximity to a first-order dissipative phase transition manifests as metastability, the theory and tools introduced in this work can be used to investigate such transitions---which occur in the large size limit---through the metastable behavior of many-body systems of moderate sizes accessible to numerics.

Quantum criticality of the Ohmic spin-boson model in a high dense spectrum: Symmetries, quantum fluctuations and correlations

Study of dissipative quantum phase transitions in the Ohmic spin-boson model is numerically challenging in a dense limit of environmental modes. In this work, large-scale numerical simulations are carried out based on the variational principle. The validity of variational calculations, spontaneous breakdown of symmetries, and quantum fluctuations and correlations in the Ohmic bath are carefully analyzed, and the critical coupling as well as exponents are accurately determined in the weak tunneling and continuum limits. In addition, quantum criticality of the Ohmic bath is uncovered both in the delocalized phase and at the transition point.

Coherence of resonant light-matter interaction in the strong-coupling limit

We explore the role of quantum fluctuations in the strong-coupling limit of the dissipative Jaynes–Cummings oscillator driven on resonance. For weak excitation, we derive analytical expressions for the spectrum and the intensity correlation function for the photons scattered by the two-state atom coupled to the coherently driven cavity mode. We do so by writing down a birth–death process adding the higher orders in the excitation strength needed to go beyond the pure-state factorization, following the method introduced in Carmichael (2008). Our results for the first and second-order correlation functions are complemented by the numerical investigation of the waiting-time distribution for the photon emissions directed sideways, and the comparison with ordinary resonance fluorescence. To close out our discussion, we increase the driving field amplitude and approach the critical point organizing a second-order dissipative quantum phase transition by depicting the excitation pathways in the intracavity field distribution for a finite system size.

Solving the Liouvillian Gap with Artificial Neural Networks.

We propose a machine-learning inspired variational method to obtain the Liouvillian gap, which plays a crucial role in characterizing the relaxation time and dissipative phase transitions of open quantum systems. By using "spin bi-base mapping," we map the density matrix to a pure restricted-Boltzmann-machine (RBM) state and transform the Liouvillian superoperator to a rank-two non-Hermitian operator. The Liouvillian gap can be obtained by a variational real-time evolution algorithm under this non-Hermitian operator. We apply our method to the dissipative Heisenberg model in both one and two dimensions. For the isotropic case, we find that the Liouvillian gap can be analytically obtained and in one dimension even the whole Liouvillian spectrum can be exactly solved using the Bethe ansatz method. By comparing our numerical results with their analytical counterparts, we show that the Liouvillian gap could be accessed by the RBM approach efficiently to a desirable accuracy, regardless of the dimensionality and entanglement properties.

Time- and frequency-domain two-particle correlations of a driven dissipative Bose-Hubbard model

We investigate theoretically the time- and frequency-domain two-particle correlations of a driven-dissipative Bose-Hubbard model at and near a dissipative phase transition (DPT). We compute the Hanbury Brown--Twiss (HBT)-type two-particle temporal correlation function ${g}^{2}(\ensuremath{\tau})$ which, as a function of time delay $\ensuremath{\tau}$, exhibits oscillations with frequencies determined by the imaginary part of the Liouvillian gap. As the gap closes near a transition point, the oscillations at that point die down. For parameters slightly away from the transition point, the HBT correlations show oscillations from superbunching to antibunching regimes. We show that the Fourier transform of HBT correlations into the frequency domain provide information about DPT and Liouvillian dynamics. We numerically solve the many-body Lindblad master equation and calculate the Wigner distribution of the system in the steady state to ascertain the DPT. Below a certain drive strength, the Fourier transform shows a two-peak structure, while above that strength it exhibits either a Lorentzian-like single-peak structure or a structure with two dips. The width of the single-peak structure is minimal at the phase-transition point and the peak of this structure always lies at zero frequency. The positions of the two symmetrical peaks in case of a two-peak structure are given by the imaginary parts of the Liouvillian gap while their half width at half maximum is given by the real part of the gap. The positions and widths of the two dips are also related to low-lying eigenvalues of the Liouvillian operator. We discuss quantum statistical properties of the model in terms of the HBT correlation function and its Fourier transform.

Limit cycle phase and Goldstone mode in driven dissipative systems

In this article, we theoretically investigate the first- and second-order quantum dissipative phase transitions of a three-mode cavity with a Hubbard interaction. In both types, there is a mean-field (MF) limit cycle phase where the local U(1) symmetry and the time-translational symmetry of the Liouvillian superoperator are spontaneously broken. In MF, this spontaneous symmetry breaking manifests itself through the appearance of an unconditionally and fully squeezed state at the cavity output, connected to the well-known Goldstone mode. By employing the Wigner function formalism, hence, properly including the quantum noise, we show that away from the thermodynamic limit and within the quantum regime, fluctuations notably limit the coherence time of the Goldstone mode due to the phase diffusion. Our theoretical predictions suggest that interacting multimode photonic systems are rich, versatile test beds for investigating the crossovers between the mean-field picture and quantum phase transitions, a problem that can be investigated in various platforms including superconducting circuits, semiconductor microcavities, atomic Rydberg polaritons, and cuprite excitons.

Application of non-Hermitian Hamiltonian model in open quantum optical systems* *Project supported by the National Natural Science Foundation of China (Grant Nos. 61925503, 11874038, and 11654002), the Key Project of the National Key R&D Program of China (Grant Nos. 2016YFA0301402 and 2020YFA0309400), the Program for the Innovative Talents of Higher Education Institutions of Shanxi, the Program for Sanjin Scholars of

Non-Hermitian systems have observed numerous novel phenomena and might lead to various applications. Unlike standard quantum physics, the conservation of energy guaranteed by the closed system is broken in the non-Hermitian system, and the energy can be exchanged between the system and the environment. Here we present a scheme for simulating the dissipative phase transition with an open quantum optical system. The competition between the coherent interaction and dissipation leads to the second-order phase transition. Furthermore, the quantum correlation in terms of squeezing is studied around the critical point. Our work may provide a new route to explore the non-Hermitian quantum physics with feasible techniques in experiments.

Dissipative phase transitions in the fully connected Ising model with p -spin interaction

In this paper, we study the driven-dissipative p-spin models for $p\geq 2$. In thermodynamics limit, the equation of motion is derived by using a semiclassical approach. The long-time asymptotic states are obtained analytically, which exhibit multi-stability in some regions of the parameter space. The steady state is unique as the number of spins is finite. But the thermodynamic limit of the steady-state magnetization displays nonanalytic behavior somewhere inside the semiclassical multi-stable region. We find both the first-order and continuous dissipative phase transitions. As the number of spins increases, both the Liouvillian gap and magnetization variance vanish according to a power law at the continuous transition. At the first-order transition, the gap vanishes exponentially accompanied by a jump of magnetization in thermodynamic limit. The properties of transitions depend on the symmetry and semiclassical multistability, being qualitatively different among $p=2$, odd $p$ ($p\geq 3$) and even $p$ ($p\geq 4$).

Bethe ansatz approach for dissipation: exact solutions of quantum many-body dynamics under loss

We develop a Bethe ansatz based approach to study dissipative systems experiencing loss. The method allows us to exactly calculate the spectra of interacting, many-body Liouvillians. We discuss how the dissipative Bethe ansatz opens the possibility of analytically calculating the dynamics of a wide range of experimentally relevant models including cold atoms subjected to one and two body losses, coupled cavity arrays with bosons escaping the cavity, and cavity quantum electrodynamics. As an example of our approach we study the relaxation properties in a boundary driven XXZ spin chain. We exactly calculate the Liouvillian gap and find different relaxation rates with a novel type of dynamical dissipative phase transition. This physically translates into the formation of a stable domain wall in the easy-axis regime despite the presence of loss. Such analytic results have previously been inaccessible for systems of this type.