Abstract:While several numerical techniques are available for predicting the dynamics of non-Markovian open quantum systems, most struggle with simulations for very long memory and propagation times, e.g., due to superlinear scaling with the number of time steps n . Here, we introduce a numerically exact algorithm to calculate process tensors—compact representations of environmental influences—which provides a scaling advantage over previous algorithms by leveraging self-similarity of the tensor networks that represent the environment. It is applicable to environments with Gaussian statistics, such as for spin-boson-type open quantum systems. Based on a divide-and-conquer strategy, our approach requires only O(n log n) singular value decompositions for environments with infinite memory. Where the memory can be truncated after nc time steps, a nominal scaling O(nc log nc) is found, which is independent of n . This improved scaling is enabled by identifying process tensors with repeatable blocks. To demonstrate the power and utility of our approach, we provide three examples. (1) We calculate the fluorescence spectra of a quantum dot under both strong driving and strong dot-phonon couplings, a task requiring simulations over millions of time steps, which we are able to perform in minutes. (2) We efficiently find process tensors describing superradiance of multiple emitters. (3) We explore the limits of our algorithm by considering coherence decay with a very strongly coupled environment. The observed computation time is not necessarily proportional to the number of singular value decompositions because the matrix dimensions also depend on the number of time steps. Nevertheless, quasilinear and sublinear scaling of computation time is found in practice for a wide range of parameters. While there are instances where existing methods can achieve comparable nominal scaling by precalculating effective propagators for time-independent or periodic system Hamiltonians, process tensors contain all the information needed to extract arbitrary multitime correlation functions of the system when driven by arbitrary time-dependent system Hamiltonians. The algorithm we present here not only significantly extends the scope of numerically exact techniques to open quantum systems with long memory times, but it also has fundamental implications for the simulation complexity of tensor network approaches. https://doi.org/10.1103/PhysRevX.14.011010 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

Taming NSPT fluctuations in $O(N)$ Non-Linear Sigma Model: simulations in the large $N$ regime

NSPT for $O(N)$ non-linear sigma model: the larger $N$ the better

Large fluctuations in NSPT computations: a lesson from $O(N)$ non-linear sigma models

Large $N$ Limit of the $O(N)$ Linear Sigma Model Via Stochastic Quantization

Large N Limit of the O ( N ) Linear Sigma Model in 3D

Large $N$ limit and $1/N$ expansion of invariant observables in $O(N)$ linear $\sigma$-model via SPDE

Large $N$ limit and $1/N$ expansion of invariant observables in $O(N)$ linear $σ$-model via SPDE

Conformal invariance constraints in the $O(N)$ models: a first study within the nonperturbative renormalization group

Studies of critical phenomena and phase transitions in large lattices with Monte-Carlo based non-perturbative approaches

Numerical Simulation of Large-Scale Nonlinear Open Quantum Mechanics

$O(N)\times O(2)$ scalar models: including $\mathcal{O}(\partial^2)$ corrections in the Functional Renormalization Group analysis

Nonlinear sigma models for monitored dynamics of free fermions

Non-local non-linear sigma models

Fluctuations in the non-linear supersymmetric hyperbolic sigma model with long-range interactions

Sublinear Scaling in Non-Markovian Open Quantum Systems Simulations

Nonlinear Geometric Optics Based Multiscale Stochastic Galerkin Methods for Highly Oscillatory Transport Equations with Random Inputs

Critical Exponents of the Nonlinear Sigma Model on a Grassmann Manifold U(N)/U(m)U(N - M) by the 1/N-Expansion

On Real Time Dynamics of Large $N$ Models

Nonlinear dynamics for the Ising model

Modeling Nonlinear Optics of Nanosystems with Sum-Over-States Model