# Duality in quantum transport models

Rouven Frassek, Cristian Giardinà, Jorge Kurchan

► To cite this version:

Rouven Frassek, Cristian Giardinà, Jorge Kurchan. Duality in quantum transport models. SciPost Physics, 2021, 10 (6), pp.135. 10.21468/SciPostPhys.10.6.135. hal-02981263

HAL Id: hal-02981263

https://hal.science/hal-02981263v1

Submitted on 23 Aug 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Duality in quantum transport models

Rouven Frassek¹*, Cristian Giardinà²† and Jorge Kurchan¹‡

1 Laboratoire de Physique de l'École Normale Supérieure, ENS, Université PSL, CNRS, Sorbonne Université, Université de Paris, 75005 Paris, France

2 University of Modena and Reggio Emilia, FIM, Via G. Campi 213/b, 41125 Modena, Italy

* rouven.frassek@phys.ens.fr, † cristian.giardina@unimore.it,

‡ kurchan.jorge@gmail.com

## Abstract

We develop the ‘duality approach’, that has been extensively studied for classical models of transport, for quantum systems in contact with a thermal ‘Lindbladian’ bath. The method provides (a) a mapping of the original model to a simpler one, containing only a few particles and (b) shows that any dynamic process of this kind with generic baths may be mapped onto one with equilibrium baths. We exemplify this through the study of a particular model: the quantum symmetric exclusion process introduced in [1]. As in the classical case, the whole construction becomes intelligible by considering the dynamical symmetries of the problem.

<figure><img src="image_2.png" /></figure>

Copyright R. Frassek *et al.*

This work is licensed under the Creative Commons

Attribution 4.0 International License.

Published by the SciPost Foundation.

Received 22-10-2020

Accepted 28-05-2021

Published 08-06-2021

doi:10.21468/SciPostPhys.10.6.135

<figure><img src="image_3.png" /></figure>

Check for
updates

## Contents

<table><tr><td>1</td><td>Introduction</td><td>2</td></tr><tr><td>2</td><td>Outline of duality technique</td><td>3</td></tr><tr><td>3</td><td>A quantum transport model</td><td>4</td></tr><tr><td>4</td><td>Real replicas and expectation values</td><td>6</td></tr><tr><td>5</td><td>Explicit construction of a bath</td><td>6</td></tr><tr><td>6</td><td>Transformations</td><td>7</td></tr><tr><td>7</td><td>Billiard pocket</td><td>8</td></tr><tr><td>8</td><td>Correlation functions & real replicas</td><td>9</td></tr><tr><td>9</td><td>Mapping the driven system into one satisfying detailed balance</td><td>11</td></tr><tr><td>10</td><td>Conclusion</td><td>12</td></tr><tr><td></td><td>A Algebra and representations</td><td>12</td></tr><tr><td></td><td>References</td><td>13</td></tr></table>
# 1 Introduction

A central role in non-equilibrium statistical mechanics of classical systems is played by stochastic processes. For instance, in the study of mass transport associated to a non-equilibrium steady state with non-zero current, a pivotal role has been played by the *simple symmetric exclusion process* (SSEP), that has been the subject of intensive investigations since its introduction. In the study of this process emerge properties that are believed to be universal signatures of a non-equilibrium stationary states, such as long-range correlations (in turn the source of non-local large deviation functionals for the density). Similarly, non-gaussian fluctuations are observed in the asymmetric version of the process (related to the Kardar-Parisi-Zhang universality class). In both settings several tools from the theory of Markov processes have been used. We shall focus here on one such a tool, which is known in the probabilistic literature as *duality*.

Dual processes were introduced in the realm of interacting particle systems at the early days of the field by Spitzer and Liggett [2, 3]. For instance, the symmetric exclusion process on the lattice turns out to be *self-dual*, and this property has been heavily used and fundamental to develop the ergodic theory of exclusion process in infinite volume. In the non-equilibrium set-up, dual processes are also useful to study the stationary measure, which is necessarily non-reversible to sustain a current. Again, the simplest example is the open exclusion process on a chain, which is coupled at its ends to reservoirs which inject and remove particles at different rates. Here the dual process simplifies the analysis by *transforming the reservoirs into absorbing boundaries* [4]. In doing so, the study of correlation functions in the open system is reduced to following the dynamics of a few dual particles that are eventually absorbed in the boundaries.

Dual counterparts of the open exclusion process played a crucial role in the construction of the so-called *hydrodynamic limit* [5], i.e. a macroscopic theory described by partial differential equations. This turned out to be true for a large class of diffusive systems (Kipnis-Marchioro-Presutti process, symmetric inclusion process, Brownian energy process [6–11]) related to the study of Fourier’s law of heat conduction. Other simplifications due to duality occurred in the study of asymmetric exclusion process on the infinite line [12–15], directly related to the height profile of interface growth models. There, duality helps the study of *fluctuations* around the hydrodynamic limit: the evolution of one dual particle is related to the microscopic version of the Cole/Hopf transform that maps the non-linear and ill-posed KPZ equation to the linear stochastic heat equation [15, 16].

Duality has also another surprising consequence. Years ago, Tailleur et al. [17] were puzzled by a crucial step associated with the solution of Bertini et al [18] of the *large deviations* (around the hydrodynamic limit) of a family of transport models: the fact that the explicit construction of a trajectory with time-reversed extremes was possible even out of equilibrium. The solution of the puzzle was that these models revert, via a non-local transformation, into problems with detailed balance. In a recent paper [19], we have shown that this is a general consequence of duality, to be expected if (and probably only if) some form of duality is present.

It is then natural to ask if duality may be introduced, and if so, which (if any) of the simplifications obtained by such a construction survive in a *quantum stochastic system* [1, 20–24]. This is precisely the question addressed in this paper. At first sight, the two problems may seem completely different, as the replacement of the Markov evolution with a quantum semigroup substantially changes the averaging procedure: quantum averages are quadratic in the wavefunction, whereas the probability density appears linearly in the averages of classical stochastic systems. However, as we shall further discuss below, there is a key argument that brings the two problems in close contact: in both settings, classical and quantum, the evolution
<figure><img src="image_4.png" /><figcaption>Figure 1: Schematic picture of a transport model on a graph. A non-equilibrium stationary state is attained as a consequence of interactions with sources at different parameters. The example shows three 'baths' with couplings $\Gamma_i$ and average densities of particles $\lambda_i$, with $i = 1, 2, 3$.</figcaption></figure>

can be described in algebraic terms using a Lie algebra. For classical systems this was remarked years ago in a pioneering paper by Schütz and Sandow [25] and then further extended in [26]. Here we show, by focusing on a specific example, that quantum stochastic systems admits the same algebraic description, provided one moves to *superoperators*. This allows to repeat the pattern that lead to the formulation of duality in Markov systems.

The example we shall discuss is a system of free fermions with a noisy dynamics which has been extensively studied in [1, 23, 24]. The system plays a paradigmatic role for quantum stochastic evolutions, similar to the exclusion process of classical systems, and indeed it has been named the *quantum exclusion process* (we will discuss bosons in a separate paper). For simplicity we restrict here our discussion of duality to the boundary-driven setting. We believe our construction of a dual process extends to other quantum models as well, particularly those of exactly solvable quantum Liouvillians (or Lindbladians), beyond the quantum symmetric exclusion (e.g. [27–30]).

