---
title: Dissipation-Assisted Operator Evolution
url: https://www.emergentmind.com/topics/dissipation-assisted-operator-evolution-daoe
type: topic
---

# Dissipation-Assisted Operator Evolution

Dissipation-Assisted Operator Evolution (DAOE) is a Heisenberg-picture method for calculating transport properties of strongly interacting lattice systems in the high temperature regime. The method modifies unitary operator evolution by inserting an artificial dissipation that reduces the weight on non-local operators, typically defined in a Pauli-string or fermionic-string basis. When the evolving observable is represented as a matrix product operator (MPO), this dissipation causes operator entanglement to turn over and decay at late times, thereby extending the accessible simulation window. Physical transport coefficients are then obtained by weakening the dissipation and extrapolating to the zero-dissipation limit [2004.05177]. Subsequent work generalized the approach to finite chemical potential, weakly interacting fermions, open-system nonequilibrium steady states, and quasiparticle decay problems in one dimension [2408.08249; 2310.16043; 2210.06494; 2508.00254].

## 1. Formal definition and operator-weight dissipation

In its original form, DAOE evolves an operator \(O(t)\) under a Hamiltonian \(H\) that conserves some charge \(Q\), but replaces exact Heisenberg evolution by
\[
\frac{d}{dt}O(t)=i[H,O(t)]-\gamma\,\mathcal D[O(t)] ,
\]
where \(\gamma>0\) is a dissipation rate and \(\mathcal D\) acts diagonally on a local operator basis such as Pauli strings. If \(S\) is a Pauli string of length \(\ell_S\), with \(\ell_S\) defined as the number of non-identity factors, and if \(\ell_*\) is a cutoff length, then
\[
\mathcal D[S]=
\begin{cases}
0,&\ell_S\le \ell_*\\
(\ell_S-\ell_*)S,&\ell_S>\ell_* \, .
\end{cases}
\]
Equivalently, strings with \(\ell_S>\ell_*\) are multiplied by \(e^{-\gamma(\ell_S-\ell_*)}\), while strings with \(\ell_S\le \ell_*\) are left unchanged. In the limit \(\gamma\to0\) or \(\ell_*\to\infty\), one recovers exact unitary evolution [2004.05177].

A discrete-time formulation alternates unitary evolution and dissipative pruning. Writing the Liouvillian as \(\mathcal L(O)=[H,O]\), one applies \(e^{i\mathcal L\Delta t}\) and then the diagonal superoperator \(e^{-\gamma\Delta t\,\mathcal D}\). This realizes DAOE as a Trotterized evolution in operator space, with the key approximation located entirely in the suppression of high-weight strings rather than in a modification of the microscopic Hamiltonian [2004.05177].

