© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2409.10458 · Physical review. E

Full distribution and large deviations of local observables in an exactly solvable current-carrying steady state of a strongly driven X X Z chain

Sandipan Manna, G. J. Sreejith

Published Apr 28, 2025• 0 citations

Abstract

Current-carrying steady states of interacting spin chains exhibit rich structures generated through an interplay of constraints from the Hamiltonian dynamics and those induced by the current. The XXZ spin chain when coupled to maximally polarizing Lindblad terms (with opposite signs on either end) admits an exact solution for the steady state in a matrix product state (MPS) form. We use this exact solution to study the correlations and distributions of local spin observables in the nonequilibrium steady state. We present exact expressions for spin correlators, entropy per site, and scaled cumulant generating functions (SCGF) for distributions of local observables in the XX limit (Ising anisotropy Δ=0). Further, we use the exact MPS solution in the Δ>0 regime to calculate numerically exact entropy, correlations, as well as full distributions of spin observables in large systems. In systems where Δ is a cosine of rational multiple of π, we can numerically exactly estimate the large system limit of the SCGF and the large deviation/rate functions of local-z magnetization. For these, we show that the deviations of the SCGF, calculated in finite systems, from the asymptotic large system size limit decay exponentially with system size; however, the decay rate is a discontinuous function of Δ. The x magnetization density shows a double peak structure at Δ≲1, suggesting short-range ferromagnetic ordering in the x direction similar to what was reported for the ground state of the XXZ chain.

DOI: 10.1103/physreve.111.044139

Related papers

Exact Nonequilibrium Steady State of a Strongly Driven Open XXZ Chain

Tomaž Prosen

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)].

arXiv:1106.2978•Physical Review Letters•Sep 19, 2011•290 citations

Exact Large-Deviation Statistics for a Nonequilibrium Quantum Spin Chain

Marko Žnidarič

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/Γ.

arXiv:1310.3670•Physical Review Letters•Jan 28, 2014•67 citations

Exact steady states of the impurity-doped XXZ spin chain coupled to dissipators

Chihiro Matsui, Naoto Tsuji

Abstract We give an exact matrix product steady state and matrix product forms of local observables for the bulk impurity-doped XXZ spin model coupled to dissipators at both ends, whose dynamics is described by the Lindblad quantum master equation. We find that local magnetization is induced at the impurity site when the spin current flows, which is contrary to the usual situation where the current suppresses magnetization due to heating. It is proved that this current-induced magnetization survives in the thermodynamic limit, and the spin current does not depend on the impurity strength. We discuss the role of bulk integrability by comparing the results with those of non-integrable models solved numerically by the quantum trajectory method.

arXiv:2303.17084•Journal of Statistical Mechanics Theory and Experiment•Mar 18, 2024•2 citations

Connected correlations, fluctuations and current of magnetization in the steady state of boundary driven XXZ spin chains

Berislav Buča, Tomaž Prosen

Abstract. We show how to exploit algebraic relations among the operators (or matrices) which constitute the non-equilibrium matrix product steady state of a boundary driven quantum spin chain to find partial differential equations determining all the m -point correlation functions in the continuum (or thermodynamic) limit. These partial differential equations, the order of which is determined by scaling of the non-equilibrium partition function, are readily solved if we also know the boundary conditions. In this way we can avoid resorting to representation theory of the matrix product algebra. We apply our methods to study the distributions, or moments, of the magnetization and the spin current observables in boundary driven open XXZ spin chains, as well as some connected correlation functions. Interestingly, we find that the transverse connected correlation functions are thermodynamically non-decaying and long range at the isotropic point .

arXiv:1509.04911•Journal of Statistical Mechanics Theory and Experiment•Feb 1, 2016•12 citations

Non-equilibrium steady states of the boundary-driven Heisenberg quantum spin chain

Gunter Markus Schutz

We consider an open Heisenberg quantum spin chain, coupled at the ends to boundary reservoirs polarized in different directions, which sets up a twisting gradient across the chain. We derive a matrix product ansatz for the non-equilibrium density matrix. For the isotropic chain we calculate the exact magnetization profiles and magnetization currents in the nonequilibrium steady steady state of a chain with $N$ sites. The magnetization profiles are harmonic functions with a frequency proportional to the twisting angle $\theta$. The currents of the magnetization components lying in the twisting plane and in the orthogonal direction behave qualitatively differently: In-plane steady state currents scale as $1/N^{2}$ for fixed and sufficiently large boundary coupling, and vanish as the coupling increases, while the transversal current increases with the coupling and saturates to $2\theta/N$.

•Jan 1, 2015•0 citations

Nonequilibrium steady states in the quantum XXZ spin chain

Thiago Sabetta, Grégoire Misguich

