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:2610.05035

Process tensor simulations of non-Markovian open quantum systems with large uncoupled subspaces

Thomas K. Bracht, Moritz Cygorek

quant-phSubmitted Oct 4, 2026• Announced Oct 6, 2026• 0 citations
arXivHTMLPDF

Abstract

Process tensors in matrix product operator representation (PT-MPOs) are a powerful tool to calculate the temporal evolution and multi-time correlations of open quantum systems. While much progress has been made in making the calculation of process tensors for larger system dimensions possible, the numerical cost of the subsequent propagation of the system typically scales as $\mathcal{O}(d^6)$, which renders PT-MPO methods prohibitively expensive for system dimensions above several dozens. Here, we show how PT-MPO methods can be used efficiently when only some of the system states couple to the environment, and how combination with a Taylor expansion method allows for a reduction of the numerical complexity down to $\mathcal{O}(d^3)$. We show that these methods enable efficient simulations of systems with several hundreds of states. This is particularly relevant for solid-state cavity-QED systems with multiple modes, or when calculating frequency-filtered correlation functions by coupling sensor systems to the system of interest.

Comments
7 pages, 3 figures

Related papers

Uniform process tensor approach for the calculation of multi-time correlation functions of non-Markovian open systems

Matteo Garbellini, Konrad Mickiewicz, Valentin Link +2 more

The process tensor framework to open quantum systems provides the most general description of multi-time correlations in non-Markovian quantum dynamics. A compressed representation of a process tensor in terms of matrix product operators (MPO) can be used for numerically exact calculations of multi-time correlation functions in systems strongly coupled to a non-Markovian reservoir. We show here that the numerical scaling for computing multi-dimensional spectra can be significantly improved using a time-translation invariant MPO representation of the process tensor obtained from the uniform time-evolving matrix product operator (uniTEMPO) method. In particular, this approach provides a spectral representation of the non-Markovian dynamics that gives direct access to correlation functions in Fourier-space, avoiding explicit real-time evolution. We calculate linear and 2D electronic spectra for an example system and discuss the performance and numerical scaling of our simulations.

quant-pharXiv:2603.04970•The Journal of Chemical Physics•Mar 6, 2026•0 citations

Efficient construction of time-invariant process tensors for simulating high-dimensional non-Markovian open quantum systems

Émile Cochin, Jonathan Keeling, Brendon W. Lovett +1 more

Numerical methods for obtaining exact dynamics of non-Markovian open quantum systems are mostly limited to either small systems or to short-time evolution only. Here, we propose a new algorithm for computing process tensors--matrix product operator (MPO) representations that capture the environment influence--which achieves greatly enhanced computational scalings with system size, while maintaining linear scaling with simulation length. We build on recent developments in the field which allow for long-time evolutions through process tensors which have a time-translational invariance. These can be built for general Gaussian environments and generic coupling operators with the system using infinite time-evolving block decimation (iTEBD). We introduce a modified iTEBD algorithm using intermediate compression steps which bring down the computation time scaling with system size $d$ from $\mathcal{O}(d^8)$ to $\mathcal{O}(d^4)$, as well as significantly lowering the required memory. To illustrate the power of this method, we apply it to the problem of dispersive qubit readout in circuit QED, which was previously out-of-reach numerically. The full treatment of the measurement resonator, which requires a large system space, combined with the long simulation times precipitated by the separation of timescales between the measurement drive and the environment dissipation, is now possible. The algorithm we introduce not only allows for capturing non-Markovian dynamics in large open quantum systems, but also further extends all the existing capabilities of process tensors, for example in quantum optimal control, or in computation of multi-time correlations or of steady states, to more complex systems with tens of levels.

quant-pharXiv:2603.06840•Mar 10, 2026•0 citations

ACE: A general-purpose non-Markovian open quantum systems simulation toolkit based on process tensors

Moritz Cygorek, Erik M. Gauger

We describe a general-purpose computational toolkit for simulating open quantum systems, which provides numerically exact solutions for composites of zero-dimensional quantum systems that may be strongly coupled to multiple, quite general non-Markovian environments. It is based on process tensor matrix product operators (PT-MPOs), which efficiently encapsulate environment influences. The code features implementations of several PT-MPO algorithms, in particular Automated Compression of Environments for general environments comprised of independent modes as well as schemes for generalized spin boson models. The latter includes a divide-and-conquer scheme for periodic PT-MPOs, which enable million time step simulations for realistic models. PT-MPOs can be precalculated and reused for efficiently probing different time-dependent system Hamiltonians. They can also be stacked together and combined to provide numerically complete solutions of small networks of open quantum systems. The code is written in C++ and is fully controllable by configuration files, for which we have developed a versatile and compact human-readable format.