The notion of “weight” is central. In spin chains it is the number of nontrivial Pauli factors in a tensor-product string, \(\ell_{\mathcal S}=\#\{j\mid \alpha_j\neq0\}\) for \(\mathcal S=\bigotimes_j \sigma_j^{\alpha_j}\). In fermionic generalizations, the relevant measure is instead a “Fermi weight,” constructed from a local basis \(o^0=I\), \(o^1=f^\dagger\), \(o^2=f\), \(o^3=[f^\dagger,f]=2f^\dagger f-I\), with weights \(0,1,1,2\) respectively [2310.16043].

The open-system literature emphasizes a closely related point: DAOE is an artificial, linear superoperator acting on operator Hilbert space, and while it is contractive and Hermitian under the Hilbert–Schmidt inner product, it does not in general preserve positivity of density matrices. In boundary-driven Lindblad problems, it therefore functions as a controlled truncation device rather than as a physical bath in the GKLS sense [2210.06494].

## 2. MPO representation, vectorization, and operator-entanglement control

DAOE is designed for tensor-network simulation. One reshapes an operator \(O(t)\) into a vectorized state \(|O(t)\rangle\) in a doubled Hilbert space and represents it as a matrix-product state (MPS) of bond dimension \(\chi\). In this language, the dissipation superoperator \(e^{-\gamma\mathcal D}\) itself becomes an MPO of bond dimension \(\ell_*+1\), because it only needs to count operator weight up to the cutoff and apply the appropriate multiplicative factor [2004.05177].

Without dissipation, generic Hamiltonian evolution produces linear growth of operator entanglement,
\[
S_{\rm ent}(t)=-\mathrm{Tr}\bigl[\rho_A(t)\log\rho_A(t)\bigr]\sim \alpha t ,
\]
where \(\rho_A(t)\) is the reduced density matrix of half the chain in the vectorized state \(|O(t)\rangle\). Under DAOE with \(\gamma>0\) and finite \(\ell_*\), \(S_{\rm ent}(t)\) initially rises but then turns over and decays at late times. This is the mechanism that permits long-time evolution with fixed finite \(\chi\): the method suppresses the operator strings responsible for unbounded entanglement growth while preserving the low-weight sectors that dominate hydrodynamic observables [2004.05177].

A practical implementation proceeds by choosing \(\ell_*\), \(\gamma\), \(\chi\), and a Trotter step \(\Delta t\); initializing \(|O(0)\rangle\) as an MPS; applying one TEBD step for \(e^{i\mathcal L\Delta t}\); applying the MPO for \(e^{-\gamma\mathcal D\Delta t}\); and monitoring truncation errors at each gate, increasing \(\chi\) until a preset tolerance is met. Typical parameter regimes in the original benchmarks are \(\Delta t\sim0.01\text{--}1\), \(\ell_*=2\text{--}5\), and \(\chi=256\text{--}768\), which allow times \(t\gtrsim20\) on system sizes \(L\gtrsim50\) with negligible finite-size or truncation effects. The computational cost scales as \(O(N_{\rm steps}\times L\times\chi^3)\) [2004.05177].

Later tensor-network studies reported the same qualitative gain. In Majorana-chain simulations, DAOE converts a strictly growing bond dimension into a profile that rises and then relaxes back, allowing times \(\sim200\) or more in 1D chains with modest \(\chi\sim64\text{--}256\). In finite-density transport studies, the generalized method reaches \(t\gg1/\rho\) with moderate \(\chi\lesssim10^3\), precisely because the dissipator strongly biases the entanglement spectrum toward low-weight sectors [2311.17148; 2408.08249].

## 3. Transport diagnostics and zero-dissipation extrapolation

DAOE was introduced primarily as a transport method. For a local conserved density \(q_j\), one computes the infinite-temperature correlation function
\[
C_j(t)=\frac{1}{\mathcal N}\mathrm{Tr}\bigl[q_j\,\tilde q_{j_0}(t)\bigr],
\]
the mean-square displacement
\[
d^2(t)=\sum_j (j-j_0)^2\,C_j(t),
\]
and the time-dependent diffusion constant
\[
D(t)=\tfrac12\,\frac{d}{dt}d^2(t).
\]
For each fixed \(\gamma>0\), \(D(t)\) approaches a late-time value \(D(\gamma)\). As \(\gamma\) is reduced, the numerical data exhibit a linear convergence
\[
D(\gamma)\simeq D_0+A\gamma \qquad (\gamma\ll1),
\]
so the physical diffusion constant is extracted as \(D_0=D(\gamma\to0)\) by fitting the last few points in \(D(\gamma)\) versus \(\gamma\). Repeating this for several \(\ell_*\) yields an estimate of systematic uncertainty through the spread in \(D_0\), typically of order a few percent [2004.05177].

The original benchmarks established the method on two generic 1D models at infinite temperature. For the tilted-field Ising chain,
\[
H=\sum_j\bigl(g_xX_j+g_zZ_j+(Z_{j-1}Z_j+Z_jZ_{j+1})/2\bigr),
\qquad g_x=1.4,\ g_z=0.9045,
\]
with conserved density \(q_j=h_j\), the extrapolated energy diffusion constant is
\[
D_{\rm energy}=1.40\pm0.02.
\]
For the two-leg XX spin ladder,
\[
H=\sum_{j,a}[X_{j,a}X_{j+1,a}+Y_{j,a}Y_{j+1,a}]+\sum_j(X_{j,1}X_{j,2}+Y_{j,1}Y_{j,2}),
\]
with conserved density \(q_j=(Z_{j,1}+Z_{j,2})/2\), the extrapolated spin diffusion constant is
\[
D_{\rm spin}=0.95\pm0.03,
\]
in excellent agreement with prior estimates [2004.05177].

An alternative transport diagnostic appears in boundary-driven open systems. There one computes the nonequilibrium steady-state current
\[
J\propto \frac{\Delta S}{N^\chi},
\]
with \(\chi=0\) for ballistic transport, \(0<\chi<1\) for superdiffusion, \(\chi=1\) for diffusion, and \(\chi>1\) for subdiffusion. This NESS formulation does not replace the MSD analysis, but it provides a complementary lens on how DAOE affects transport scaling through its action on conserved operators [2210.06494].

## 4. Finite chemical potential, BBGKY-style truncation, and hydrodynamic crossover

The \(\mu\neq0\) generalization modifies the basis in which operator weight is defined. At finite chemical potential one introduces the inner product
\[
(A|B)_\mu \equiv \mathrm{Tr}[A^\dagger \rho_\mu^{1/2} B \rho_\mu^{1/2}],\qquad
\rho_\mu=e^{-\mu\sum_i \sigma_i^z}/Z,
\]
and Gram–Schmidt orthonormalizes the local Pauli basis to \(\{\mathbb 1,\tilde\sigma^z,\tilde\sigma^x,\tilde\sigma^y\}\), with
\[
\tilde\sigma^z\equiv c_\mu(\sigma^z-\langle \sigma^z\rangle_\mu\,\mathbb 1),\qquad
\tilde\sigma^{x,y}=\sqrt{c_\mu}\,\sigma^{x,y},\qquad
c_\mu\equiv \cosh(\mu/2).
\]
The same cutoff-and-damping rule is then applied in this orthonormalized basis, through
\[
\mathcal D_{\ell_*,\gamma}^\mu[O]\equiv M_\mu\,\mathcal D_{\ell_*,\gamma}[M_\mu^{-1}O],
\]
with \(M_\mu\) the local basis change. As in the original construction, the resulting dissipator is encoded as an MPO of bond dimension \(\ell_*+1\) [2408.08249].

The rationale changes in an important way. At \(\mu=0\), the dominant Feynman paths in \(C_{ij}(t)=\mathrm{Tr}[\sigma_i^z(t)\sigma_j^z]/\mathrm{Tr}[1]\) start and end on low-weight operators, and ergodicity suggests that once a string grows beyond \(\ell_*\) it is unlikely to return. At finite \(\mu\), long \(\sigma^z\)-type strings carry nonzero overlap with \(\sigma^z\) and cannot simply be discarded. In the \(\tilde\sigma\) basis, throwing away a long \(\tilde\sigma\)-string is equivalent to replacing each \(\tilde\sigma^z\to\langle \sigma^z\rangle_\mu\,\mathbb 1\), so high-order strings are replaced by disconnected ensemble averages in the spirit of the BBGKY hierarchy. This preserves the long-string physics relevant for density correlations while still dramatically reducing operator entanglement [2408.08249].

This finite-density generalization resolves the ballistic-to-diffusive crossover at low filling. Defining the equilibrium density
\[
\rho(\mu)=\tfrac12[1-\tanh(\mu/2)],
\]
one finds numerically
\[
D(\mu)\simeq \frac{c_\mu^2}{2\lambda}=\frac{\cosh^2(\mu/2)}{2\lambda},
\]
so for \(\mu\gg1\), \(D\sim e^\mu\sim1/\rho\). The same scaling follows from a memory-matrix estimate with fast-mode decay rate \(\Gamma=\lambda c_\mu^{-2}\sim \lambda e^{-\mu}\), which gives \(D=(2\Gamma)^{-1}\). The interpretation is that low-density transport has a long ballistic regime with mean free path or scattering length \(\ell_{\rm sc}\sim1/\rho\) before crossing over to diffusion [2408.08249].

The real-space correlator exhibits this crossover explicitly. For distances and times \(\ell,t\ll \ell_{\rm sc}\approx1/\rho\), one has the free-particle Bessel form
\[
C(r,t)\approx J_{|r|}(t)^2,
\]
whereas for \(t\gg1/\rho\), the Fourier-space correlator becomes diffusive,
\[
C(k,t)\approx e^{-Dk^2 t}.
\]
A minimal memory-matrix model recovers the same structure by separating slow operators from an \(\ell\ge3\) fast subspace, approximating the fast autocorrelator as white noise, and deriving a reduced Liouvillian whose poles describe a damped ballistic mode and a diffusive contribution. A single extra fit parameter \(\alpha\) introduced through \(\lambda(k)=\lambda_0+\alpha\sin^2(k/2)\) quantitatively reproduces the DAOE\(^\mu\) data for all \(k\lesssim\pi/2\) [2408.08249].

## 5. Fermionic DAOE, exact free limits, and weak-integrability-breaking regimes

The fermionic formulation adapts DAOE to problems where Pauli weight is not the physically natural notion of operator size. In one dimension, fermions are mapped to spins via Jordan–Wigner,
\[
f_x=\Bigl(\prod_{r<x}\sigma_r^z\Bigr)\sigma_x^-,
\qquad
f_x^\dagger=\Bigl(\prod_{r<x}\sigma_r^z\Bigr)\sigma_x^+,
\]
and the dissipator is built in a Fermi-string basis. The fermionic superoperator acts diagonally as
\[
\mathcal G_{\ell_*,\gamma}|o^\mu\rangle
=e^{-\gamma\,\max(0,\ell(\mu)-\ell_*)}|o^\mu\rangle ,
\]
so strings with Fermi weight above \(\ell_*\) are exponentially suppressed [2310.16043].

A decisive property of fDAOE is exact recovery of free dynamics. For a quadratic Hamiltonian, any single-body operator \(f_x\) or \(f_x^\dagger\) evolves only within the single-body sector, and any density \(n_x(t)=f_x^\dagger(t)f_x(t)\) remains weight-2. If \(\ell_*\ge2\), the dissipator never acts on this sector. Consequently, fDAOE reproduces ballistic spreading with no approximation error in the free limit and provides an ideal perturbative starting point when interactions \(\Delta\) are small [2310.16043].

Algorithmically, the fermionic dissipator is implemented as a finite-state automaton MPO of bond dimension \(\ell_*+1\), split into parity-even and parity-odd blocks. The automaton tracks both the current Fermi weight and the fermion parity, so that the Jordan–Wigner \(\sigma^z\) strings required by anticommutation do not artificially count toward Fermi weight when they lie between \(f\) and \(f^\dagger\). Time evolution then alternates Trotterized two-site gates with applications of \(\mathcal G_{\ell_*,\gamma}\), followed by standard MPS truncation [2310.16043].

In a weakly interacting Fermi gas, this framework reveals a ballistic-to-diffusive crossover controlled by interaction strength. Numerically one finds \(D(t)\sim t\) at early times, saturation at \(t_D\sim1/\Delta^2\), and a late-time diffusion constant \(D\sim1/\Delta^2\). A memory-kernel derivation and a Fermi’s golden rule estimate identify \(t_D\) as the fermion–fermion scattering time and the lifetime of the single-particle Green’s function [2310.16043].

In an interacting one-dimensional Majorana chain,
\[
H=\sum_n i\,\eta_n\eta_{n+1}-U\sum_n\eta_{n-1}\eta_n\eta_{n+1}\eta_{n+2},
\]
energy transport instead yields
\[
D_\infty\propto U^{-4}\qquad (U\ll1),
\]
contrary to naive expectations based on Fermi’s Golden rule but consistent with recent predictions based on weak integrability breaking. In the weak-interaction regime, where the fermionic nature of the system is most relevant, FDAOE is found to simulate the system more efficiently than DAOE. The same work contrasts DAOE/FDAOE with density matrix truncation (DMT): DMT preserves all operators up to a support diameter \(\ell_{\rm pres}\) but has no small parameter like \(\gamma\), whereas DAOE/FDAOE provide a perturbative control parameter and a clear operator-weight criterion [2311.17148].

DAOE has also been used in a more delicate 1D finite-temperature setting where naive Fermi’s Golden Rule breaks down. Combining DAOE with particle–particle ladder resummations and a melonic memory-matrix resummation, one study predicts a logarithmic enhancement of the quasiparticle decay rate,
\[
\tau^{-1}\sim \Delta^2\log \Delta^{-2},
\]
for interacting 1D lattice fermions at non-zero temperature. The DAOE simulations in that work use system sizes up to \(L=600\), bond dimension \(\chi\le120\), cutoff \(l^*=3\), dissipation period \(\delta t_\gamma=1.0\), and \(\gamma\in[0.005,0.03]\), with linear extrapolation in \(\gamma\to0\) [2508.00254].

## 6. Conserved quantities, limitations, and nomenclature

A recurring conclusion across applications is that DAOE’s effect on transport is controlled by its effect on the system’s conserved quantities. In boundary-driven XXZ-family chains and in the disordered XY model, the scaling exponent \(\chi\) is preserved when DAOE preserves the relevant local currents and is driven toward diffusion when DAOE breaks the operators responsible for ballistic, superdiffusive, or localized behavior [2210.06494].

The model dependence is explicit. In the chaotic staggered-anisotropy XXZ chain, DAOE preserves the total \(S^z\) and the two-site spin current for \(\ell_*\ge2\), and the system remains diffusive with \(\chi\approx1\). In clean XXZ at \(\Delta<1\), ballistic transport survives for sufficiently small \(\gamma\) when \(\ell_*>1\), because the total current operator is itself a sum of weight-2 Pauli strings; if \(\ell_*=1\), DAOE damps the current and the system becomes diffusive. At the \(SU(2)\)-symmetric point \(\Delta=1\), \(\ell_*=1\) restores generic diffusion, while \(\ell_*\ge2\) leaves a remnant superdiffusion on simulated sizes \(N\sim256\), with an expected crossover to \(\chi=1\) only for \(N\gg\ell_*/\gamma\). In the disordered XY chain, DAOE acts as a dephasing bath for localized orbitals and produces a crossover to diffusive scaling at lengths \(N\gg \ell_*/\xi\) [2210.06494].

These results clarify a frequent misconception. DAOE is not a generic “noise” that necessarily forces diffusion; rather, it selectively damps high-weight sectors, and the hydrodynamic outcome depends on whether the physically relevant conserved densities and currents remain inside the retained low-weight subspace. A plausible implication is that operator-weight truncation is best interpreted as a symmetry-sensitive coarse graining rather than as a uniform decoherence mechanism [2210.06494].

The method also has clear limitations. It is best suited to high temperatures; at low \(T\), the thermal density matrix itself has growing operator entanglement. The parameters \(\ell_*\) and \(\Delta t\) must be chosen so that the cutoff does not project away hydrodynamic operators, although the final answer is insensitive to their precise value as long as \(\ell_*\) exceeds the support of the conserved density and \(\Delta t\) is of order the microscopic timescale. If nonlocal strings carry important physics, DAOE must be pushed to larger \(\ell_*\) or smaller \(\gamma\), increasing cost. In finite-density memory-matrix treatments, the white-noise approximation for the fast sector is uncontrolled [2004.05177; 2408.08249; 2508.00254].

The acronym “DAOE” also has a separate, unrelated usage on arXiv. “Energy-Dissipative Evolutionary Deep Operator Neural Networks” uses DAOE to denote a DeepONet-based operator-learning framework for gradient-flow-type PDEs with a scalar auxiliary variable construction enforcing unconditional energy dissipation at the discrete level. That nomenclature refers to a neural-operator method for dissipative PDEs, not to Heisenberg-picture operator-weight truncation in quantum lattice dynamics [2306.06281].

Source: https://www.emergentmind.com/topics/dissipation-assisted-operator-evolution-daoe