Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2004.05177

Dissipation-assisted operator evolution method for capturing hydrodynamic transport

Tibor Rakovszky, Curt W. von Keyserlingk, Frank Pollmann

Published Apr 10, 2020• 17 citations
arXivHTMLPDF

Abstract

We introduce the dissipation-assisted operator evolution (DAOE) method for calculating transport properties of strongly interacting lattice systems in the high temperature regime. DAOE is based on evolving observables in the Heisenberg picture, and applying an artificial dissipation that reduces the weight on non-local operators. We represent the observable as a matrix product operator, and show that the dissipation leads to a decay of operator entanglement, allowing us to capture the dynamics to long times. We test this scheme by calculating spin and energy diffusion constants in a variety of physical models. By gradually weakening the dissipation, we are able to consistently extrapolate our results to the case of zero dissipation, thus estimating the physical diffusion constant with high precision.

DOI
10.48550/arxiv.2004.05177

Related papers

Energy diffusion in weakly interacting chains with fermionic dissipation assisted operator evolution

En-Jui Kuo, Brayden Ware, Peter Lunts +2 more

Interacting lattice Hamiltonians at high temperature generically give rise to energy transport governed by the classical diffusion equation; however, predicting the rate of diffusion requires numerical simulation of the microscopic quantum dynamics. For the purpose of predicting such transport properties, computational time evolution methods must be paired with schemes to control the growth of entanglement to tractably simulate for sufficiently long times. One such truncation scheme---dissipation-assisted operator evolution (DAOE)---controls entanglement by damping out components of operators with large Pauli weight. In this paper, we generalize DAOE to treat fermionic systems. Our method instead damps out components of operators with large fermionic weight. We investigate the performance of DAOE, the new fermionic DAOE (FDAOE), and another simulation method, density matrix truncation (DMT), in simulating energy transport in an interacting one-dimensional Majorana chain. The chain is found to have a diffusion coefficient scaling like interaction strength to the fourth power, contrary to naive expectations based on Fermi's Golden rule---but consistent with recent predictions based on the theory of 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.

arXiv:2311.17148•Physical review. B./Physical review. B•Aug 26, 2024•6 citations

Dissipation-assisted operator evolution method for capturing hydrodynamic transport

Tibor Rakovszky, Curt W. von Keyserlingk, Frank Pollmann

Transport properties are among the most practically important features of materials, yet their precise calculation in strongly interacting systems poses a formidable challenge. One might even expect that it requires a quantum computer. The authors introduce an algorithm that tackles this challenge with only moderate resources on a classical computer. It relies on introducing an artificial dissipation process that only affects nonlocal observables, thus reducing the computational cost without strongly influencing physically relevant quantities. The algorithm is tested on various paradigmatic models with good convergence results.

Physical review. B./Physical review. B•Feb 16, 2022•105 citations

The ballistic to diffusive crossover in a weakly-interacting Fermi gas

Jerome Lloyd, Tibor Rakovszky, Frank Pollmann +1 more

Charge and energy are expected to diffuse in interacting systems of fermions at finite temperatures, even in the absence of disorder, with the interactions inducing a crossover from the coherent and ballistic streaming of quasi-particles at early times, to incoherent diffusive behavior at late times. The relevant crossover timescales and the transport coefficients are both controlled by the strength of interactions. In this work we develop a numerical method to simulate such systems at high temperatures, applicable in a wide range of interaction strengths, by adapting Dissipation-assisted Operator Evolution (DAOE) to fermions. Our fermion DAOE, which approximates the exact dynamics by systematically discarding information from high $n$-point functions, is tailored to capture non-interacting dynamics exactly, thus providing a good starting point for the weakly interacting problem. Applying our method to a microscopic model of weakly interacting fermions, we numerically demonstrate that the crossover from ballistic to diffusive transport happens at a time $t_D\sim1/Δ^{2}$ and that the diffusion constant similarly scales as $D \sim 1/Δ^2$, where $Δ$ is the interaction strength. We substantiate this scaling with a Fermi's golden rule calculation in the operator spreading picture, interpreting $t_D$ as the fermion-fermion scattering time and lifetime of the single-particle Green's function.

arXiv:2310.16043•Oct 24, 2023•0 citations

Probing hydrodynamic crossovers with dissipation-assisted operator evolution