arXiv:2405.19319•The Journal of Chemical Physics•Aug 19, 2024•21 citations

Open Quantum System Dynamics from Infinite Tensor Network Contraction

Valentin Link, Hong-Hao Tu, Walter T. Strunz

Approaching the long-time dynamics of non-Markovian open quantum systems presents a challenging task if the bath is strongly coupled. Recent proposals address this problem through a representation of the so-called process tensor in terms of a tensor network. We show that for Gaussian environments highly efficient contraction to a matrix product operator (MPO) form can be achieved with infinite MPO evolution methods, leading to significant computational speed-up over existing proposals. The result structurally resembles open system evolution with carefully designed auxiliary degrees of freedom, as in hierarchical or pseudomode methods. Here, however, these degrees of freedom are generated automatically by the MPO evolution algorithm. Moreover, the semigroup form of the resulting propagator enables us to explore steady-state physics, such as phase transitions.

arXiv:2307.01802•Physical Review Letters•May 16, 2024•44 citations

Understanding and utilizing the inner bonds of process tensors

Moritz Cygorek, Erik M. Gauger

Process tensor matrix product operators (PT-MPOs) enable numerically exact simulations for an unprecedentedly broad range of open quantum systems. By representing environment influences in MPO form, they can be efficiently compressed using established algorithms. The dimensions of inner bonds of the compressed PT-MPO may be viewed as an indicator of the complexity of the environment. Here, we show that the inner bonds themselves, not only their dimensions, have a concrete physical meaning: They represent the subspace of the full environment Liouville space which hosts environment excitations that may influence the subsequent open quantum systems dynamics the most. This connection can be expressed in terms of lossy linear transformations, whose pseudoinverses facilitate the extraction of environment observables. We demonstrate this by extracting the environment spin of a central spin problem, the current through a quantum system coupled to two leads, the number of photons emitted from quantum emitters into a structured environment, and the distribution of the total absorbed energy in a driven non-Markovian quantum system into system, environment, and interaction energy terms. Numerical tests further indicate that different PT-MPO algorithms compress environments to similar subspaces. Thus, the physical interpretation of inner bonds of PT-MPOs both provides a conceptional understanding and it enables new practical applications.

arXiv:2404.01287•SciPost Physics•Jan 20, 2025•5 citations

Unifying methods for optimal control in non-Markovian quantum systems via process tensors

Carlos Ortega-Taberner, Eoin O'Neill, Eoin P. Butler +2 more

The large dimensionality of environments is the limiting factor in applying optimal control to open quantum systems beyond the Markovian approximation. Various methods exist to simulate non-Markovian systems, which effectively reduce the environment to a number of active degrees of freedom. Here, we show that several of these methods can be expressed in terms of a process tensor in the form of a matrix-product-operator, which serves as a unifying framework to show how they can be used in optimal control and to compare their performance. The matrix-product-operator form provides a general scheme for computing gradients using back propagation and allows the efficiency of the different methods to be compared via the bond dimensions of their respective process tensors.

