Exact nonequilibrium steady state of a strongly driven open XXZ chain

Stable quantum-correlated many-body states through engineered dissipation

Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-range quantum correlations and a ground-state fidelity of 0.86 for 18 qubits at the critical point. In two dimensions, we found mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we explored transport in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors.

Stochastic Thermodynamics of a Finite Quantum System Coupled to Two Heat Baths

We consider a situation where an N-level system (NLS) is coupled successively to two heat baths with different temperatures without being necessarily thermalized and approaches a steady state. For this situation we apply a general Jarzynski-type equation and conclude that heat and entropy is flowing from the hot bath to the cold one. The Clausius relation between increase of entropy and transfer of heat divided by a suitable temperature assumes the form of two inequalities. Our approach is illustrated by an analytical example. For the linear regime, i.e., for small temperature differences between the two heat baths, we derive an expression for the heat conduction coefficient.

Establishing non-zero energy currents with the one-way street phenomenon and other symmetry properties in boundary driven spin systems

The aim of this paper is the investigation of properties of the energy current of usual (frequently found in the literature) boundary driven spin systems. For inhomogeneous Heisenberg (XXZandXXX) spin chains, we numerically compute the steady state, in the absence of an external magnetic field, and confirm the previously shown occurrence of the one-way street phenomenon, precisely, there is a nonzero energy current that preserves its magnitude and direction as we invert the baths at the edges, an effect stronger than the perfect rectification (which means current in a direction and zero current as we invert the baths). The consideration of several different polarizations at the edges reestablishes that it is a ubiquitous phenomenon. And, even for these inhomogeneous versions of spin chains, we also establish, by analytical methods, other symmetry properties of the energy current and confirm them by numerical computations.

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.

Transport and spectral properties of the XX+XXZ diode and stability to dephasing

We study the transport and spectral property of a segmented diode formed by an XX+XXZ spin chain. This system has been shown to become an ideal rectifier for spin current for large enough anisotropy. Here we show numerical evidence that the system in reverse bias has signatures pointing toward the existence of three different transport regimes depending on the value of the anisotropy: ballistic, diffusive, and insulating. In forward bias we observe two regimes, ballistic and diffusive. The system in forward and reverse bias shows significantly different spectral properties, with distribution of rapidities converging toward different functions. In the presence of dephasing the system becomes diffusive, rectification is significantly reduced, the relaxation gap increases, and the spectral properties in forward and reverse bias tend to converge. For large dephasing the relaxation gap decreases again as a result of quantum Zeno physics.

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.

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.

Giant Magnetoresistance in Boundary-Driven Spin Chains

In solid state physics, giant magnetoresistance is the large change in electrical resistance due to an external magnetic field. Here we show that giant magnetoresistance is possible in a spin chain composed of weakly interacting layers of strongly coupled spins. This is found for all system sizes even down to a minimal system of four spins. The mechanism driving the effect is a mismatch in the energy spectrum resulting in spin excitations being reflected at the boundaries between layers. This mismatch, and thus the current, can be controlled by external magnetic fields resulting in giant magnetoresistance. A simple rule for determining the behavior of the spin transport under the influence of a magnetic field is presented based on the energy levels of the strongly coupled spins.

Dynamics of large deviations in the hydrodynamic limit: Noninteracting systems

We study the dynamics of the statistics of the energy transferred across a point along a quantum chain which is prepared in the inhomogeneous initial state obtained by joining two identical semi-infinite parts thermalized at two different temperatures. In particular, we consider the transverse field Ising and harmonic chains as prototypical models of noninteracting fermionic and bosonic excitations, respectively. Within the so-called hydrodynamic limit of large space-time scales we first discuss the mean values of the energy density and current, and then, aiming at the statistics of fluctuations, we calculate exactly the scaled cumulant generating function of the transferred energy. From the latter, the evolution of the associated large deviation function is obtained. A natural interpretation of our results is provided in terms of a semiclassical picture of quasiparticles moving ballistically along classical trajectories. Similarities and differences between the transferred energy scaled cumulant and the large deviation functions in the cases of noninteracting fermions and bosons are discussed.