N. S. Srivatsa, Oliver Lunt, Tibor Rakovszky +1 more

Using artificial dissipation to tame entanglement growth, we chart the emergence of diffusion in a generic interacting lattice model for varying U(1) charge densities. We follow the crossover from ballistic to diffusive transport above a scale set by the scattering length, finding the intuitive result that the diffusion constant scales as D ∝ 1 / ρ at low densities ρ . Our numerical approach generalizes the Dissipation-Assisted Operator Evolution algorithm: in the spirit of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, we effectively approximate nonlocal operators by their ensemble averages, rather than discarding them entirely. This greatly reduces the operator entanglement entropy, while still giving accurate predictions for diffusion constants across all density scales. We further construct a minimal model for the transport crossover, yielding charge correlation functions which agree well with our numerical data. Our results clarify the dominant contributions to hydrodynamic correlation functions of conserved densities, and serve as a guide for generalizations to low-temperature transport.

arXiv:2408.08249•Physical review. B./Physical review. B•Apr 27, 2026•0 citations

A Universal Approach to Hydrodynamics in Higher Dimensions via Tensor Networks

Sajant Anand, Christopher D. White, Jack Kemp +3 more

Extracting late-time hydrodynamic transport properties directly from microscopic quantum models, especially in dimensions beyond one, is believed to be classically intractable. This owes, in part, to the extensive growth of entanglement in higher dimensions and often a mismatch between geometric and tensor network locality. Here, we extend the Density Matrix Truncation (DMT) algorithm to accurately and quantitatively predict high-temperature transport coefficients (e.g. diffusion constants) in up to three dimensions. Additionally, we demonstrate the ability to capture the slow crossover dynamics characteristic of weakly broken integrability. The intuition underlying our approach is to preserve the geometrically local information crucial for hydrodynamics throughout the time evolution, despite using a one-dimensional matrix product state ansatz for higher-dimensional calculations. Our work pushes the boundary of classical computability and opens the door to benchmarking the high-temperature, late-time dynamics of generic quantum simulators.

quant-pharXiv:2610.40776•Nov 2, 2026

Open-system spin transport and operator weight dissipation in spin chains

Yongchan Yoo, Christopher David White, Brian G. Swingle

We use nonequilibrium steady states to study the effect of dissipation-assisted operator evolution (DAOE) on the scaling behavior of transport in one-dimensional spin chains. We consider three models in the XXZ family: the XXZ model with staggered anisotropy, which is chaotic; XXZ model with no external field and tunable interaction, which is Bethe-ansatz integrable and (in the zero interaction limit) free-fermion integrable; and the disordered XY model, which is free-fermion integrable and Anderson localized. We find evidence that DAOE's effect on transport is controlled by its effect on the system's conserved quantities. To the extent that DAOE preserves those symmetries, it preserves the scaling of the system's transport properties; to the extent it breaks those conserved quantities, it pushes the system towards diffusive scaling of transport.

