Exact Bethe Ansatz Spectrum of a Tight-Binding Chain with Dephasing Noise

Transport and Nonreciprocity in Monitored Quantum Devices: An Exact Study

We study noninteracting fermionic systems undergoing continuous monitoring and driven by biased reservoirs. Averaging over the measurement outcomes, we derive exact formulas for the particle and heat flows in the system. We show that these currents feature competing elastic and inelastic components, which depend nontrivially on the monitoring strength γ. We highlight that monitor-induced inelastic processes lead to nonreciprocal currents, allowing one to extract work from measurements without active feedback control. We illustrate our formalism with two distinct monitoring schemes providing measurement-induced power or cooling. Optimal performances are found for values of the monitoring strength γ, which are hard to address with perturbative approaches.

Lieb-Schultz-Mattis Theorem in Open Quantum Systems

The Lieb-Schultz-Mattis (LSM) theorem provides a general constraint on quantum many-body systems and plays a significant role in the Haldane gap phenomena and topological phases of matter. Here, we extend the LSM theorem to open quantum systems and establish a general theorem that restricts the steady state and spectral gap of Liouvillians based solely on symmetry. Specifically, we demonstrate that the unique gapped steady state is prohibited when translation invariance and U(1) symmetry are simultaneously present for noninteger filling numbers. As an illustrative example, we find that no dissipative gap is open in the spin-1/2 dissipative Heisenberg model, while a dissipative gap can be open in the spin-1 counterpart-an analog of the Haldane gap phenomena in open quantum systems. Furthermore, we show that the LSM constraint manifests itself in a quantum anomaly of the dissipative form factor of Liouvillians. We also find the LSM constraints due to symmetry intrinsic to open quantum systems, such as Kubo-Martin-Schwinger symmetry. Our work leads to a unified understanding of phases and phenomena in open quantum systems.

Anomalous Diffusion in the Long-Range Haken-Strobl-Reineker Model

We analyze the propagation of excitons in a d-dimensional lattice with power-law hopping ∝1/r^{α} in the presence of dephasing, described by a generalized Haken-Strobl-Reineker model. We show that in the strong dephasing (quantum Zeno) regime the dynamics is described by a classical master equation for an exclusion process with long jumps. In this limit, we analytically compute the spatial distribution, whose shape changes at a critical value of the decay exponent α_{cr}=(d+2)/2. The exciton always diffuses anomalously: a superdiffusive motion is associated to a Lévy stable distribution with long-range algebraic tails for α≤α_{cr}, while for α>α_{cr} the distribution corresponds to a surprising mixed Gaussian profile with long-range algebraic tails, leading to the coexistence of short-range diffusion and long-range Lévy flights. In the many-exciton case, we demonstrate that, starting from a domain-wall exciton profile, algebraic tails appear in the distributions for any α, which affects thermalization: the longer the hopping range, the faster equilibrium is reached. Our results are directly relevant to experiments with cold trapped ions, Rydberg atoms, and supramolecular dye aggregates. They provide a way to realize an exclusion process with long jumps experimentally.

Relaxation to a Parity-Time Symmetric Generalized Gibbs Ensemble after a Quantum Quench in a Driven-Dissipative Kitaev Chain

The construction of the generalized Gibbs ensemble, to which isolated integrable quantum many-body systems relax after a quantum quench, is based upon the principle of maximum entropy. In contrast, there are no universal and model-independent laws that govern the relaxation dynamics and stationary states of open quantum systems, which are subjected to Markovian drive and dissipation. Yet, as we show, relaxation of driven-dissipative systems after a quantum quench can, in fact, be determined by a maximum entropy ensemble, if the Liouvillian that generates the dynamics of the system has parity-time symmetry. Focusing on the specific example of a driven-dissipative Kitaev chain, we show that, similar to isolated integrable systems, the approach to a parity-time symmetric generalized Gibbs ensemble becomes manifest in the relaxation of local observables and the dynamics of subsystem entropies. In contrast, the directional pumping of fermion parity, which is induced by nontrivial non-Hermitian topology of the Kitaev chain, represents a phenomenon that is unique to relaxation dynamics in driven-dissipative systems. Upon increasing the strength of dissipation, parity-time symmetry is broken at a finite critical value, which thus constitutes a sharp dynamical transition that delimits the applicability of the principle of maximum entropy. We show that these results, which we obtain for the specific example of the Kitaev chain, apply to broad classes of noninteracting fermionic models, and we discuss their generalization to a noninteracting bosonic model and an interacting spin chain.

