Solvable quantum nonequilibrium model exhibiting a phase transition and a matrix product representation

Enhancing violations of Leggett-Garg inequalities in nonequilibrium correlated many-body systems by interactions and decoherence

We identify different schemes to enhance the violation of Leggett-Garg inequalities in open many-body systems. Considering a nonequilibrium archetypical setup of quantum transport, we show that particle interactions control the direction and amplitude of maximal violation, and that in the strongly-interacting and strongly-driven regime bulk dephasing enhances the violation. Through an analytical study of a minimal model we unravel the basic ingredients to explain this decoherence-enhanced quantumness, illustrating that such an effect emerges in a wide variety of systems.

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.

Uhlmann curvature in dissipative phase transitions

A novel approach based on the Uhlmann curvature is introduced for the investigation of non-equilibrium steady-state quantum phase transitions (NESS-QPTs). Equilibrium phase transitions fall invariably into two markedly non-overlapping categories: classical phase transitions and quantum phase transitions. NESS-QPTs offer a unique arena where such a distinction fades off. We propose a method to reveal and quantitatively assess the quantum character of such critical phenomena. We apply this tool to a paradigmatic class of lattice fermion systems with local reservoirs, characterised by Gaussian non-equilibrium steady states. The relations between the behaviour of the Uhlmann curvature, the divergence of the correlation length, the character of the criticality and the dissipative gap are demonstrated. We argue that this tool can shade light upon the nature of non equilibrium steady state criticality in particular with regard to the role played by quantum vs classical fluctuations.

Solvable Family of Driven-Dissipative Many-Body Systems

Exactly solvable models have played an important role in establishing the sophisticated modern understanding of equilibrium many-body physics. Conversely, the relative scarcity of solutions for nonequilibrium models greatly limits our understanding of systems away from thermal equilibrium. We study a family of nonequilibrium models, some of which can be viewed as dissipative analogues of the transverse-field Ising model, in that an effectively classical Hamiltonian is frustrated by dissipative processes that drive the system toward states that do not commute with the Hamiltonian. Surprisingly, a broad and experimentally relevant subset of these models can be solved efficiently. We leverage these solutions to compute the effects of decoherence on a canonical trapped-ion-based quantum computation architecture, and to prove a no-go theorem on steady-state phase transitions in a many-body model that can be realized naturally with Rydberg atoms or trapped ions.

Fluctuating hydrodynamics, current fluctuations, and hyperuniformity in boundary-driven open quantum chains

We consider a class of either fermionic or bosonic noninteracting open quantum chains driven by dissipative interactions at the boundaries and study the interplay of coherent transport and dissipative processes, such as bulk dephasing and diffusion. Starting from the microscopic formulation, we show that the dynamics on large scales can be described in terms of fluctuating hydrodynamics. This is an important simplification as it allows us to apply the methods of macroscopic fluctuation theory to compute the large deviation (LD) statistics of time-integrated currents. In particular, this permits us to show that fermionic open chains display a third-order dynamical phase transition in LD functions. We show that this transition is manifested in a singular change in the structure of trajectories: while typical trajectories are diffusive, rare trajectories associated with atypical currents are ballistic and hyperuniform in their spatial structure. We confirm these results by numerically simulating ensembles of rare trajectories via the cloning method, and by exact numerical diagonalization of the microscopic quantum generator.

Solution of the Lindblad equation for spin helix states

Using Lindblad dynamics we study quantum spin systems with dissipative boundary dynamics that generate a stationary nonequilibrium state with a nonvanishing spin current that is locally conserved except at the boundaries. We demonstrate that with suitably chosen boundary target states one can solve the many-body Lindblad equation exactly in any dimension. As solution we obtain pure states at any finite value of the dissipation strength and any system size. They are characterized by a helical stationary magnetization profile and a ballistic spin current which is independent of system size, even when the quantum spin system is not integrable. These results are derived in explicit form for the one-dimensional spin-1/2 Heisenberg chain and its higher-spin generalizations, which include the integrable spin-1 Zamolodchikov-Fateev model and the biquadratic Heisenberg chain.

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