arXiv:2210.06494•Physical review. B./Physical review. B•Mar 9, 2023•7 citations
We introduce the dissipation-assisted operator evolution (DAOE) method for calculating transport properties of strongly interacting lattice systems in the high temperature regime. DAOE is based on evolving observables in the Heisenberg picture, and applying an artificial dissipation that reduces the weight on non-local operators. We represent the observable as a matrix product operator, and show that the dissipation leads to a decay of operator entanglement, allowing us to capture the dynamics to long times. We test this scheme by calculating spin and energy diffusion constants in a variety of physical models. By gradually weakening the dissipation, we are able to consistently extrapolate our results to the case of zero dissipation, thus estimating the physical diffusion constant with high precision.
Interacting lattice Hamiltonians at high temperature generically give rise to energy transport governed by the classical diffusion equation; however, predicting the rate of diffusion requires numerical simulation of the microscopic quantum dynamics. For the purpose of predicting such transport properties, computational time evolution methods must be paired with schemes to control the growth of entanglement to tractably simulate for sufficiently long times. One such truncation scheme---dissipation-assisted operator evolution (DAOE)---controls entanglement by damping out components of operators with large Pauli weight. In this paper, we generalize DAOE to treat fermionic systems. Our method instead damps out components of operators with large fermionic weight. We investigate the performance of DAOE, the new fermionic DAOE (FDAOE), and another simulation method, density matrix truncation (DMT), in simulating energy transport in an interacting one-dimensional Majorana chain. The chain is found to have a diffusion coefficient scaling like interaction strength to the fourth power, contrary to naive expectations based on Fermi's Golden rule---but consistent with recent predictions based on the theory of 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.
Transport properties are among the most practically important features of materials, yet their precise calculation in strongly interacting systems poses a formidable challenge. One might even expect that it requires a quantum computer. The authors introduce an algorithm that tackles this challenge with only moderate resources on a classical computer. It relies on introducing an artificial dissipation process that only affects nonlocal observables, thus reducing the computational cost without strongly influencing physically relevant quantities. The algorithm is tested on various paradigmatic models with good convergence results.
Charge and energy are expected to diffuse in interacting systems of fermions at finite temperatures, even in the absence of disorder, with the interactions inducing a crossover from the coherent and ballistic streaming of quasi-particles at early times, to incoherent diffusive behavior at late times. The relevant crossover timescales and the transport coefficients are both controlled by the strength of interactions. In this work we develop a numerical method to simulate such systems at high temperatures, applicable in a wide range of interaction strengths, by adapting Dissipation-assisted Operator Evolution (DAOE) to fermions. Our fermion DAOE, which approximates the exact dynamics by systematically discarding information from high nnn-point functions, is tailored to capture non-interacting dynamics exactly, thus providing a good starting point for the weakly interacting problem. Applying our method to a microscopic model of weakly interacting fermions, we numerically demonstrate that the crossover from ballistic to diffusive transport happens at a time tD∼1/Δ2t_D\sim1/Δ^{2}tD​∼1/Δ2 and that the diffusion constant similarly scales as D∼1/Δ2D \sim 1/Δ^2D∼1/Δ2, where ΔΔΔ is the interaction strength. We substantiate this scaling with a Fermi's golden rule calculation in the operator spreading picture, interpreting tDt_DtD​ as the fermion-fermion scattering time and lifetime of the single-particle Green's function.
Using artificial dissipation to tame entanglement growth, we chart the emergence of diffusion in a generic interacting lattice model for varying U(1) charge densities. We follow the crossover from ballistic to diffusive transport above a scale set by the scattering length, finding the intuitive result that the diffusion constant scales as D ∝ 1 / ρ at low densities ρ . Our numerical approach generalizes the Dissipation-Assisted Operator Evolution algorithm: in the spirit of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, we effectively approximate nonlocal operators by their ensemble averages, rather than discarding them entirely. This greatly reduces the operator entanglement entropy, while still giving accurate predictions for diffusion constants across all density scales. We further construct a minimal model for the transport crossover, yielding charge correlation functions which agree well with our numerical data. Our results clarify the dominant contributions to hydrodynamic correlation functions of conserved densities, and serve as a guide for generalizations to low-temperature transport.
Extracting late-time hydrodynamic transport properties directly from microscopic quantum models, especially in dimensions beyond one, is believed to be classically intractable. This owes, in part, to the extensive growth of entanglement in higher dimensions and often a mismatch between geometric and tensor network locality. Here, we extend the Density Matrix Truncation (DMT) algorithm to accurately and quantitatively predict high-temperature transport coefficients (e.g. diffusion constants) in up to three dimensions. Additionally, we demonstrate the ability to capture the slow crossover dynamics characteristic of weakly broken integrability. The intuition underlying our approach is to preserve the geometrically local information crucial for hydrodynamics throughout the time evolution, despite using a one-dimensional matrix product state ansatz for higher-dimensional calculations. Our work pushes the boundary of classical computability and opens the door to benchmarking the high-temperature, late-time dynamics of generic quantum simulators.
We use nonequilibrium steady states to study the effect of dissipation-assisted operator evolution (DAOE) on the scaling behavior of transport in one-dimensional spin chains. We consider three models in the XXZ family: the XXZ model with staggered anisotropy, which is chaotic; XXZ model with no external field and tunable interaction, which is Bethe-ansatz integrable and (in the zero interaction limit) free-fermion integrable; and the disordered XY model, which is free-fermion integrable and Anderson localized. We find evidence that DAOE's effect on transport is controlled by its effect on the system's conserved quantities. To the extent that DAOE preserves those symmetries, it preserves the scaling of the system's transport properties; to the extent it breaks those conserved quantities, it pushes the system towards diffusive scaling of transport.