In a quantum stochastic systems two sources of randomness coexist: *quantum fluctuations* (due to the quantum evolution) and *dynamical fluctuations* (due to a noisy dynamics). These play a role similar to thermodynamic fluctuations and quenched disorder, in disordered systems.

## 2 Outline of duality technique

The duality technique may be described through the following steps:

* We start with a general graph, not necessarily (although often) a chain. Between the vertices of the graph particles or energy are transported, the amount per unit time is a stochastic process (Figure 1).

* We connect 'leads' of the graph to sources of heat or particles: the 'baths'. These may be constructed from first principles by considering each bath as a very large (in fact, infinite) equilibrium system. Transport occurs when the equilibria of the baths are incompatible (Figure 2).

* In all systems where duality has been introduced, it happens that the Hamiltonian generating the dynamics may be written in terms of the generators of a (non-abelian, Lie)
group. The bulk is invariant under the group operations, while the baths are not. Acting with the group hence only modifies the baths' properties.

* Like in any stochastic system, we may compute expectation values of observables. We may then switch to the 'adjoint', in which one evolves the observables rather than the state of the system. (cf. going from Schrödinger to Heisenberg pictures).

* Crucially: there is a judicious combination of group operation and passing to the adjoint that maps the system into a system with purely absorbing baths, and only a few particles. This is the 'dual' setting.

* Stationary results for the original system are retrieved from the long-time results, when the baths have emptied completely the system, by counting how many particles (or how much energy) was absorbed by each lead.

We shall outline in what follows how these steps are implemented in the quantum case.

### 3 A quantum transport model

On a graph $G = (V, E)$ with vertex set $V$ and edge set $E$, we consider a set of fermionic creation/annihilation operators and the random Hamiltonian

$$H_{\eta}(t) = \sum_{(k,l) \in E} \sqrt{c_{kl}} [a_{k}^{\dagger} a_{l} \eta_{kl}(t) + a_{l}^{\dagger} a_{k} \bar{\eta}_{kl}(t)], \quad (1)$$

where $\{a_i, a_j\} = 0$, $\{a_i^{\dagger}, a_j^{\dagger}\} = 0$, $\{a_i, a_j^{\dagger}\} = \delta_{i,j}$. This describes a system of free fermions jumping on the graph $G$, the $c_{kl}$ are coupling constants. The jump terms are noisy: the external quenched noise is given by pairs of independent and identical distributed complex conjugated Gaussian white noise (one pair for each edge of the graph) with covariances