We construct an exact map between a tight-binding model on any bipartite lattice in the presence of dephasing noise and a Hubbard model with imaginary interaction strength. In one dimension, the exact many-body Liouvillian spectrum can be obtained by application of the Bethe ansatz method. We find that both the nonequilibrium steady state and the leading decay modes describing the relaxation at late times are related to the η-pairing symmetry of the Hubbard model. We show that there is a remarkable relation between the time evolution of an arbitrary k-point correlation function in the dissipative system and k-particle states of the corresponding Hubbard model.

Relaxation times of dissipative many-body quantum systems

We study relaxation times, also called mixing times, of quantum many-body systems described by a Lindblad master equation. We in particular study the scaling of the spectral gap with the system length, the so-called dynamical exponent, identifying a number of transitions in the scaling. For systems with bulk dissipation we generically observe different scaling for small and for strong dissipation strength, with a critical transition strength going to zero in the thermodynamic limit. We also study a related phase transition in the largest decay mode. For systems with only boundary dissipation we show a generic bound that the gap cannot be larger than ∼1/L. In integrable systems with boundary dissipation one typically observes scaling of ∼1/L(3), while in chaotic ones one can have faster relaxation with the gap scaling as ∼1/L and thus saturating the generic bound. We also observe transition from exponential to algebraic gap in systems with localized modes.

Large-deviation statistics of a diffusive quantum spin chain and the additivity principle

Using the large-deviation formalism, we study the statistics of current fluctuations in a diffusive nonequilibrium quantum spin chain. The boundary-driven XX chain with dephasing consists of a coherent bulk hopping and a local dissipative dephasing. We analytically calculate the exact expression for the second current moment in a system of any length and then numerically demonstrate that in the thermodynamic limit, higher-order cumulants and the large-deviation function can be calculated using the additivity principle or macroscopic hydrodynamic theory. This shows that the additivity principle can also hold in systems that are not purely stochastic, and can in particular be valid in quantum systems. We also show that in large systems, the current fluctuations are the same as in the classical symmetric simple exclusion process.

Quantum information-geometry of dissipative quantum phase transitions

A general framework for analyzing the recently discovered phase transitions in the steady state of dissipation-driven open quantum systems is still lacking. To fill this gap, we extend the so-called fidelity approach to quantum phase transitions to open systems whose steady state is a Gaussian fermionic state. We endow the manifold of correlation matrices of steady states with a metric tensor g measuring the distinguishability distance between solutions corresponding to a different set of control parameters. The phase diagram can then be mapped out in terms of the scaling behavior of g and connections with the Liouvillean gap and the model correlation functions unveiled. We argue that the fidelity approach, thanks to its differential-geometric and information-theoretic nature, provides insights into dissipative quantum critical phenomena as well as a general and powerful strategy to explore them.

Exact large-deviation statistics for a nonequilibrium quantum spin chain

We consider a one-dimensional XX spin chain in a nonequilibrium setting with a Lindblad-type boundary driving. By calculating large-deviation rate function in the thermodynamic limit, a generalization of free energy to a nonequilibrium setting, we obtain a complete distribution of current, including closed expressions for lower-order cumulants. We also identify two phase-transition-like behaviors in either the thermodynamic limit, at which the current probability distribution becomes discontinuous, or at maximal driving, when the range of possible current values changes discontinuously. In the thermodynamic limit the current has a finite upper and lower bound. We also explicitly confirm nonequilibrium fluctuation relation and show that the current distribution is the same under mapping of the coupling strength Γ→1/Γ.

Eigenvalue statistics as an indicator of integrability of nonequilibrium density operators

We propose to quantify the complexity of nonequilibrium steady state density operators, as well as of long-lived Liouvillian decay modes, in terms of the level spacing distribution of their spectra. Based on extensive numerical studies in a variety of models, some solvable and some unsolved, we conjecture that the integrability of density operators (e.g., the existence of an algebraic procedure for their construction in finitely many steps) is signaled by a Poissonian level statistics, whereas in the generic nonintegrable cases one finds level statistics of a Gaussian unitary ensemble of random matrices. Eigenvalue statistics can therefore be used as an efficient tool to identify integrable quantum nonequilibrium systems.

Heat transport through lattices of quantum harmonic oscillators in arbitrary dimensions

