Zhi-Yuan Wei, Joel Rajakumar, Jon Nelson, Daniel Malz, Michael J. Gullans, Alexey V. Gorshkov (Mar 24 2026).
Abstract: We study how matrix-product-operator (MPO) truncation errors evolve when simulating two setups: (1) 1D Haar-random circuits under either depolarizing noise or amplitude-damping noise, and (2) 1D Lindbladian dynamics of a non-integrable quantum Ising model under either depolarizing or amplitude-damping noise. We first show that the average purity of the system density matrix relaxes to a steady value on a timescale that scales inversely with the noise rate. We then show that truncation errors contract exponentially in both system size
N and the evolution time
t, as the noisy dynamics maps different density matrices toward the same steady state. This yields an empirical bound on the
L1​ truncation error that is exponentially tighter in
N than the existing bound. Together, these results provide empirical evidence that MPO simulation algorithms may efficiently sample from the output of 1D noisy random circuits [setup (1)] at arbitrary circuit depth, and from the steady state of 1D Lindbladian dynamics [setup (2)].