We investigate the real-time dynamics of a critical spin-1/2 chain (XXZ model) prepared in an inhomogeneous initial state with different magnetizations on the left and right halves. Using the time-evolving block decimation method, we follow the front propagation by measuring the magnetization and entanglement entropy profiles. At long times, as in the free fermion case [Antal et al., Phys. Rev. E 59, 4912 (1999)], a large central region develops where correlations become time independent and translation invariant. The shape and speed of the fronts is studied numerically and we evaluate the stationary current as a function of initial magnetic field and as a function of the anisotropy $\ensuremath{\Delta}$. We compare the results with the conductance of a Tomonaga-Luttinger liquid, and with the exact free-fermion solution at $\ensuremath{\Delta}=0$. We also investigate the two-point correlations in the stationary region and find a good agreement with the ``twisted'' form obtained by Lancaster and Mitra [Phys. Rev. E 81, 061134 (2010)] using bosonization. Some deviations are nevertheless observed for strong currents.

arXiv:1308.4851•Physical Review B•Dec 11, 2013•112 citations

Current-carrying steady states of interacting spin chains exhibit rich structures generated through an interplay of constraints from the Hamiltonian dynamics and those induced by the current. The XXZ spin chain when coupled to maximally polarizing Lindblad terms (with opposite signs on either end) admits an exact solution for the steady state in a matrix product state (MPS) form. We use this exact solution to study the correlations and distributions of local spin observables in the nonequilibrium steady state. We present exact expressions for spin correlators, entropy per site, and scaled cumulant generating functions (SCGF) for distributions of local observables in the XX limit (Ising anisotropy Δ=0). Further, we use the exact MPS solution in the Δ>0 regime to calculate numerically exact entropy, correlations, as well as full distributions of spin observables in large systems. In systems where Δ is a cosine of rational multiple of π, we can numerically exactly estimate the large system limit of the SCGF and the large deviation/rate functions of local-z magnetization. For these, we show that the deviations of the SCGF, calculated in finite systems, from the asymptotic large system size limit decay exponentially with system size; however, the decay rate is a discontinuous function of Δ. The x magnetization density shows a double peak structure at Δ≲1, suggesting short-range ferromagnetic ordering in the x direction similar to what was reported for the ground state of the 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)].

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/Γ.

Abstract We give an exact matrix product steady state and matrix product forms of local observables for the bulk impurity-doped XXZ spin model coupled to dissipators at both ends, whose dynamics is described by the Lindblad quantum master equation. We find that local magnetization is induced at the impurity site when the spin current flows, which is contrary to the usual situation where the current suppresses magnetization due to heating. It is proved that this current-induced magnetization survives in the thermodynamic limit, and the spin current does not depend on the impurity strength. We discuss the role of bulk integrability by comparing the results with those of non-integrable models solved numerically by the quantum trajectory method.

Abstract. We show how to exploit algebraic relations among the operators (or matrices) which constitute the non-equilibrium matrix product steady state of a boundary driven quantum spin chain to find partial differential equations determining all the m -point correlation functions in the continuum (or thermodynamic) limit. These partial differential equations, the order of which is determined by scaling of the non-equilibrium partition function, are readily solved if we also know the boundary conditions. In this way we can avoid resorting to representation theory of the matrix product algebra. We apply our methods to study the distributions, or moments, of the magnetization and the spin current observables in boundary driven open XXZ spin chains, as well as some connected correlation functions. Interestingly, we find that the transverse connected correlation functions are thermodynamically non-decaying and long range at the isotropic point .

We consider an open Heisenberg quantum spin chain, coupled at the ends to boundary reservoirs polarized in different directions, which sets up a twisting gradient across the chain. We derive a matrix product ansatz for the non-equilibrium density matrix. For the isotropic chain we calculate the exact magnetization profiles and magnetization currents in the nonequilibrium steady steady state of a chain with

N

sites. The magnetization profiles are harmonic functions with a frequency proportional to the twisting angle

\theta

. The currents of the magnetization components lying in the twisting plane and in the orthogonal direction behave qualitatively differently: In-plane steady state currents scale as

1/N^{2}

for fixed and sufficiently large boundary coupling, and vanish as the coupling increases, while the transversal current increases with the coupling and saturates to

2\theta/N

.

We investigate the real-time dynamics of a critical spin-1/2 chain (XXZ model) prepared in an inhomogeneous initial state with different magnetizations on the left and right halves. Using the time-evolving block decimation method, we follow the front propagation by measuring the magnetization and entanglement entropy profiles. At long times, as in the free fermion case [Antal et al., Phys. Rev. E 59, 4912 (1999)], a large central region develops where correlations become time independent and translation invariant. The shape and speed of the fronts is studied numerically and we evaluate the stationary current as a function of initial magnetic field and as a function of the anisotropy $\ensuremath{\Delta}$. We compare the results with the conductance of a Tomonaga-Luttinger liquid, and with the exact free-fermion solution at $\ensuremath{\Delta}=0$. We also investigate the two-point correlations in the stationary region and find a good agreement with the ``twisted'' form obtained by Lancaster and Mitra [Phys. Rev. E 81, 061134 (2010)] using bosonization. Some deviations are nevertheless observed for strong currents.