Rise and Fall, and Slow Rise Again, of Operator Entanglement under Dephasing

The operator space entanglement entropy, or simply "operator entanglement" (OE), is an indicator of the complexity of quantum operators and of their approximability by matrix product operators (MPOs). We study the OE of the density matrix of 1D many-body models undergoing dissipative evolution. It is expected that, after an initial linear growth reminiscent of unitary quench dynamics, the OE should be suppressed by dissipative processes as the system evolves to a simple stationary state. Surprisingly, we find that this scenario breaks down for one of the most fundamental dissipative mechanisms: dephasing. Under dephasing, after the initial "rise and fall," the OE can rise again, increasing logarithmically at long times. Using a combination of MPO simulations for chains of infinite length and analytical arguments valid for strong dephasing, we demonstrate that this growth is inherent to a U(1) conservation law. We argue that in an XXZ spin model and a Bose-Hubbard model the OE grows universally as 1/4log_{2}t at long times and as 1/2log_{2}t for a Fermi-Hubbard model. We trace this behavior back to anomalous classical diffusion processes.

Impact of Dissipation on Universal Fluctuation Dynamics in Open Quantum Systems

Recent theoretical and experimental works have explored universal dynamics related to surface growth physics in isolated quantum systems. In this Letter, we theoretically elucidate that dissipation drastically alters universal particle-number-fluctuation dynamics associated with surface-roughness growth in one-dimensional free fermions and bosons. In a system under dephasing that causes loss of spatial coherence, we numerically find that a universality class of surface-roughness dynamics changes from the ballistic class to a class with the Edwards-Wilkinson scaling exponents and an unconventional scaling function. We provide the analytical derivation of the diffusion equation from the dephasing Lindblad equation via a renormalization-group technique and succeed in explaining the drastic change. Furthermore, we numerically find the same change of the universality class under a more nontrivial dissipation, i.e., symmetric incoherent hopping.

Exact Solutions of Interacting Dissipative Systems via Weak Symmetries

We demonstrate how the presence of continuous weak symmetry can be used to analytically diagonalize the Liouvillian of a class of Markovian dissipative systems with strong interactions or nonlinearity. This enables an exact description of the full dynamics and dissipative spectrum. Our method can be viewed as implementing an exact, sector-dependent mean-field decoupling, or alternatively, as a kind of quantum-to-classical mapping. We focus on two canonical examples: a nonlinear bosonic mode subject to incoherent loss and pumping, and an inhomogeneous quantum Ising model with arbitrary connectivity and local dissipation. In both cases, we calculate and analyze the full dissipation spectrum. Our method is applicable to a variety of other systems, and could provide a powerful new tool for the study of complex driven-dissipative quantum systems.

Constructing Integrable Lindblad Superoperators

We develop a new method for the construction of one-dimensional integrable Lindblad systems, which describe quantum many body models in contact with a Markovian environment. We find several new models with interesting features, such as annihilation-diffusion processes, a mixture of coherent and classical particle propagation, and a rectified steady state current. We also find new ways to represent known classical integrable stochastic equations by integrable Lindblad operators. Our method can be extended to various other situations and it establishes a structured approach to the study of solvable open quantum systems.

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.

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.

Exact Liouvillian Spectrum of a One-Dimensional Dissipative Hubbard Model

