Exterior integrability : Yang-Baxter form of non-equilibrium steady-state density operator

Prosen, T., Enej Ilievski, Popkov, V.

Published Jul 30, 2013• 27 citations

Abstract

A new type of quantum transfer matrix, arising as a Cholesky factor for the steady-state density matrix of a dissipative Markovian process associated with the boundary-driven Lindblad equation for the isotropic spin-1/2 Heisenberg (XXX) chain, is presented. The transfer matrix forms a commuting family of non-Hermitian operators depending on the spectral parameter, which is essentially the strength of dissipative coupling at the boundaries. The intertwining of the corresponding Lax and monodromy matrices is performed by an infinitely dimensional Yang–Baxter R-matrix, which we construct explicitly and is essentially different from the standard 4 × 4 XXX R-matrix. We also discuss a possibility to construct Bethe ansatz for the spectrum and eigenstates of the non-equilibrium steady-state density operator. Furthermore, we indicate the existence of a deformed R-matrix in the infinite dimensional auxiliary space for the anisotropic XXZ spin-1/2 chain, which in general provides a sequence of new, possibly quasi-local, conserved quantities of the bulk XXZ dynamics.

Related papers

Exterior integrability: Yang–Baxter form of non-equilibrium steady-state density operator

Tomaž Prosen, Enej Ilievski, Vladislav Popkov

A new type of quantum transfer matrix, arising as a Cholesky factor for the steady-state density matrix of a dissipative Markovian process associated with the boundary-driven Lindblad equation for the isotropic spin-1/2 Heisenberg ( XXX ) chain, is presented. The transfer matrix forms a commuting family of non-Hermitian operators depending on the spectral parameter, which is essentially the strength of dissipative coupling at the boundaries. The intertwining of the corresponding Lax and monodromy matrices is performed by an infinitely dimensional Yang–Baxter R -matrix, which we construct explicitly and is essentially different from the standard 4 × 4 XXX R -matrix. We also discuss a possibility to construct Bethe ansatz for the spectrum and eigenstates of the non-equilibrium steady-state density operator. Furthermore, we indicate the existence of a deformed R -matrix in the infinite dimensional auxiliary space for the anisotropic XXZ spin-1/2 chain, which in general provides a sequence of new, possibly quasi-local, conserved quantities of the bulk XXZ dynamics.

arXiv:1304.7944New Journal of PhysicsJul 30, 201327 citations

Matrix product solutions of boundary driven quantum chains

Tomaž Prosen

We review recent progress on constructing non-equilibrium steady state density operators of boundary driven locally interacting quantum chains, where driving is implemented via Markovian dissipation channels attached to the chain's ends. We discuss explicit solutions in three different classes of quantum chains, specifically, the paradigmatic (anisotropic) Heisenberg spin- chain, the Fermi–Hubbard chain, and the Lai–Sutherland spin-1 chain, and discuss universal concepts which characterize these solutions, such as matrix product ansatz and a more structured walking graph state ansatz. The central theme is the connection between the matrix product form of nonequilibrium states and the integrability structures of the bulk Hamiltonian, such as the Lax operators and the Yang–Baxter equation. However, there is a remarkable distinction with respect to the conventional quantum inverse scattering method, namely addressing nonequilibrium steady state density operators requires non-unitary irreducible representations of Yang–Baxter algebra which are typically of infinite dimensionality. Such constructions result in non-Hermitian, and often also non-diagonalisable families of commuting transfer operators which in turn result in novel conservation laws of the integrable bulk Hamiltonians. For example, in the case of the anisotropic Heisenberg model, quasi-local conserved operators which are odd under spin reversal (or spin flip) can be constructed, whereas the conserved operators stemming from orthodox Hermitian transfer operators (via logarithmic differentiation) are all even under spin reversal.

arXiv:1504.00783Journal of Physics A Mathematical and TheoreticalAug 19, 2015136 citations

Quantum group approach to steady states of boundary-driven open quantum systems

Enej Ilievski, Bojan Žunkovič

We present a systematic approach for constructing steady state density operators of Markovian dissipative evolution for open quantum chain models with integrable bulk interaction and boundary incoherent driving. The construction is based on fundamental solutions of the quantum Yang–Baxter equation pertaining to quantum algebra symmetries and their quantizations ( q -deformations). In particular, we facilitate a matrix-product state description, by resorting to generic spin- s infinite-dimensional solutions associated with non-compact spins, serving as ancillary degrees of freedom. After formally deriving already known steady states for the anisotropic spin-1/2 Heisenberg chain from first symmetry principles, we obtain a class of solutions belonging to interacting quantum gases with SU ( N )-symmetric Hamiltonians, using a restricted set of incoherent boundary jump processes, and point out how new non-trivial generalizations emerge from twists of quantum group structures. Finally, we discuss the possibilities for the analytical calculation of observables by employing the algebraic properties of associated auxiliary vertex operators.

arXiv:1307.5546Journal of Statistical Mechanics Theory and ExperimentJan 16, 201430 citations

Exact steady state manifold of a boundary driven spin-1 Lai–Sutherland chain

Enej Ilievski, Tomaž Prosen

We present an explicit construction of a family of steady state density matrices for an open integrable spin-1 chain with bilinear and biquadratic interactions, also known as the Lai–Sutherland model, driven far from equilibrium by means of two oppositely polarizing Markovian dissipation channels localized at the boundary. The steady state solution exhibits n+1 fold degeneracy, for a chain of length n, due to existence of (strong) Liouvillian U(1) symmetry. The latter can be exploited to introduce a chemical potential and define a grand canonical nonequilibrium steady state ensemble. The matrix product form of the solution entails an infinitely-dimensional representation of a non-trivial Lie algebra (semidirect product of sl2 and a non-nilpotent radical) and hints to a novel Yang–Baxter integrability structure.

arXiv:1402.0342Nuclear Physics BMar 23, 201434 citations

Exact Matrix Product Solution for the Boundary-Driven Lindblad XXZ Chain

Dragi Karevski, Vladislav Popkov, Gunter Markus Schutz

We demonstrate that the exact nonequilibrium steady state of the one-dimensional Heisenberg XXZ spin chain driven by boundary Lindblad operators can be constructed explicitly with a matrix product ansatz for the nonequilibrium density matrix where the matrices satisfy a quadratic algebra. This algebra turns out to be related to the quantum algebra U(q)[SU(2)]. Coherent state techniques are introduced for the exact solution of the isotropic Heisenberg chain with and without quantum boundary fields and Lindblad terms that correspond to two different completely polarized boundary states. We show that this boundary twist leads to nonvanishing stationary currents of all spin components. Our results suggest that the matrix product ansatz can be extended to more general quantum systems kept far from equilibrium by Lindblad boundary terms.

arXiv:1211.7010Physical Review LettersJan 24, 2013162 citations

Construction of the steady state density matrix and quasilocal charges for the spin-1/2 XXZ chain with boundary magnetic fields

Chihiro Matsui, Tomaž Prosen

We construct the nonequilibrium steady state (NESS) density operator of the spin-1/2 XXZ chain with non-diagonal boundary magnetic fields coupled to boundary dissipators. The Markovian boundary dis- sipation is found with which the NESS density operator is expressed in terms of the product of the Lax operators by relating the dissipation parameters to the boundary parameters of the spin chain. The NESS density operator can be expressed in terms of a non-Hermitian transfer operator (NHTO) which forms a commuting family of quasilocal charges. The optimization of the Mazur bound for the high temperature Drude weight is discussed by using the quasilocal charges and the conventional local charges constructed through the Bethe ansatz.

Mar 2, 20187 citations