arXiv:2406.17719•The Journal of Chemical Physics•Sep 28, 2024•7 citations
Process tensors in matrix product operator representation (PT-MPOs) are a powerful tool to calculate the temporal evolution and multi-time correlations of open quantum systems. While much progress has been made in making the calculation of process tensors for larger system dimensions possible, the numerical cost of the subsequent propagation of the system typically scales as O(d6)\mathcal{O}(d^6)O(d6), which renders PT-MPO methods prohibitively expensive for system dimensions above several dozens. Here, we show how PT-MPO methods can be used efficiently when only some of the system states couple to the environment, and how combination with a Taylor expansion method allows for a reduction of the numerical complexity down to O(d3)\mathcal{O}(d^3)O(d3). We show that these methods enable efficient simulations of systems with several hundreds of states. This is particularly relevant for solid-state cavity-QED systems with multiple modes, or when calculating frequency-filtered correlation functions by coupling sensor systems to the system of interest.
The process tensor framework to open quantum systems provides the most general description of multi-time correlations in non-Markovian quantum dynamics. A compressed representation of a process tensor in terms of matrix product operators (MPO) can be used for numerically exact calculations of multi-time correlation functions in systems strongly coupled to a non-Markovian reservoir. We show here that the numerical scaling for computing multi-dimensional spectra can be significantly improved using a time-translation invariant MPO representation of the process tensor obtained from the uniform time-evolving matrix product operator (uniTEMPO) method. In particular, this approach provides a spectral representation of the non-Markovian dynamics that gives direct access to correlation functions in Fourier-space, avoiding explicit real-time evolution. We calculate linear and 2D electronic spectra for an example system and discuss the performance and numerical scaling of our simulations.
Numerical methods for obtaining exact dynamics of non-Markovian open quantum systems are mostly limited to either small systems or to short-time evolution only. Here, we propose a new algorithm for computing process tensors--matrix product operator (MPO) representations that capture the environment influence--which achieves greatly enhanced computational scalings with system size, while maintaining linear scaling with simulation length. We build on recent developments in the field which allow for long-time evolutions through process tensors which have a time-translational invariance. These can be built for general Gaussian environments and generic coupling operators with the system using infinite time-evolving block decimation (iTEBD). We introduce a modified iTEBD algorithm using intermediate compression steps which bring down the computation time scaling with system size ddd from O(d8)\mathcal{O}(d^8)O(d8) to O(d4)\mathcal{O}(d^4)O(d4), as well as significantly lowering the required memory. To illustrate the power of this method, we apply it to the problem of dispersive qubit readout in circuit QED, which was previously out-of-reach numerically. The full treatment of the measurement resonator, which requires a large system space, combined with the long simulation times precipitated by the separation of timescales between the measurement drive and the environment dissipation, is now possible. The algorithm we introduce not only allows for capturing non-Markovian dynamics in large open quantum systems, but also further extends all the existing capabilities of process tensors, for example in quantum optimal control, or in computation of multi-time correlations or of steady states, to more complex systems with tens of levels.
We describe a general-purpose computational toolkit for simulating open quantum systems, which provides numerically exact solutions for composites of zero-dimensional quantum systems that may be strongly coupled to multiple, quite general non-Markovian environments. It is based on process tensor matrix product operators (PT-MPOs), which efficiently encapsulate environment influences. The code features implementations of several PT-MPO algorithms, in particular Automated Compression of Environments for general environments comprised of independent modes as well as schemes for generalized spin boson models. The latter includes a divide-and-conquer scheme for periodic PT-MPOs, which enable million time step simulations for realistic models. PT-MPOs can be precalculated and reused for efficiently probing different time-dependent system Hamiltonians. They can also be stacked together and combined to provide numerically complete solutions of small networks of open quantum systems. The code is written in C++ and is fully controllable by configuration files, for which we have developed a versatile and compact human-readable format.
Approaching the long-time dynamics of non-Markovian open quantum systems presents a challenging task if the bath is strongly coupled. Recent proposals address this problem through a representation of the so-called process tensor in terms of a tensor network. We show that for Gaussian environments highly efficient contraction to a matrix product operator (MPO) form can be achieved with infinite MPO evolution methods, leading to significant computational speed-up over existing proposals. The result structurally resembles open system evolution with carefully designed auxiliary degrees of freedom, as in hierarchical or pseudomode methods. Here, however, these degrees of freedom are generated automatically by the MPO evolution algorithm. Moreover, the semigroup form of the resulting propagator enables us to explore steady-state physics, such as phase transitions.
Process tensor matrix product operators (PT-MPOs) enable numerically exact simulations for an unprecedentedly broad range of open quantum systems. By representing environment influences in MPO form, they can be efficiently compressed using established algorithms. The dimensions of inner bonds of the compressed PT-MPO may be viewed as an indicator of the complexity of the environment. Here, we show that the inner bonds themselves, not only their dimensions, have a concrete physical meaning: They represent the subspace of the full environment Liouville space which hosts environment excitations that may influence the subsequent open quantum systems dynamics the most. This connection can be expressed in terms of lossy linear transformations, whose pseudoinverses facilitate the extraction of environment observables. We demonstrate this by extracting the environment spin of a central spin problem, the current through a quantum system coupled to two leads, the number of photons emitted from quantum emitters into a structured environment, and the distribution of the total absorbed energy in a driven non-Markovian quantum system into system, environment, and interaction energy terms. Numerical tests further indicate that different PT-MPO algorithms compress environments to similar subspaces. Thus, the physical interpretation of inner bonds of PT-MPOs both provides a conceptional understanding and it enables new practical applications.
The large dimensionality of environments is the limiting factor in applying optimal control to open quantum systems beyond the Markovian approximation. Various methods exist to simulate non-Markovian systems, which effectively reduce the environment to a number of active degrees of freedom. Here, we show that several of these methods can be expressed in terms of a process tensor in the form of a matrix-product-operator, which serves as a unifying framework to show how they can be used in optimal control and to compare their performance. The matrix-product-operator form provides a general scheme for computing gradients using back propagation and allows the efficiency of the different methods to be compared via the bond dimensions of their respective process tensors.