Giant Spin Current Rectification Due to the Interplay of Negative Differential Conductance and a Non-Uniform Magnetic Field

In XXZ chains with large enough interactions, spin transport can be significantly suppressed when the bias of the dissipative driving becomes large enough. This phenomenon of negative differential conductance is caused by the formation of two oppositely polarized ferromagnetic domains at the edges of the chain. Here, we show that this many-body effect, combined with a non-uniform magnetic field, can allow for a high degree of control of the spin current. In particular, by studying all of the possible shapes of local magnetic fields potentials, we find that a configuration in which the magnetic field points up for half of the chain and down for the other half, can result in giant spin-current rectification, for example, up to 108 for a system with only 8 spins. Our results show clear indications that the rectification can increase with the system size.

Using Matrix-Product States for Open Quantum Many-Body Systems: Efficient Algorithms for Markovian and Non-Markovian Time-Evolution

This paper presents an efficient algorithm for the time evolution of open quantum many-body systems using matrix-product states (MPS) proposing a convenient structure of the MPS-architecture, which exploits the initial state of system and reservoir. By doing so, numerically expensive re-ordering protocols are circumvented. It is applicable to systems with a Markovian type of interaction, where only the present state of the reservoir needs to be taken into account. Its adaption to a non-Markovian type of interaction between the many-body system and the reservoir is demonstrated, where the information backflow from the reservoir needs to be included in the computation. Also, the derivation of the basis in the quantum stochastic Schrödinger picture is shown. As a paradigmatic model, the Heisenberg spin chain with nearest-neighbor interaction is used. It is demonstrated that the algorithm allows for the access of large systems sizes. As an example for a non-Markovian type of interaction, the generation of highly unusual steady states in the many-body system with coherent feedback control is demonstrated for a chain length of N=30.

Entanglement statistics in Markovian open quantum systems: A matter of mutation and selection

Controlling dynamical fluctuations in open quantum systems is essential both for our comprehension of quantum nonequilibrium behavior and for its possible application in near-term quantum technologies. However, understanding these fluctuations is extremely challenging due, to a large extent, to a lack of efficient important sampling methods for quantum systems. Here, we devise a unified framework-based on population-dynamics methods-for the evaluation of the full probability distribution of generic time-integrated observables in Markovian quantum jump processes. These include quantities carrying information about genuine quantum features, such as quantum superposition or entanglement, not accessible with existing numerical techniques. The algorithm we propose provides dynamical free-energy and entropy functionals which, akin to their equilibrium counterpart, permit one to unveil intriguing phase-transition behavior in quantum trajectories. We discuss some applications and further disclose coexistence and hysteresis, between a highly entangled phase and a low entangled one, in large fluctuations of a strongly interacting few-body system.

Uncoupled Majorana fermions in open quantum systems: on the efficient simulation of non-equilibrium stationary states of quadratic Fermi models

A decomposition of the non-equilibrium stationary state of a quadratic Fermi system influenced by linear baths is obtained and used to establish a simulation protocol in terms of tensor states. The scheme is then applied to examine the occurrence of uncoupled Majorana fermions in Kitaev chains subject to baths on the ends. The resulting phase diagram is compared against the topological characterization of the equilibrium chain and the protocol efficiency is studied with respect to this model.

Thermalization in open many-body systems based on eigenstate thermalization hypothesis

We investigate steady states of macroscopic quantum systems under dissipation not obeying the detailed balance condition. We argue that the Gibbs state at an effective temperature gives a good description of the steady state provided that the system Hamiltonian obeys the eigenstate thermalization hypothesis (ETH) and the perturbation theory in the weak system-environment coupling is valid in the thermodynamic limit. We derive a criterion to guarantee the validity of the perturbation theory, which is satisfied in the thermodynamic limit for sufficiently weak dissipation when the Liouvillian is gapped for bulk-dissipated systems, while the perturbation theory breaks down in boundary-dissipated chaotic systems due to the presence of diffusive transports. We numerically confirm these theoretical predictions. This paper suggests a connection between steady states of macroscopic open quantum systems and the ETH.