A one-dimensional dissipative Hubbard model with two-body loss is shown to be exactly solvable. We obtain an exact eigenspectrum of a Liouvillian superoperator by employing a non-Hermitian extension of the Bethe-ansatz method. We find steady states, the Liouvillian gap, and an exceptional point that is accompanied by the divergence of the correlation length. A dissipative version of spin-charge separation induced by the quantum Zeno effect is also demonstrated. Our result presents a new class of exactly solvable Liouvillians of open quantum many-body systems, which can be tested with ultracold atoms subject to inelastic collisions.

Vaporization Dynamics of a Dissipative Quantum Liquid

We investigate the stability of a Luttinger liquid, upon suddenly coupling it to a dissipative environment. Within the Lindblad equation, the environment couples to local currents and heats the quantum liquid up to infinite temperatures. The single particle density matrix reveals the fractionalization of fermionic excitations in the spatial correlations by retaining the initial noninteger power law exponents, accompanied by an exponential decay in time with an interaction dependent rate. The spectrum of the time evolved density matrix is gapped, which collapses gradually as -ln(t). The von Neumann entropy crosses over from the early time -tln(t) behavior to ln(t) growth for late times. The early time dynamics is captured numerically by performing simulations on spinless interacting fermions, using several numerically exact methods. Our results could be tested experimentally in bosonic Luttinger liquids.

Integrability of one-dimensional Lindbladians from operator-space fragmentation

We introduce families of one-dimensional Lindblad equations describing open many-particle quantum systems that are exactly solvable in the following sense: (i) The space of operators splits into exponentially many (in system size) subspaces that are left invariant under the dissipative evolution; (ii) the time evolution of the density matrix on each invariant subspace is described by an integrable Hamiltonian. The prototypical example is the quantum version of the asymmetric simple exclusion process (ASEP) which we analyze in some detail. We show that in each invariant subspace the dynamics is described in terms of an integrable spin-1/2 XXZ Heisenberg chain with either open or twisted boundary conditions. We further demonstrate that Lindbladians featuring integrable operator-space fragmentation can be found in spin chains with arbitrary local physical dimensions.

Exact Nonequilibrium Steady State of Open XXZ/XYZ Spin-1/2 Chain with Dirichlet Boundary Conditions

We investigate a dissipatively driven XYZ spin-1/2 chain in the Zeno limit of strong dissipation, described by the Lindblad master equation. The nonequilibrium steady state is expressed in terms of a matrix product ansatz using novel site-dependent Lax operators. The components of Lax operators satisfy a simple set of linear recurrence equations that generalize the defining algebraic relations of the quantum group U_{q}(sl_{2}). We reveal connection between the nonequilibrium steady state of the nonunitary dynamics and the respective integrable model with edge magnetic fields, described by coherent unitary dynamics.

Dissipation-Induced Luttinger Liquid Correlations in a One-Dimensional Fermi Gas

We study a one-dimensional Fermi gas in the presence of dissipative coupling to environment through the Lindblad equation. The dissipation involves energy exchange with the environment and favours the relaxation of electrons to excitations. After switching on the dissipation, the system approaches a steady state, which is described by a generalized Gibbs ensemble. The fermionic single particle density matrix resembles deceivingly to that in a hermitian interaction quench. It decays inversely with the distance for short times due to the fermionic correlations in the initial state, which changes into a noninteger power law decay for late times, representing dissipation-induced Luttinger liquid behavior. However, the crossover between the two regions occurs due to dissipation-induced damping, and is unrelated to the propagation of excitations. The velocity of information spreading is set by the dissipative coupling, and differs significantly from the original sound velocity. The thermodynamic entropy grows as -t ln t initially, and saturates to an extensive value. Our results can be tested experimentally in one-dimensional Dirac systems.

Inhomogeneous matrix product ansatz and exact steady states of boundary-driven spin chains at large dissipation

