We introduce a general numerical method to compute the dynamics and multitime correlations of chains of quantum systems, where each system may couple strongly to a structured environment. The method combines the process tensor formalism for general (possibly non-Markovian) open quantum systems with time-evolving block decimation for one-dimensional chains. It systematically reduces the numerical complexity originating from system-environment correlations before integrating them into the full many-body problem, making a wide range of applications numerically feasible. We illustrate the power of this method by studying two examples. First, we study the thermalization of individual spins of a short XYZ Heisenberg chain with strongly coupled thermal leads. Our results confirm the complete thermalization of the chain when coupled to a single bath, and they reveal distinct effective temperatures in low-, mid-, and high-frequency regimes when the chain is placed between a hot and a cold bath. Second, we study the dynamics of diffusion in a longer XY chain, when each site couples to its own bath.

• Received 14 January 2022
• Revised 25 April 2022
• Accepted 12 July 2023

DOI:https://doi.org/10.1103/PhysRevResearch.5.033078

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.

Published by the American Physical Society