Large deviations and fluctuation theorem for the quantum heat current in the spin-boson model

We study the heat current flowing between two baths consisting of harmonic oscillators interacting with a qubit through a spin-boson coupling. An explicit expression for the generating function of the total heat flowing between the right and left baths is derived by evaluating the corresponding Feynman-Vernon path integral by performing the noninteracting blip approximation (NIBA). We recover the known expression, obtained by using the polaron transform. This generating function satisfies the Gallavotti-Cohen fluctuation theorem, both before and after performing the NIBA. We also verify that the heat conductance is proportional to the variance of the heat current, retrieving the well-known fluctuation dissipation relation. Finally, we present numerical results for the heat current.

Quantum diffusion in spin chains with phase space methods

Connecting short-time microscopic dynamics with long-time hydrodynamics in strongly correlated quantum systems is one of the outstanding questions. In particular, it is hard to determine various hydrodynamic coefficients such as the diffusion constant or viscosity starting from a microscopic model: exact quantum simulations are limited to either small system sizes or to short times, which are insufficient to reach asymptotic behavior and so various approximations must be applied. We show that these difficulties, at least for particular models, can be circumvented by using the cluster truncated Wigner approximation (CTWA), which maps quantum Hamiltonian dynamics into classical Hamiltonian dynamics in auxiliary high-dimensional phase space. We apply CTWA to XXZ next-nearest-neighbor spin-1/2 chains and XY spin ladders, and find behavior consisting of short-time spin relaxation which gradually crosses over to emergent diffusive behavior at long times. For a random initial state, we show that CTWA correctly reproduces the whole spin spectral function. Necessary in this construction is sampling from properly fluctuating initial conditions: the Dirac mean-field (variational) ansatz, which neglects such fluctuations, leads to incorrect predictions.

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.

Energy current manipulation and reversal of rectification in graded XXZ spin chains

This work is devoted to the investigation of nontrivial transport properties in many-body quantum systems. Precisely, we study transport in the steady state of spin-1/2 Heisenberg XXZ chains, driven out of equilibrium by two magnetic baths at their end points. We take graded versions of the model, i.e. asymmetric chains in which some structure gradually changes in space. We investigate how we can manipulate and control the energy and spin currents of such chains by tuning external and/or inner parameters. In particular, we describe the occurrence of energy current rectification and its reversal due to the application of external magnetic fields. We show that, after carefully chosen inner parameters for the system, by turning on an external magnetic field we can find spin and energy currents propagating in different directions. More interestingly, we may find cases in which rectifications of energy and spin currents occur in opposite directions, i.e. if the energy current is larger when flowing from left to right side, then the spin current is larger if it flows from right to left side. We still describe situations with inversion of the energy current direction as we increase the system asymmetry. We stress that our work aims the development of theoretical knowledge as well as the stimulation of future experimental applications.

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.

Quantitative Impact of Integrals of Motion on the Eigenstate Thermalization Hypothesis

Even though the eigenstate thermalization hypothesis (ETH) may be introduced as an extension of the random matrix theory, physical Hamiltonians and observables differ from random operators. One of the challenges is to embed local integrals of motion (LIOMs) within the ETH. Here we make steps towards a unified treatment of the ETH in integrable and nonintegrable models with translational invariance. Specifically, we focus on the impact of LIOMs on the fluctuations and structure of the diagonal matrix elements of local observables. We first show that nonvanishing fluctuations entail the presence of LIOMs. Then we introduce a generic protocol to construct observables, subtracted by their projections on LIOMs as well as products of LIOMs. The protocol systematically reduces fluctuations and/or the structure of the diagonal matrix elements. We verify our arguments by numerical results for integrable and nonintegrable models.

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.