$$\mathbb{E}[\eta_{kl}(t) \bar{\eta}_{k'l'}(t')] = \delta(t-t') \delta_{kl,k'l'}. \quad (2)$$

Using the Trotter's product formula, it then follows for the density matrix

$$\rho_{\eta}(t) = T \left[ \exp \left\{ -i \int_{t_0}^{t} [H_{\eta}(t'), \cdot] dt' \right\} \right] \rho(0), \quad (3)$$

where $T$ denotes "time-order". Expanding we have

$$\rho_{\eta}(t) = T \left[ \left( 1 - i \int_{t_0}^{t} dt_1 [H_{\eta}(t_1), \cdot] - \frac{1}{2} \int_{t_0}^{t} dt_1 \int_{t_0}^{t} dt_2 [H_{\eta}(t_2), [H_{\eta}(t_1), \cdot]] + \dots \right) \right] \rho(0). \quad (4)$$

We are interested in the noise-dependent expectation

$$\langle A(t) \rangle_{\eta} = \mathrm{Tr}[A \rho_{\eta}(t)], \quad (5)$$

where $A$ is a generic operator (which may be expressed in terms of $a_i$ and $a_i^{\dagger}$) as well as the *averaged-quenched* expectation

$$\mathbb{E}[\langle A(t) \rangle_{\eta}] = \mathrm{Tr}[A \rho(t)], \quad (6)$$
where we have defined the quenched-averaged density matrix

$$ \rho(t) = \mathbb{E}[\rho_{\eta}(t)]. \qquad (7) $$

In the formulas above 'Tr' denotes the trace operation yielding the quantum expectations. Averaging (4) with respect to the external noise and using the covariances (2) we get the evolution equation [1]

$$ \begin{aligned} \frac{d}{dt}\rho &= -\frac{1}{2} \sum_{(k,\ell)\in E} c_{kl}([a_k^\dagger a_\ell, [a_\ell^\dagger a_k, \rho]] + [a_\ell^\dagger a_k, [a_k^\dagger a_\ell, \rho]]) \\ &\equiv -\mathcal{H}(\rho). \end{aligned} \qquad (8) $$

Developing the commutators $\mathcal{H}$ may be written as

$$ \mathcal{H} = - \sum_{(k, \ell) \in E} c_{kl} \left( \mathcal{J}_k^+ \mathcal{J}_\ell^- + \mathcal{J}_k^- \mathcal{J}_\ell^+ + 2C_k^+ C_\ell^+ + 2C_k^- C_\ell^- - \frac{1}{2} \right), \qquad (9) $$

where we have defined the following on-site superoperators acting as follows on an operator A:

$$ \mathcal{J}_i^+(A) = a_i^\dagger A a_i, \qquad (10) $$

$$ \mathcal{J}_i^-(A) = a_i A a_i^\dagger, $$

$$ C_i^+(A) = \frac{1}{2}(a_i^\dagger a_i A - A a_i^\dagger a_i), $$

$$ C_i^-(A) = \frac{1}{2}(a_i^\dagger a_i A + A a_i^\dagger a_i - A). $$

Just by applying two of these in succession (in both orders), one can directly check they satisfy (on each vertex i) an u(2) algebra, decomposable as su(2)

$$ [C_i^-, \mathcal{J}_i^+] = \pm \mathcal{J}_i^+, \quad [\mathcal{J}_i^+, \mathcal{J}_i^-] = 2C_i^-, \qquad (11) $$

and $C_i^+$, which commutes with everything $[C_i^+, \mathcal{J}_i^+] = 0 = [C_i^+, C_i^-]$ and are hence constants of motion for every i. It is easy to see that they count the number of creation minus the number of destruction operators in every site. Throughout this paper, this number is zero: there is always an equal number of creation and destruction operators in all sites. Other representations are of course possible.

We get, for the averaged-quenched expectation value of an operator A at time t:

$$ \mathbb{E}[\langle A(t) \rangle_{\eta}] = \mathrm{Tr}[Ae^{-t\mathcal{H}}\rho(0)] = \mathrm{Tr}[\rho^{\dagger}(0)e^{-t\mathcal{H}^{\dagger}}A^{\dagger}]^{*}, \qquad (12) $$

where we have introduced the adjoint $\mathcal{D}^\dagger$ of a superoperator $\mathcal{D}$ as:

$$ \mathrm{Tr}[A^{\dagger}\mathcal{D}(B)]^{*} = \mathrm{Tr}[B^{\dagger}\mathcal{D}^{\dagger}(A)], \quad \forall A, B. \qquad (13) $$

Note that the same symbol is used both for the adjoint $A^\dagger$ of an operator A on the Hilbert space, and for the adjoint of superoperators.

Using the cyclic property of the trace, it is easy to see that:

$$ [\mathcal{J}_i^{\pm}]^{\dagger} = \mathcal{J}_i^{\mp} \quad \text{and} \quad [C_i^{\pm}]^{\dagger} = C_i^{\pm}. \qquad (14) $$
<figure><img src="image_5.png" /><figcaption>Figure 2: Schematic picture of a reservoir connected to site 1 of the graph. The bath is made of infinitely many vertices in a papyrus fashion (in the figure we represent some of them as the points connected from the left to site 1) with the same rule of transport of the bulk system. Each bath is an equilibrium system.</figcaption></figure>

## 4 Real replicas and expectation values

What we have done up to now is not enough [20, 21]. Indeed, suppose we wish to calculate the following expectation: $E[\langle A(t) \rangle_{\eta} \langle B(t) \rangle_{\eta}] = E\{\text{Tr}[A\rho_{\eta}(t)]\text{Tr}[B\rho_{\eta}(t)]\}$. A way to do this is to replicate the system twice (all superoperators, operators and states) and use the fact that the trace of a tensor product is the product of the traces:

$$E[\langle A(t) \rangle_{\eta} \langle B(t) \rangle_{\eta}] = E\{\text{Tr}[A^{\alpha}B^{\beta}\rho_{\eta}^{(\alpha)} \otimes \rho_{\eta}^{(\beta)}(t)]\}, \quad (15)$$

where the 'replica index' $\alpha, \beta = 1, 2$. This will lead us to product density functions $\rho_2 = E[\rho_{\eta}^{(1)} \otimes \rho_{\eta}^{(2)}]$ of two real replicas, and an averaged operator $\mathcal{H}_2(a_{i,\alpha}^{\dagger}, a_{i,\alpha}, a_{i,\beta}^{\dagger}, a_{i,\beta})$ acting on it. The noise is the same for both replicas. Note the close analogy with disordered systems: here the fermions are the replicated variables (like the spins of a spin-glass), and the noise is the disorder, playing the role of the disordered interactions $J_{ij}$. For simplicity, we shall run the steps in detail for a single replica, and then present in a more compact way the same steps for the n-times replicated case.

## 5 Explicit construction of a bath

We shall construct the bath as an ensemble of $s$ independent sites represented by fermionic operators $b_m^{\dagger}, b_m$ connected to site number 1 in a papyrus fashion like in Figure 2. The bath is larger than any other number in the problem ($s \to \infty$), and its particle density is fixed on average by its initial density matrix:

$$\rho_{\text{bath}} = \frac{e^{-\mu \sum_{m=1}^{s} b_m^{\dagger} b_m}}{(1 + e^{-\mu})^s}, \quad (16)$$

which we should tensor-product with the initial one of the chain.

Thus we consider a sum of $s$ interaction terms, which we shall assume are of intensity $c_{1,m} = \Gamma_1/s$, with $\Gamma_1 > 0$ the coupling constant to the reservoir. Denoting $\mathcal{K}^{\pm} = (\sum_{m \in \text{bath}} \mathcal{J}_{m}^{\pm})/s$ and $\mathcal{D}^{\pm} = (\sum_{m \in \text{bath}} \mathcal{C}_{m}^{\pm})/s$, the combined jump superoperators to all the leaves in the bath (which obviously obey the same algebra), we may write:

$$\mathcal{H}_{1}^{(s)} = -\Gamma_{1} \left( \mathcal{K}^{+} \mathcal{J}_{1}^{-} + \mathcal{K}^{-} \mathcal{J}_{1}^{+} + 2\mathcal{D}^{+} \mathcal{C}_{1}^{+} + 2\mathcal{D}^{-} \mathcal{C}_{1}^{-} - \frac{1}{2} \right). \quad (17)$$
Now, let us consider the action of $\mathcal{K}^{\pm}$, $\mathcal{D}^{\pm}$ on $\rho_{bath}$. For example,

$$ \mathcal{K}^{+}(\rho_{bath}) = \frac{1}{s} \sum_{m} b_{m}^{\dagger} \rho_{bath} b_{m} = \frac{1}{s} e^{\mu} \sum_{m} b_{m}^{\dagger} b_{m} \rho_{bath}. $$

Because $s$ is of large, by the law of large numbers we may substitute the number of fermions in the bath by its average: $\frac{1}{s} \sum_{m=1}^{s} b_{m}^{\dagger} b_{m} \rightarrow \frac{e^{-\mu}}{1+e^{-\mu}} \equiv \lambda$ as $s \rightarrow \infty$. This implies that we may substitute, for large $s$

$$ \mathcal{K}^{+} \rightarrow 1 - \lambda, \quad \mathcal{K}^{-} \rightarrow \lambda, \qquad (18) $$

$$ \mathcal{D}^{-} \rightarrow \frac{2\lambda - 1}{2}, \quad \mathcal{D}^{+} \rightarrow 0. $$

The net effect is that (now allowing for several leads $i$):

$$ \mathcal{H}_{bath} = - \sum_{i \in V} \Gamma_i \left[ \lambda_i \left( \mathcal{J}_i^+ + \mathcal{C}_i^- - \frac{1}{2} \right) + (1 - \lambda_i) \left( \mathcal{J}_i^- - \mathcal{C}_i^- - \frac{1}{2} \right) \right], \qquad (19) $$

which is clearly of the Lindbladian form.

# 6 Transformations

Bearing the group structure in mind, the transformations leading to a dual process are now readable from the analogous ones in the classical system [19]. To simplify we start considering just one lead to a bath $i = 1$

$$ (\mathcal{H}_1)^{\dagger} = - \left[ \lambda(\mathcal{J}_1^- + \mathcal{C}_1^- - \frac{1}{2}) + (1-\lambda)(\mathcal{J}_1^+ - \mathcal{C}_1^- - \frac{1}{2}) \right], \qquad (20) $$

where we have used (14). Now conjugate with $e^{\mathcal{J}_1^+}$ and find

$$ e^{-\mathcal{J}_1^+} (\mathcal{H}_1)^{\dagger} e^{\mathcal{J}_1^+} = -\lambda \mathcal{J}_1^- + \mathcal{C}_1^- + \frac{1}{2}. $$

We now consider an extended system and write $\mathcal{H} = \mathcal{H}_{bulk} + \mathcal{H}_{baths}$. Doing the same with every lead $i$, and using that the bulk superoperator $\mathcal{H}_{bulk}$ commutes with $\mathcal{J}_{tot}^+ = \sum_{i \in V} \mathcal{J}_i^+$ we get

$$ \begin{aligned} \mathcal{H}^{\dagger} &= e^{\mathcal{J}_{tot}^{+}} \left( \mathcal{H}_{bulk} - \sum_{i \in V} \Gamma_i \left\{ \lambda_i \mathcal{J}_i^{-} - \mathcal{C}_i^{-} - \frac{1}{2} \right\} \right) e^{-\mathcal{J}_{tot}^{+}} \\ &= e^{\mathcal{J}_{tot}^{+}} \mathcal{H}' e^{-\mathcal{J}_{tot}^{+}}. \end{aligned} \qquad (21) $$

Let us now reinstate baths, undoing the step we did before. We thus attach to each vertex $i \in V$ an extra site with superoperators $\mathcal{K}_i^{\pm}$ and $\mathcal{D}_i^{\pm}$ describing now empty baths (i.e. all having $\lambda = 0$, cf. Eq. (18))

$$ \lambda_i \mathcal{J}_i^{-} - \mathcal{C}_i^{-} - \frac{1}{2} \rightarrow \lambda_i \mathcal{K}_i^{+} \mathcal{J}_i^{-} + 2 \mathcal{D}_i^{-} \mathcal{C}_i^{-} - \frac{1}{2}. \qquad (22) $$

Comparing to (17) we notice that the terms $\mathcal{K}_i^{-} \mathcal{J}_i^{+}$ and $\mathcal{D}_i^{+} \mathcal{C}_i^{-}$ are absent because the baths are now empty. We stress once more that, although labeled with the indexes $i$ of the vertices of the graphs, the superoperators describing the empty bath live in an additional extra space (i.e. they act on vertices attached to the bath). Finally another conjugation allows to eliminate
the disturbing factors $\lambda_i$. Indeed, using the fact that $e^{\ln \lambda_i (\mathcal{D}_i^- + 1/2)} \mathcal{K}_i^+ e^{-\ln \lambda_i (\mathcal{D}_i^- + 1/2)} = \lambda_i \mathcal{K}_i^+$ we may write

$$ \mathcal{H}' = e^{\sum_i \ln \lambda_i (\mathcal{D}_i^- + 1/2)} \mathcal{H}^{\text{dual}} e^{-\sum_i \ln \lambda_i (\mathcal{D}_i^- + 1/2)}, \qquad (23) $$

where

$$ \mathcal{H}^{\text{dual}} \equiv \mathcal{H}_{\text{bulk}} - \sum_{i \in V} \Gamma_i \left\{ \mathcal{K}_i^+ \mathcal{J}_i^- + 2 \mathcal{D}_i^- \mathcal{C}_i^- - \frac{1}{2} \right\}. \qquad (24) $$

All in all, the conjugation in (21) can thus be used to write the expectation of an observable $O$ as:

$$ \begin{aligned} \mathrm{Tr}[Oe^{-t\mathcal{H}}\rho(0)] &= \mathrm{Tr}[\rho(0)^{\dagger}e^{-t\mathcal{H}^{\dagger}}O^{\dagger}]^{*} \\ &= \mathrm{Tr}[\rho(0)^{\dagger}e^{\mathcal{J}_{\mathrm{tot}}^{+}}e^{-t\mathcal{H}'}e^{-\mathcal{J}_{\mathrm{tot}}^{+}}O^{\dagger}]^{*}. \end{aligned} \qquad (25) $$

As we shall later see, it is convenient to define

$$ \hat{O} \equiv e^{-\mathcal{J}_{\mathrm{tot}}^{+}} O^{\dagger} = \Pi_i e^{-\mathcal{J}_i^{+}} O^{\dagger}. \qquad (26) $$

Note that because of the product form, $\hat{O}$ depends exclusively on the same sites as $O$. We have to make a choice of how we extend $O$ in the product space. We shall choose (without loss of generality) that it is $\tilde{O} = \hat{O} \otimes |0_b\rangle\langle 0_b|$: the bath sites are completely empty at time zero. Instead, the operator $\rho(0)$ is now understood as $\rho(0) \otimes 1_b$, it acts on the bath sites as the identity. Suppose for example that $O = a_k^\dagger a_k$ for a chain of length $N$, then

$$ \begin{aligned} \Pi_i e^{-\mathcal{J}_i^+} (a_k^\dagger a_k)^\dagger &= (1 - a_1^\dagger a_1) \dots a_k^\dagger a_k \dots (1 - a_N^\dagger a_N) \\ &= |0\rangle \langle 0|_1 \dots \otimes |1\rangle \langle 1|_k \dots \otimes |0\rangle \langle 0|_N, \end{aligned} \qquad (27) $$

i.e. a chain with one fermion in site $k$, and otherwise empty.

Going back to (25) and using once again the adjoint in (14), we get

$$ \mathrm{Tr}[Oe^{-t\mathcal{H}}\rho(0)] = \mathrm{Tr}\left[(e^{\mathcal{J}_{\mathrm{tot}}^{-}}(\rho(0)))^{\dagger} e^{-t\mathcal{H}'}(\hat{O})\right]^{*}, \qquad (28) $$

and inserting (23) into (28) we obtain the final result

$$ \mathbb{E}[\langle O(t) \rangle_{\eta}] = \mathrm{Tr}_{\mathrm{bulk}} \mathrm{Tr}_{\mathrm{baths}} \left[ \left\{ \left( e^{\mathcal{J}_{\mathrm{tot}}^{-}}(\rho(0)) \right)^{\dagger} e^{\sum_{i} \ln \lambda_{i} (\mathcal{D}_{i}^{-} + 1/2)} \right\} e^{-t \mathcal{H}^{\mathrm{dual}}} (\tilde{O}) \right]^{*}. \qquad (29) $$

We have used the fact that $\tilde{O}$ has no particles in the bath sites, so that $(\mathcal{D}_i^- + 1/2)(\tilde{O}) = 0$ (cf. Eq. (10)). This is the duality result: to compute the expectation of the observable $O$ evolving with the original process and starting from the density matrix $\rho(0)$, we consider instead an initial density matrix $\tilde{O}$ which evolves through the dual process, and we evaluate the ‘observable’ $(e^{\mathcal{J}_{\mathrm{tot}}^{-}}(\rho(0)))^{\dagger} e^{\sum_i (\ln \lambda_i) (\mathcal{D}_i^- + 1/2)}$ at the end. Note that if $\tilde{O}$ has not the requirements of unit trace and positivity of a density matrix, it can always be brought into one that has those properties, by addition of a term proportional to the invariant measure, and multiplication by a constant (to normalize).

# 7 Billiard pocket

Duality relations show their power when considering infinite time evolution. When $t \to \infty$ the evolution voids the chain. We may thus expect that $\lim_{t \to \infty} e^{-t\mathcal{H}^{\mathrm{dual}}}(\tilde{O}) = [\text{empty bulk}] \otimes O_\infty$,
where $O_\infty$ lives in the space of bath sites exclusively. Equation (29) becomes, then:

$$ \lim_{t \to \infty} \mathbb{E}[\langle O(t) \rangle_\eta] = \langle 0 | e^{\mathcal{J}_{\text{tot}}}(\rho(0)) | 0 \rangle_{\text{bulk}} \text{Tr}_{\text{baths}} \left[ e^{\sum_i \ln \lambda_i (\mathcal{D}_i^{-1/2})} O_\infty \right]^* \\ = \text{Tr}_{\text{baths}} \left[ e^{\sum_i \ln \lambda_i (\mathcal{D}_i^{-1/2})} O_\infty \right]^* \quad (30) $$

where we have used $\langle 0 | (e^{\mathcal{J}_{\text{tot}}}(\rho(0)))^\dagger | 0 \rangle_{\text{bulk}} = \text{Tr}_{\text{bulk}}[\rho(0)] = 1$. All the interesting information is stored in the matrix $O_\infty$. Because $O_\infty$ is a combination of bath sites that are either void or have one fermion, and denoting $n_i$ the total number of fermions in the bath $i$, we get the simple expression:

$$ \lim_{t \to \infty} \text{Tr}[Oe^{-t\mathcal{H}}\rho(0)] = \sum_{\{n_i\}} c_{\{n_i\}}^{[0]} \Pi_i \lambda_i^{n_i}, \quad (31) $$

where the coefficients $c_{\{n_i\}}^{[0]}$ depend on the observable $O$. We shall see a concrete example of this presently.

**Example.** Consider a linear chain of length $N$ connected at the extremes (denoted by 1 and $N$) to two reservoirs at densities $\lambda_L$ on the left, respectively $\lambda_R$ on the right. Suppose we are interested in the average number of “particles” at site $k$ in the stationary state. For this problem we now put $O_k = a_k^\dagger a_k$ (see Eq. (27)) and $i \in \{L, R\}$. We start the dual evolution from the density operator associated to the pure state with only one fermion at site $k$. In the long-time limit we will have either a fermion in the left or a fermion in the right bath, with probabilities $c_L^{[O_k]}$ and $c_R^{[O_k]} = 1 - c_L^{[O_k]}$, respectively. Since for a symmetric random walk

$$ c_L^{[O_k]} = 1 - \frac{k}{N+1} \quad (32) $$

we conclude from (31) the stationary linear profile

$$ \langle a_k^\dagger a_k \rangle = \lambda_L + \frac{\lambda_R - \lambda_L}{N+1} k, \quad k \in \{1, \dots, N\}. \quad (33) $$

A similar result can be obtained on a generic graph $G$, replacing the $c_L^{[O_k]}$ with the harmonic function of the symmetric random walk on the graph $G$.

# 8 Correlation functions & real replicas

We introduce replicas (labeled by 1, 2, ...) by considering copies of the system subject to same realization of the external noise $\eta$ and characterized by fermionic operators $a_{i,1}, a_{i,1}^\dagger, a_{i,2}, a_{i,2}^\dagger, \dots$ (in the bulk) and $b_{im,1}, b_{im,1}^\dagger, b_{im,2}, b_{im,2}^\dagger, \dots$ (in the bath).

A generic number $n$ of copies is described by the noise-dependent density operator given by the tensor product and we are interested in its average:

$$ \rho_n = \mathbb{E}[\rho_{\eta}^{(1)} \otimes \rho_{\eta}^{(2)} \otimes \dots \otimes \rho_{\eta}^{(n)}], \quad (34) $$

where each $\rho_{\eta}^{(\alpha)}$ acts on $a_{i,\alpha}, a_{i,\alpha}^\dagger, b_{im,\alpha}, b_{im,\alpha}^\dagger$'s, which are coupled by the same realization of the external noise. This system evolves with the replicated noisy Hamiltonian $H_\eta(t) = \sum_{\alpha=1}^n (\sum_{(k,l) \in E} [a_{k,\alpha}^\dagger a_{l,\alpha} \eta_{kl}(t) + a_{l,\alpha}^\dagger a_{k,\alpha} \bar{\eta}_{kl}(t)])$. Developing up to second order, and averaging over the noise we get the evolution equation $\frac{d}{dt}\rho_n = -\mathcal{H}_n \rho_n$ with

$$ \mathcal{H}_n = - \sum_{\alpha, \beta=1}^{n} \sum_{(k,l) \in E} c_{kl} (\mathcal{J}_{k,\alpha\beta}^{+} \mathcal{J}_{l,\beta\alpha}^{-} + \mathcal{J}_{k,\alpha\beta}^{-} \mathcal{J}_{l,\beta\alpha}^{+}) + 2C_{k,\alpha\beta}^{+} C_{l,\beta\alpha}^{+} + 2C_{k,\alpha\beta}^{-} C_{k,\beta\alpha}^{-} - \frac{\delta_{\alpha\beta}}{2} \quad (35) $$
obtained by expanding commutators just as above. The operators are now:

$$ \mathcal{J}_{i,\alpha\beta}^{+}(A) = a_{i,\alpha}^{\dagger} A a_{i,\beta}, \qquad (36) $$

$$ \mathcal{J}_{i,\alpha\beta}^{-}(A) = a_{i,\beta} A a_{i,\alpha}^{\dagger}, $$

$$ C_{i,\alpha\beta}^{+}(A) = \frac{1}{2}(a_{i,\alpha}^{\dagger} a_{i,\beta} A - A a_{i,\alpha}^{\dagger} a_{i,\beta}), $$

$$ C_{i,\alpha\beta}^{-}(A) = \frac{1}{2}(a_{i,\alpha}^{\dagger} a_{i,\beta} A + A a_{i,\alpha}^{\dagger} a_{i,\beta} - A \delta_{\alpha\beta}), $$

which satisfy a u(2n) algebra. The operators $\sum_{\alpha} C_{i,\alpha\alpha}^{+}$ commute with everything and thus do not evolve, they count the total number of creators minus annihilators per site. We shall only be interested in the case they are all zero. $\mathcal{H}_n$ is a nearest-neighbor quantum su(2n) chain.

We may now distinguish bath terms, and we obtain in an entirely analogous manner:

$$ \begin{aligned} \mathcal{H}_{n,bath} = & -\sum_{\alpha,\beta=1}^{n} \sum_{i \in V} \Gamma_i [\mathcal{K}_{i,\alpha\beta}^{+} \mathcal{J}_{i,\alpha\beta}^{-} + \mathcal{K}_{i,\alpha\beta}^{-} \mathcal{J}_{i,\alpha\beta}^{+} + 2\mathcal{D}_{i,\alpha\beta}^{+} C_{i,\alpha\beta}^{+} + 2\mathcal{D}_{i,\alpha\beta}^{-} C_{i,\alpha\beta}^{-} - \frac{1}{2}\delta_{\alpha,\beta}] \\ & -\sum_{\alpha,\beta=1}^{n} \sum_{i \in V} \Gamma_i [\lambda_{i,\alpha\beta}(\mathcal{J}_{i,\alpha\beta}^{+} + C_{i,\alpha\beta}^{-} - \frac{\delta_{\alpha,\beta}}{2}) + (\delta_{\alpha\beta} - \lambda_{i,\alpha\beta})(\mathcal{J}_{i,\alpha\beta}^{-} - C_{i,\alpha\beta}^{-} - \frac{\delta_{\alpha,\beta}}{2})], \end{aligned} \qquad (37) $$

where we have substituted as before bath operators by their expectation values:

$$ \begin{aligned} \mathcal{K}_{i,\alpha\beta}^{+} &\rightarrow \delta_{\alpha\beta} - \lambda_{i,\alpha\beta}, & \mathcal{K}_{i,\alpha\beta}^{-} &\rightarrow \lambda_{i,\alpha\beta}, \\ \mathcal{D}_{i,\alpha\beta}^{-} &\rightarrow \frac{2\lambda_{i,\alpha\beta} - \delta_{\alpha\beta}}{2}, & \mathcal{D}_{i,\alpha\beta}^{+} &\rightarrow 0. \end{aligned} \qquad (38) $$

The $\lambda_{i,\alpha\beta}$ define the bath. The most general replica-symmetric form is $\lambda_{i,\alpha\beta} = \lambda_i \delta_{\alpha\beta} + \lambda'_i$. Matrices with different $(\lambda_i, \lambda'_i)$ commute and may be diagonalized simultaneously for all baths by a rotation in the fermion space. We shall in fact not need to do this here for the following reason: one can easily show that for $\alpha \neq \beta$, $\text{Tr}[a_{i\alpha}^{\dagger} a_{i\beta} \rho] = \mathbb{E}[\langle a_{i}^{\dagger} \rangle_{\rho} \langle a_{i} \rangle_{\rho}] \propto \lambda'_i$, but this expectation must vanish. Hence $\lambda'_i = 0$ and the $\lambda$ matrix is proportional to the identity. We get

$$ \mathcal{H}_{n,baths} = -\sum_{i \in V} \sum_{\alpha=1}^{n} \Gamma_i \left[ \lambda_i (\mathcal{J}_{i,\alpha\alpha}^{+} + C_{i,\alpha\alpha}^{-} - \frac{1}{2}) + (1-\lambda_i)(\mathcal{J}_{i,\alpha\alpha}^{-} - C_{i,\alpha\alpha}^{-} - \frac{1}{2}) \right]. \qquad (39) $$

The set of generators $(\mathcal{J}_{i,\alpha\alpha}^{\pm}, C_{i,\alpha\alpha}^{-})$ build an su(2) algebra for each $1 \le \alpha \le n$ and commute for different $\alpha$: from the point of view of the baths we are back to the single replica problem. We must just repeat the transformations for every replica:

$$ \mathcal{H}' \rightarrow \mathcal{H}_{\text{bulk}} - \sum_{i} \Gamma_{i} \sum_{\alpha} \left\{ \lambda_{i} \mathcal{J}_{i,\alpha\alpha}^{-} - C_{i,\alpha\alpha}^{-} - \frac{1}{2} \right\} \qquad (40) $$

and then continue introducing the empty bath as before, to get:

$$ \mathcal{H}^{\text{dual}} = \mathcal{H}_{\text{bulk}} - \sum_{i\alpha} \Gamma_i \left\{ \mathcal{K}_{i,\alpha\alpha}^{+} \mathcal{J}_{i,\alpha\alpha}^{-} + 2\mathcal{D}_{i,\alpha\alpha}^{-} \mathcal{C}_{i,\alpha\alpha}^{-} - \frac{1}{2} \right\}. \qquad (41) $$

Because of the diagonal nature of the hopping into the bath operator, the bath sites can only exchange creation and destruction operators of the same replica in pairs. Mathematically, this comes from the fact that a diagonal bath respects an extra symmetry, generated by the operators $C_{aa}^{+\text{tot}} = \sum_{i \in V} C_{i,\alpha\alpha}^{+}$ which in fact count the total number of creators minus annihilators
<figure><img src="image_6.png" /><figcaption>Figure 3: The configurations that are relevant in the calculation of two-point correlation functions in the test. Red and black denote replicas one and two, full points are creators and empty points are annihilators.</figcaption></figure>

in the whole system for each species $\alpha$ separately. These numbers are conserved, and equal to zero for all $\alpha$ with reasonable initial conditions, and for reasonable operators $O$ with non-zero expectation. We shall restrict to this situation here. This is shown in Figure 3 where all the possibilities are shown for a two-point operator. We shall then have, for the dual chain at infinite times

$$ \lim_{t \to \infty} \mathrm{Tr}[Oe^{-t\mathcal{H}}\rho(0)] = c_{\{n_{L,\alpha}\}}^{[O]} \lambda_L^{\sum_{\alpha} n_{L,\alpha}} + c_{\{n_{R,\alpha}\}}^{[O]} \lambda_R^{\sum_{\alpha} n_{R,\alpha}}, \quad (42) $$

where we have used the fact that creation and annihilation operators exit in pairs of the same species.

For example, when $O = a_{j,1}^{\dagger}a_{i,1}a_{i,2}^{\dagger}a_{j,2}$ (the case with $n = 2$ replicas), for the stationary state of a chain we have:

$$ \begin{aligned} \mathbb{E}[G_{ij}G_{ji}] &= \mathbb{E}[\mathrm{Tr}(a_{j,1}^{\dagger}a_{i,1}a_{i,2}^{\dagger}a_{j,2} \rho^{(1)} \otimes \rho^{(2)})] \\ &= c_{LL}^{[ij]}\lambda_L^2 + c_{RL}^{[ij]}\lambda_L\lambda_R + c_{LR}^{[ij]}\lambda_R\lambda_L + c_{RR}^{[ij]}\lambda_L^2. \end{aligned} \quad (43) $$

Because the trace is conserved, $c_{LL}^{[ij]} + c_{RL}^{[ij]} + c_{LR}^{[ij]} + c_{RR}^{[ij]} = 0$. In particular, this means that for the stationary equilibrium $\lambda_L = \lambda_R = \lambda$ the expectation $\mathbb{E}[G_{ij}G_{ji}] = 0 \forall i, j$. Similarly, putting $O = a_{i,1}^{\dagger}a_{i,1}a_{j,2}^{\dagger}a_{j,2}$ we have

$$ \mathbb{E}[G_{ii}G_{jj}] = \mathbb{E}[\mathrm{Tr}(a_{i,1}^{\dagger}a_{i,1}a_{j,2}^{\dagger}a_{j,2} \rho^{(1)} \otimes \rho^{(2)}))], \quad (44) $$

and we find, at equilibrium, $\mathbb{E}[G_{ii}G_{jj}] = \lambda^2$, because the trace is now one.

In the non-equilibrium steady state with $\lambda_L \neq \lambda_R$, the expressions for the coefficients appearing in (43) have been calculated in [1]. We conclude this section by mentioning that the invariance of the super-Hamiltonian with respect to unitary rotations of the fermions in replica space may be used to find useful relations for the non-equilibrium stationary state of a chain. For example, in the case with $n = 2$ replicas, the rotation $d_i^{\pm} = (a_{i,1} \pm a_{i,2})/\sqrt{2}$ transforms the expectation (44) and allows us to show that:

$$ \mathbb{E}[G_{ij}G_{ji}] = \mathbb{E}[G_{ii}G_{jj}]_c - \mathbb{E}[G_{ijjj}]_c, \quad (45) $$

where $\mathbb{E}[G_{ijjj}]_c = \langle(d_i^{\dagger})^{\dagger}d_i^{\dagger}(d_j^{\dagger})^{\dagger}d_j^{\dagger}\rangle_c = \langle n_i n_j \rangle_c^{\mathrm{class}}$, the label ‘classical’ denotes the connected correlation of the classical open symmetric exclusion process (see, e.g., Eq(2.4) in [31]). This identity may be verified in the expressions of Reference [1].

## 9 Mapping the driven system into one satisfying detailed balance

Recently it was shown in [19] that a large family of classical diffusive models of transport that has been considered in the past years admits a transformation into the same model in contact
with an equilibrium bath. This mapping holds for any initial conditions and is independent of the topology of the network. Such a construction provides a framework to discuss questions of time reversal in out of equilibrium contexts. Let us now briefly see that this relation holds at the quantum level with the same generality.

Consider a system in which the $\lambda_i = \lambda$ are all the same, thus making it possible to equilibrate. Transform now (40) as

$$ \begin{aligned} \mathcal{H}'' &= e^{-\frac{\lambda}{2}J_{aa}^{-\text{tot}}} \mathcal{H}' e^{\frac{\lambda}{2}J_{aa}^{-\text{tot}}} \\ &= \mathcal{H}_{\text{bulk}} - \sum_{i} \Gamma_{i} \left\{ 2C_{aa}^{-\text{tot}} - \frac{1}{2} \right\} = (\mathcal{H}'')^{\dagger}. \end{aligned} \quad (46) $$

This, when $\mathcal{H}''$ is written in terms of $\mathcal{H}$, is a detailed balance property.

Now, let us go back to a general $\mathcal{H}'$ of (40) and develop the operators on which it acts in subspaces defined by the value $m$ in $C_{aa}^{-\text{tot}}(O) = mO$. Because $[C_{aa}^{-\text{tot}}, \mathcal{J}_{i,aa}^{-}] = -\mathcal{J}_{i,aa}^{-}$ it is clear that the terms in $\mathcal{H}'$ proportional to $\mathcal{J}_{i,aa}^{-}$ (the only ones depending on the $\lambda_i$) are lower ‘triangular’ in this (super) matrix, do not affect the spectrum. Thus, all possible bath combinations are mappable by a similarity transformation into one another, and in particular to the situation with detailed balance (see [19, 32] for an extensive discussion of the classical case).

## 10 Conclusion

We have shown how to construct a dual model for a quantum transport system. We have also shown that there is a transformation that maps all possible evolutions of the system into one in which the bath ensemble satisfies detailed balance. Because these properties are directly deduced from the group structure, the generalization to other models (e.g. the quantum Kac model [33, 34] or quantum KMP [6]) should be immediate – just changing group.

## Acknowledgments

We thank Denis Bernard for helpful discussions and the anonymous referees for their comments. JK is supported by the Simons Foundation Grant No. 454943. RF is supported by the German research foundation (Deutsche Forschungsgemeinschaft DFG) Research Fellowships Programme 416527151. CG is supported by the UniMoRe project ‘ALEA’.

## A Algebra and representations

The generators can be arranged into a $u(2n)$ algebra

$$ \mathcal{E}_{AB} = \begin{pmatrix} C_{\alpha\beta}^{+} + C_{\alpha\beta}^{-} + \frac{\delta_{\alpha\beta}}{2} & \mathcal{J}_{\alpha\beta}^{+} \\ \mathcal{J}_{\alpha\beta}^{-} & C_{\alpha\beta}^{+} - C_{\alpha\beta}^{-} + \frac{\delta_{\alpha\beta}}{2} \end{pmatrix}, \quad (47) $$

satisfying the commutation relations

$$ [\mathcal{E}_{AB}, \mathcal{E}_{CD}] = \delta_{BC}\mathcal{E}_{AD} - \delta_{DA}\mathcal{E}_{CB}. \quad (48) $$

The space of states consists of arbitrary combinations of the fermionic oscillators $(a_\alpha, a_\beta^\dagger)$, with $\alpha, \beta = 1, \dots, n$. From the nilpotency of fermionic oscillators it is then clear that we are dealing with finite-dimensional representations. For a given number of replicas $n$ the space of
states can be decomposed into irreducible representations on which the first Casimir $\sum_{\alpha} C_{\alpha\alpha}^{+}$, which counts the number of creation minus annihilation operators, is proportional to the identity. These reducible representations are the antisymmetric (fundamental) representations of dimension $\binom{2n}{k}$ with Dynkin weights

$$ \lambda^{[k]} = \underbrace{(1, \dots, 1)}_{k} \underbrace{0, \dots, 0)}_{2n-k}, \qquad (49) $$

with $k = 0, 1, \dots, 2n$. The corresponding highest weight states are

$$ v_{hws}^{[k]} = \begin{cases} \prod_{i=1}^{k} a_i^{\dagger} \prod_{j=1}^{n} a_j, & 1 \le k \le n \\ \prod_{i=1}^{n} a_i^{\dagger} \prod_{j=k-n+1}^{n} a_j, & n < k \le 2n \end{cases}. \qquad (50) $$

We verify that $\mathcal{E}_{AB} v_{hws}^{[k]} = 0$ for $1 \le A < B \le 2n$ and $\mathcal{E}_{AA} v_{hws}^{[k]} = \lambda_{A}^{[k]} v_{hws}^{[k]}$ and note that the number of creation minus annihilation operators is related to the index $k$ in (49) via $k-n$.

The Hamiltonian density, cf. (35), is then mapped to the coproduct of the second Casimir

$$ \mathcal{H}_n = - \sum_{A,B=1}^{2n} \mathcal{E}_{AB}^{1} \mathcal{E}_{BA}^{2} + \frac{1}{2} \left( \sum_{A=1}^{2n} \mathcal{E}_{AA}^{1} + \sum_{A=1}^{2n} \mathcal{E}_{AA}^{2} \right). \qquad (51) $$

The first Casimir $\sum_{A=1}^{2n} \mathcal{E}_{AA}$ is local and proportional to the identity for a given irreducible representation.

From the representation theory it is clear that the quantum model can be realised “classically” where $\mathcal{E}_{AB}$ are matrices acting on a vector space. The matrices are block diagonal with block sizes corresponding to the irreducible representations $\lambda^{[k]}$ where $k = 0, 1, \dots, 2n$. Fixing the first Casimir we can restrict to the corresponding block. As described in the main text the boundary terms can be understood in the same framework. The algebra reduces to $u(2)$. It further becomes clear that there is a basis where $\mathcal{E}_{AA}$ is diagonal, $\mathcal{E}_{AB}$ with $A < B$ is upper triangular and $\mathcal{E}_{AB}$ with $A > B$ is lower triangular. This is essential for the proposed duality. We further remark that algebraically the models with $k = 1$ are equivalent to the standard multi-color SSEP, cf. [35], where the generators are replaced by the elementary matrices $\mathcal{E}_{AB} \to E_{AB}$. In the example above, see Figure 3, we have $n = 2$ and $k = 2$ such that the representation is of dimension 6. The matrix realisation of the super operator $\mathcal{H}$ can be obtained by writing the generators in the basis

$$ \begin{aligned} v_1 &= 1, & v_2 &= a_1^\dagger a_1, & v_3 &= a_1^\dagger a_2, \\ v_4 &= a_2^\dagger a_1, & v_5 &= a_2^\dagger a_2, & v_6 &= a_1^\dagger a_1 a_2^\dagger a_2. \end{aligned} $$

## References

[1] D. Bernard and T. Jin, *Open quantum symmetric simple exclusion process*, Phys. Rev. Lett. **123**, 080601 (2019), doi:10.1103/PhysRevLett.123.080601.

[2] F. Spitzer, *Interaction of Markov processes*, Adv. Math. **5**, 246 (1970), doi:10.1016/0001-8708(70)90034-4.

[3] T. M. Liggett, *Interacting particle systems*, Springer Berlin Heidelberg, ISBN 9783540226178 (2005), doi:10.1007/b138374.

[4] H. Spohn, *Long range correlations for stochastic lattice gases in a non-equilibrium steady state*, J. Phys. A: Math. Gen. **16**, 4275 (1983), doi:10.1088/0305-4470/16/18/029.
[5] A. De Masi and E. Presutti, *Mathematical methods for hydrodynamic limits*, Springer Berlin Heidelberg ISBN 9783540550044 (1991), doi:10.1007/BFb0086457.

[6] C. Kipnis, C. Marchioro and E. Presutti, *Heat flow in an exactly solvable model*, J. Stat. Phys. **27**, 65 (1982), doi:10.1007/BF01011740.

[7] C. Giardinà, F. Redig and K. Vafayi, *Correlation inequalities for interacting particle systems with duality*, J. Stat. Phys. **141**, 242 (2010), doi:10.1007/s10955-010-0055-0.

[8] C. Giardinà, J. Kurchan, F. Redig and K. Vafayi, *Duality and hidden symmetries in interacting particle systems*, J. Stat. Phys. **135**, 25 (2009), doi:10.1007/s10955-009-9716-2.

[9] C. Franceschini and C. Giardinà, *Stochastic duality and orthogonal polynomials*, in Springer Proceedings in Mathematics & Statistics, Springer Singapore, ISBN 9789811503016 (2019), doi:10.1007/978-981-15-0302-3_7.

[10] S. Floreani, F. Redig and F. Sau, *Orthogonal polynomial duality of boundary driven particle systems and non-equilibrium correlations*, (2020), arXiv:2007.08272.

[11] C. Franceschini, P. Gonçalves and F. Sau, *Symmetric inclusion process with slow boundary: hydrodynamics and hydrostatics*, (2020), arXiv:2007.11998.

[12] G. M. Schütz, *Duality relations for asymmetric exclusion processes*, J. Stat. Phys. **86**, 1265 (1997), doi:10.1007/BF02183623.

[13] L.-H. Gwa and H. Spohn, *Bethe solution for the dynamical-scaling exponent of the noisy Burgers equation*, Phys. Rev. A **46**, 844 (1992), doi:10.1103/physreva.46.844.

[14] T. Imamura and T. Sasamoto, *Current moments of 1D ASEP by duality*, J. Stat. Phys. **142**, 919 (2011), doi:10.1007/s10955-011-0149-3.

[15] A. Borodin, I. Corwin and T. Sasamoto, *From duality to determinants for q-TASEP and ASEP*, Ann. Probab. **42**, 2314 (2014), doi:10.1214/13-AOP868.

[16] I. Corwin, P. Ghosal, H. Shen and L.-C. Tsai, *Stochastic PDE limit of the six vertex model*, Commun. Math. Phys. **375**, 1945 (2020), doi:10.1007/s00220-019-03678-z.

[17] J. Tailleur, J. Kurchan and V. Lecomte, *Mapping out-of-equilibrium into equilibrium in one-dimensional transport models*, J. Phys. A: Math. Theor. **41**, 505001 (2008), doi:10.1088/1751-8113/41/50/505001.

[18] L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio and C. Landim, *Macroscopic fluctuation theory*, Rev. Mod. Phys. **87**, 593 (2015), doi:10.1103/RevModPhys.87.593.

[19] R. Frassek, C. Giardina and J. Kurchan, *Duality and hidden equilibrium in transport models*, SciPost Phys. **9**, 054 (2020), doi:10.21468/SciPostPhys.9.4.054.

[20] D. A. Rowlands and A. Lamacraft, *Noisy coupled qubits: Operator spreading and the Fredrickson-Andersen model*, Phys. Rev. B **98**, 195125 (2018), doi:10.1103/physrevb.98.195125.

[21] D. A. Rowlands and A. Lamacraft, *Noisy spins and the Richardson-Gaudin model*, Phys. Rev. Lett. **120**, 090401 (2018), doi:10.1103/PhysRevLett.120.090401.

[22] T. Zhou and A. Nahum, *Emergent statistical mechanics of entanglement in random unitary circuits*, Phys. Rev. B **99**, 174205 (2019), doi:10.1103/PhysRevB.99.174205.
[23] M. Bauer, D. Bernard and T. Jin, *Stochastic dissipative quantum spin chains (I) : Quantum fluctuating discrete hydrodynamics*, SciPost Phys. **3**, 033 (2017), doi:10.21468/SciPostPhys.3.5.033.

[24] M. Bauer, D. Bernard and T. Jin, *Equilibrium fluctuations in maximally noisy extended quantum systems*, SciPost Phys. **6**, 045 (2019), doi:10.21468/SciPostPhys.6.4.045.

[25] G. Schütz and S. Sandow, *Non-Abelian symmetries of stochastic processes: Derivation of correlation functions for random-vertex models and disordered-interacting-particle systems*, Phys. Rev. E **49**, 2726 (1994), doi:10.1103/PhysRevE.49.2726.

[26] C. Giardinà, J. Kurchan and F. Redig, *Duality and exact correlations for a model of heat conduction*, J. Math. Phys. **48**, 033301 (2007), doi:10.1063/1.2711373.

[27] F. H. L. Essler, H. Frahm, F. Göhmann, A. Klümper and V. E. Korepin, *The one-dimensional Hubbard model*, Cambridge University Press, Cambridge, United Kingdom, ISBN 9780521802628 (2005), doi:10.1017/CBO9780511534843.

[28] E. Ilievski, J. De Nardis, B. Wouters, J.-S. Caux, F. H. L. Essler and T. Prosen, *Complete generalized Gibbs ensembles in an interacting theory*, Phys. Rev. Lett. **115**, 157201 (2015), doi:10.1103/PhysRevLett.115.157201.

[29] J.-S. Caux and F. H. L. Essler, *Time evolution of local observables after quenching to an integrable model*, Phys. Rev. Lett. **110**, 257203 (2013), doi:10.1103/PhysRevLett.110.257203.

[30] P. Calabrese, F. H. L. Essler and M. Fagotti, *Quantum quench in the transverse-field ising chain*, Phys. Rev. Lett. **106**, 227203 (2011), doi:10.1103/PhysRevLett.106.227203.

[31] B. Derrida, J. L. Lebowitz and E. R. Speer, *Entropy of open lattice systems*, J. Stat. Phys. **126**, 1083 (2007), doi:10.1007/s10955-006-9160-5.

[32] R. Frassek, *Eigenstates of triangularisable open XXX spin chains and closed-form solutions for the steady state of the open SSEP*, J. Stat. Mech. 053104 (2020), doi:10.1088/1742-5468/ab7af3.

[33] M. Kac, *Foundations of kinetic theory*, In Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability **3**, 171 (1956).

[34] E. A. Carlen, M. C. Carvalho and M. P. Loss, *Chaos, ergodicity, and equilibria in a quantum Kac model*, Adv. Math. **358**, 106827 (2019), doi:10.1016/j.aim.2019.106827.

[35] C. Arita, A. Kuniba, K. Sakai and T. Sawabe, *Spectrum of a multi-species asymmetric simple exclusion process on a ring*, J. Phys. A: Math. Theor. **42**, 345002 (2009), doi:10.1088/1751-8113/42/34/345002.