We find novel site-dependent Lax operators in terms of which we demonstrate exact solvability of a dissipatively driven XYZ spin-1/2 chain in the Zeno limit of strong dissipation, with jump operators polarizing the boundary spins in arbitrary directions. We write the corresponding nonequilibrium steady state using an inhomogeneous matrix product ansatz, where the constituent matrices satisfy a simple set of linear recurrence relations. Although these matrices can be embedded into an infinite-dimensional auxiliary space, we have verified that they cannot be simultaneously put into a tridiagonal form, not even in the case of axially symmetric (XXZ) bulk interactions and general nonlongitudinal boundary dissipation. We expect our results to have further fundamental applications for the construction of nonlocal integrals of motion for the open XYZ model with arbitrary boundary fields, or the eight-vertex model.

Universal Signature from Integrability to Chaos in Dissipative Open Quantum Systems

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.

Non-Hermitian Many-Body Localization

Many-body localization is shown to suppress the imaginary parts of complex eigenenergies for general non-Hermitian Hamiltonians having time-reversal symmetry. We demonstrate that a real-complex transition, which we conjecture occurs upon many-body localization, profoundly affects the dynamical stability of non-Hermitian interacting systems with asymmetric hopping that respects time-reversal symmetry. Moreover, the real-complex transition is shown to be absent in non-Hermitian many-body systems with gain and/or loss that breaks time-reversal symmetry, even though the many-body localization transition still persists.

Open Quantum Symmetric Simple Exclusion Process

We present the solution to a model of fermions hopping between neighboring sites on a line with random Brownian amplitudes and open boundary conditions driving the system out of equilibrium. The average dynamics reduces to that of the symmetric simple exclusion process. However, the full distribution encodes for a richer behavior, entailing fluctuating quantum coherences which survive in the steady limit. We determine exactly the steady statistical distribution of the system state. We show that the out-of-equilibrium quantum coherence fluctuations satisfy a large-deviation principle, and we present a method to recursively compute exactly the large-deviation function. As a by-product, our approach gives a solution of the classical symmetric simple exclusion process based on fermion technology. Our results open the route towards the extension of the macroscopic fluctuation theory to many-body quantum systems.

Orthogonality Catastrophe in Dissipative Quantum Many-Body Systems

We present an analog of the phenomenon of orthogonality catastrophe in quantum many-body systems subject to a local dissipative impurity. We show that the fidelity F(t), giving a measure for distance of the time-evolved state from the initial one, displays a universal scaling form F(t)∝t^{θ}e^{-γt}, when the system supports long-range correlations, in a fashion reminiscent of traditional instances of orthogonality catastrophe in condensed matter. An exponential falloff at rate γ signals the onset of environmental decoherence, which is critically slowed down by the additional algebraic contribution to the fidelity. This picture is derived within a second-order cumulant expansion suited for Liouvillian dynamics, and substantiated for the one-dimensional transverse field quantum Ising model subject to a local dephasing jump operator, as well as for XY and XX quantum spin chains, and for the two-dimensional Bose gas deep in the superfluid phase with local particle heating. Our results hint that local sources of dissipation can be used to inspect real-time correlations and to induce a delay of decoherence in open quantum many-body systems.

Integrable Quantum Dynamics of Open Collective Spin Models

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.

Perfect Diode in Quantum Spin Chains

We study the rectification of the spin current in XXZ chains segmented in two parts, each with a different anisotropy parameter. Using exact diagonalization and a matrix product state algorithm, we find that a large rectification (of the order of 10^{4}) is attainable even using a short chain of N=8 spins, when one-half of the chain is gapless while the other has a large enough anisotropy. We present evidence of diffusive transport when the current is driven in one direction and of a transition to an insulating behavior of the system when driven in the opposite direction, leading to a perfect diode in the thermodynamic limit. The above results are explained in terms of matching of the spectrum of magnon excitations between the two halves of the chain.

Noisy Spins and the Richardson-Gaudin Model

We study a system of spins (qubits) coupled to a common noisy environment, each precessing at its own frequency. The correlated noise experienced by the spins implies long-lived correlations that relax only due to the differing frequencies. We use a mapping to a non-Hermitian integrable Richardson-Gaudin model to find the exact spectrum of the quantum master equation in the high-temperature limit and, hence, determine the decay rate. Our solution can be used to evaluate the effect of inhomogeneous splittings on a system of qubits coupled to a common bath.