Critical Behavior of the Quantum Contact Process in One Dimension

The contact process is a paradigmatic classical stochastic system displaying critical behavior even in one dimension. It features a nonequilibrium phase transition into an absorbing state that has been widely investigated and shown to belong to the directed percolation universality class. When the same process is considered in a quantum setting, much less is known. So far, mainly semiclassical studies have been conducted and the nature of the transition in low dimensions is still a matter of debate. Also, from a numerical point of view, from which the system may look fairly simple-especially in one dimension-results are lacking. In particular, the presence of the absorbing state poses a substantial challenge, which appears to affect the reliability of algorithms targeting directly the steady state. Here we perform real-time numerical simulations of the open dynamics of the quantum contact process and shed light on the existence and on the nature of an absorbing state phase transition in one dimension. We find evidence for the transition being continuous and provide first estimates for the critical exponents. Beyond the conceptual interest, the simplicity of the quantum contact process makes it an ideal benchmark problem for scrutinizing numerical methods for open quantum nonequilibrium systems.

Thermalization with detailed-balanced two-site Lindblad dissipators

The use of two-site Lindblad dissipators to generate thermal states and study heat transport was raised to prominence by Prosen and Žnidarič [J. Stat. Mech. (2009) P020351742-546810.1088/1742-5468/2009/02/P02035]. Here we propose a variant of this method based on detailed balance of internal levels of the two-site Hamiltonian and characterize its performance. We study the thermalization profile in the chain, the effective temperatures achieved by different single- and two-site observables, and we also investigate the decay of two-time correlations. We find that at a large enough temperature, the steady state approaches closely a thermal state, with a relative error below 1% for the inverse temperature estimated from different observables.

Variational Neural-Network Ansatz for Steady States in Open Quantum Systems

We present a general variational approach to determine the steady state of open quantum lattice systems via a neural-network approach. The steady-state density matrix of the lattice system is constructed via a purified neural-network Ansatz in an extended Hilbert space with ancillary degrees of freedom. The variational minimization of cost functions associated to the master equation can be performed using a Markov chain Monte Carlo sampling. As a first application and proof of principle, we apply the method to the dissipative quantum transverse Ising model.

Neural-Network Approach to Dissipative Quantum Many-Body Dynamics

In experimentally realistic situations, quantum systems are never perfectly isolated and the coupling to their environment needs to be taken into account. Often, the effect of the environment can be well approximated by a Markovian master equation. However, solving this master equation for quantum many-body systems becomes exceedingly hard due to the high dimension of the Hilbert space. Here we present an approach to the effective simulation of the dynamics of open quantum many-body systems based on machine-learning techniques. We represent the mixed many-body quantum states with neural networks in the form of restricted Boltzmann machines and derive a variational Monte Carlo algorithm for their time evolution and stationary states. We document the accuracy of the approach with numerical examples for a dissipative spin lattice system.

Variational Quantum Monte Carlo Method with a Neural-Network Ansatz for Open Quantum Systems

The possibility to simulate the properties of many-body open quantum systems with a large number of degrees of freedom (d.o.f.) is the premise to the solution of several outstanding problems in quantum science and quantum information. The challenge posed by this task lies in the complexity of the density matrix increasing exponentially with the system size. Here, we develop a variational method to efficiently simulate the nonequilibrium steady state of Markovian open quantum systems based on variational Monte Carlo methods and on a neural network representation of the density matrix. Thanks to the stochastic reconfiguration scheme, the application of the variational principle is translated into the actual integration of the quantum master equation. We test the effectiveness of the method by modeling the two-dimensional dissipative XYZ spin model on a lattice.

Mixed-Order Symmetry-Breaking Quantum Phase Transition Far from Equilibrium

