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