In d-dimensional lattices of coupled quantum harmonic oscillators, we analyze the heat current caused by two thermal baths of different temperatures, which are coupled to opposite ends of the lattice, with a focus on the validity of Fourier's law of heat conduction. We provide analytical solutions of the heat current through the quantum system in the nonequilibrium steady state using the rotating-wave approximation and bath interactions described by a master equation of Lindblad form. The influence of local dephasing in the transition of ballistic to diffusive transport is investigated.

Quantum transport efficiency and Fourier's law

We analyze the steady-state energy transfer in a chain of coupled two-level systems connecting two thermal reservoirs. Through an analytic treatment we find that the energy current is independent of the system size, hence violating Fourier's law of heat conduction. The classical diffusive behavior in Fourier's law of heat conduction can be recovered by introducing decoherence to the quantum systems constituting the chain. We relate these results to recent discussions of energy transport in biological light-harvesting systems, and discuss the role of quantum coherence and entanglement.

Transport properties of a boundary-driven one-dimensional gas of spinless fermions

We analytically study a system of spinless fermions driven at the boundary with an oscillating chemical potential. Various transport regimes can be observed: At zero driving frequency the particle current through the system is independent of the system's length; at the phase-transition frequency, being equal to the bandwidth, the current decays as ~n(-α) with the chain length n, α being either 2 or 3; below the transition the scaling of the current is ~n(-1/2), indicating anomalous transport, while it is exponentially small ~exp(-n/2ξ) above the transition. Therefore, by a simple change of frequency of the a.c. driving of the system one can vary transport from ballistic, anomalous to insulating.

Exact nonequilibrium steady state of a strongly driven open XXZ chain

An exact and explicit ladder-tensor-network ansatz is presented for the nonequilibrium steady state of an anisotropic Heisenberg XXZ spin-1/2 chain which is driven far from equilibrium with a pair of Lindblad operators acting on the edges of the chain only. We show that the steady-state density operator of a finite system of size n is-apart from a normalization constant-a polynomial of degree 2n - 2 in the coupling constant. Efficient computation of physical observables is facilitated in terms of a transfer operator reminiscent of a classical Markov process. In the isotropic case we find cosine spin profiles, 1/n(2) scaling of the spin current, and long-range correlations in the steady state. This is a fully nonperturbative extension of a recent result [Phys. Rev. Lett. 106, 217206 (2011)].

Nonequilibrium phase transition in a periodically driven XY spin chain

We present a general formulation of Floquet states of periodically time-dependent open Markovian quasifree fermionic many-body systems in terms of a discrete Lyapunov equation. Illustrating the technique, we analyze periodically kicked XY spin-½ chain which is coupled to a pair of Lindblad reservoirs at its ends. A complex phase diagram is reported with reentrant phases of long range and exponentially decaying spin-spin correlations as some of the system's parameters are varied. The structure of phase diagram is reproduced in terms of counting nontrivial stationary points of Floquet quasiparticle dispersion relation.

Spin transport in a one-dimensional anisotropic Heisenberg model

We analytically and numerically study spin transport in a one-dimensional Heisenberg model in linear-response regime at infinite temperature. It is shown that as the anisotropy parameter Δ is varied spin transport changes from ballistic for Δ<1 to anomalous at the isotropic point Δ=1, to diffusive for finite Δ>1, ending up as a perfect isolator in the Ising limit of infinite Δ. Using perturbation theory for large Δ a quantitative prediction is made for the dependence of diffusion constant on Δ.

Open XXZ spin chain: nonequilibrium steady state and a strict bound on ballistic transport

An explicit matrix product ansatz is presented, in the first two orders in the (weak) coupling parameter, for the nonequilibrium steady state of the homogeneous, nearest neighbor Heisenberg XXZ spin 1/2 chain driven by Lindblad operators which act only at the edges of the chain. The first order of the density operator becomes, in the thermodynamic limit, an exact pseudolocal conservation law and yields-via the Mazur inequality-a rigorous lower bound on the high-temperature spin Drude weight. Such a Mazur bound is a nonvanishing fractal function of the anisotropy parameter Δ for |Δ|<1.