We study the current-carrying steady state of a transverse field Ising chain coupled to magnetic thermal reservoirs and obtain the nonequilibrium phase diagram as a function of the magnetization potential of the reservoirs. Upon increasing the magnetization bias we observe a discontinuous jump of the magnetic order parameter that coincides with a divergence of the correlation length. For steady states with a nonvanishing conductance, the entanglement entropy at zero temperature displays a bias dependent logarithmic correction that violates the area law and differs from the well-known equilibrium case. Our findings show that out-of-equilibrium conditions allow for novel critical phenomena not possible at equilibrium.

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.

Integrable Trotterization: Local Conservation Laws and Boundary Driving

We discuss a general procedure to construct an integrable real-time Trotterization of interacting lattice models. As an illustrative example, we consider a spin-1/2 chain, with continuous time dynamics described by the isotropic (XXX) Heisenberg Hamiltonian. For periodic boundary conditions, local conservation laws are derived from an inhomogeneous transfer matrix, and a boost operator is constructed. In the continuous time limit, these local charges reduce to the known integrals of motion of the Heisenberg chain. In a simple Kraus representation, we also examine the nonequilibrium setting, where our integrable cellular automaton is driven by stochastic processes at the boundaries. We show explicitly how an exact nonequilibrium steady-state density matrix can be written in terms of a staggered matrix product ansatz, and we propose quasilocal conservation laws for the model with periodic boundary conditions. This simple Trotterization scheme, in particular in the open system framework, could prove to be a useful tool for experimental simulations of the lattice models in terms of trapped ion and atom optics setups.

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.

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.

Convergence radius of perturbative Lindblad-driven nonequilibrium steady states

We address the problem of analyzing the radius of convergence of perturbative expansion of nonequilibrium steady states of Lindblad-driven spin chains. A simple formal approach is developed for systematically computing the perturbative expansion of small driven systems. We consider the paradigmatic model of an open XXZ spin-1/2 chain with boundary-supported ultralocal Lindblad dissipators and treat two different perturbative cases: (i) expansion in the system-bath coupling parameter and (ii) expansion in the driving (bias) parameter. In the first case (i) we find that the radius of convergence quickly shrinks with increasing the system size, while in the second case (ii) we find that the convergence radius is always larger than 1, and in particular it approaches 1 from above as we change the anisotropy from an easy-plane (XY) to an easy-axis (Ising) regime.

Charge and spin current statistics of the open Hubbard model with weak coupling to the environment

Based on generalization and extension of our previous work [Phys. Rev. Lett. 112, 067201 (2014)PRLTAO0031-900710.1103/PhysRevLett.112.067201] to multiple independent Markovian baths we will compute the charge and spin current statistics of the open Hubbard model with weak system-bath coupling up to next-to-leading order in the coupling parameter. Only the next-to-leading and higher orders depend on the Hubbard interaction parameter. The physical results are related to those for the XXZ model in the analogous setup implying a certain universality, which potentially holds in this class of nonequilibrium models.

Microscopic theory of a nonequilibrium open bosonic chain

Quantum master equations form an important tool in the description of transport problems in open quantum systems. However, they suffer from the difficulty that the shape of the Lindblad dissipator depends sensibly on the system Hamiltonian. Consequently, most of the work done in this field has focused on phenomenological dissipators which act locally on different parts of the system. In this paper we show how to construct Lindblad dissipators to model a one-dimensional bosonic tight-binding chain connected to two baths at the first and last site, kept at different temperatures and chemical potentials. We show that even though the bath coupling is local, the effective Lindblad dissipator stemming from this interaction is inherently nonlocal, affecting all normal modes of the system. We then use this formalism to study the current of particles and energy through the system and find that they have the structure of Landauer's formula, with the bath spectral density playing the role of the transfer integral. Finally, we consider infinitesimal temperature and chemical potential gradients and show that the currents satisfy Onsager's reciprocal relations, which is a consequence of the fact that the microscopic quantum dynamics obeys detailed balance.

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.

Nonequilibrium quantum chains under multisite Lindblad baths

We study a quantum XX chain coupled to two heat reservoirs that act on multiple sites and are kept at different temperatures and chemical potentials. The baths are described by Lindblad dissipators, which are constructed by direct coupling to the fermionic normal modes of the chain. Using a perturbative method, we are able to find analytical formulas for all steady-state properties of the system. We compute both the particle or magnetization current and the energy current, both of which are found to have the structure of Landauer's formula. We also obtain exact formulas for the Onsager coefficients. All properties are found to differ substantially from those of a single-site bath. In particular, we find a strong dependence on the intensity of the bath couplings. In the weak-coupling regime, we show that the Onsager reciprocal relations are satisfied.

Positive Tensor Network Approach for Simulating Open Quantum Many-Body Systems

Open quantum many-body systems play an important role in quantum optics and condensed matter physics, and capture phenomena like transport, the interplay between Hamiltonian and incoherent dynamics, and topological order generated by dissipation. We introduce a versatile and practical method to numerically simulate one-dimensional open quantum many-body dynamics using tensor networks. It is based on representing mixed quantum states in a locally purified form, which guarantees that positivity is preserved at all times. Moreover, the approximation error is controlled with respect to the trace norm. Hence, this scheme overcomes various obstacles of the known numerical open-system evolution schemes. To exemplify the functioning of the approach, we study both stationary states and transient dissipative behavior, for various open quantum systems ranging from few to many bodies.

Nonequilibrium steady states of ideal bosonic and fermionic quantum gases

We investigate nonequilibrium steady states of driven-dissipative ideal quantum gases of both bosons and fermions. We focus on systems of sharp particle number that are driven out of equilibrium either by the coupling to several heat baths of different temperature or by time-periodic driving in combination with the coupling to a heat bath. Within the framework of (Floquet-)Born-Markov theory, several analytical and numerical methods are described in detail. This includes a mean-field theory in terms of occupation numbers, an augmented mean-field theory taking into account also nontrivial two-particle correlations, and quantum-jump-type Monte Carlo simulations. For the case of the ideal Fermi gas, these methods are applied to simple lattice models and the possibility of achieving exotic states via bath engineering is pointed out. The largest part of this work is devoted to bosonic quantum gases and the phenomenon of Bose selection, a nonequilibrium generalization of Bose condensation, where multiple single-particle states are selected to acquire a large occupation [Phys. Rev. Lett. 111, 240405 (2013)]. In this context, among others, we provide a theory for transitions where the set of selected states changes, describe an efficient algorithm for finding the set of selected states, investigate beyond-mean-field effects, and identify the dominant mechanisms for heat transport in the Bose-selected state.

The thermodynamic cost of driving quantum systems by their boundaries

The laws of thermodynamics put limits to the efficiencies of thermal machines. Analogues of these laws are now established for quantum engines weakly and passively coupled to the environment providing a framework to find improvements to their performance. Systems whose interaction with the environment is actively controlled do not fall in that framework. Here we consider systems actively and locally coupled to the environment, evolving with a so-called boundary-driven Lindblad equation. Starting from a unitary description of the system plus the environment we simultaneously obtain the Lindblad equation and the appropriate expressions for heat, work and entropy-production of the system extending the framework for the analysis of new, and some already proposed, quantum heat engines. We illustrate our findings in spin 1/2 chains and explain why an XX chain coupled in this way to a single heat bath relaxes to thermodynamic-equilibrium while and XY chain does not. Additionally, we show that an XX chain coupled to a left and a right heat baths behaves as a quantum engine, a heater or refrigerator depending on the parameters, with efficiencies bounded by Carnot efficiencies.

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.

Identifying local and quasilocal conserved quantities in integrable systems

We outline a procedure for counting and identifying a complete set of local and quasilocal conserved operators in integrable lattice systems. The method yields a systematic generation of all independent, conserved quasilocal operators related to the time average of local operators with a support on up to M consecutive sites. As an example, we study the anisotropic Heisenberg spin-1/2 chain and show that the number of independent conserved operators grows linearly with M. In addition to the known local operators, there exist novel quasilocal conserved quantities in all the parity sectors. The existence of quasilocal conserved operators is shown also for the isotropic Heisenberg model. Implications for the anomalous relaxation of quenched systems are discussed as well.

Infinitely dimensional lax structure for the one-dimensional Hubbard model

We report a two-parametric irreducible infinitely dimensional representation of the Lax integrability condition for the Fermi Hubbard chain. In addition to being of fundamental interest, hinting at possible novel quantum symmetry of the model, our construction allows for an explicit representation of an exact steady state many-body density operator for a nonequilibrium boundary-driven Hubbard chain with arbitrary (asymmetric) particle source (sink) rates at the left (right) end of the chain and with arbitrary boundary values of chemical potentials.

Flux rectification in the quantum XXZ chain

Thermal rectification is the phenomenon by which the flux of heat depends on the direction of the flow. It has attracted much interest in recent years due to the possibility of devising thermal diodes. In this paper, we consider the rectification phenomenon in the quantum XXZ chain subject to an inhomogeneous field. The chain is driven out of equilibrium by the contact at its boundaries with two different reservoirs, leading to a constant flow of magnetization from one bath to the other. The nonunitary dynamics of this system, which is modeled by a Lindblad master equation, is treated exactly for small sizes and numerically for larger ones. The functional dependence of the rectification coefficient on the model parameters (anisotropy, field amplitude, and out of equilibrium driving strength) is investigated in full detail. Close to the XX point and at small inhomogeneity and low driving, we have found an explicit expression for the rectification coefficient that is valid at all system sizes. In particular, it shows that the phenomenon of rectification persists even in the thermodynamic limit. Finally, we prove that in the case of the XX chain, there is no rectification.

Breakdown of the generalized Gibbs ensemble for current-generating quenches

We establish a relation between two hallmarks of integrable systems: the relaxation towards the generalized Gibbs ensemble (GGE) and the dissipationless charge transport. We show that the former one is possible only if the so-called Mazur bound on the charge stiffness is saturated by local conserved quantities. As an example we show how a non-GGE steady state with a current can be generated in the one-dimensional model of interacting spinless fermions with a flux quench. Moreover, an extended GGE involving the quasilocal conserved quantities can be formulated for this case.

Exactly solvable counting statistics in open weakly coupled interacting spin systems

We study the full counting statistics for interacting quantum many-body spin systems weakly coupled to the environment. In the leading order in the system-bath coupling, we derive exact spin current statistics for a large class of parity symmetric spin-1/2 systems driven by a pair of Markovian baths with local coupling operators. Interestingly, in this class of systems the leading-order current statistics are universal and do not depend on details of the Hamiltonian. Furthermore, in the specific case of a symmetrically boundary driven anisotropic Heisenberg (XXZ) spin-1/2 chain, we explicitly derive the third-order nonlinear corrections to the current statistics.

Exact nonequilibrium steady state of an open Hubbard chain

We discuss the current carrying nonequilibrium steady state of an open fermionic Hubbard chain that is strongly driven by Markovian incoherent processes localized at the chain ends. An explicit form of an exact many-body density operator for any value of the coupling parameter is presented. The structure of a matrix product form of the solution is encoded in terms of a novel diagrammatic technique that should allow for generalization to other integrable nonequilibrium models.

Lower bounds on high-temperature diffusion constants from quadratically extensive almost-conserved operators

We prove a general theorem which provides a strict lower bound on high-temperature Green-Kubo diffusion constants in locally interacting quantum lattice systems, under the assumption of existence of a quadratically extensive almost conserved quantity: an operator whose commutator with the lattice Hamiltonian is localized on the boundary sites only. We explicitly demonstrate and compute such a bound in two important models in one dimension: in the (isotropic) Heisenberg spin 1/2 chain and in the fermionic Hubbard chain.