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

Entropy

Using Matrix-Product States for Open Quantum Many-Body Systems: Efficient Algorithms for Markovian and Non-Markovian Time-Evolution

Regina Finsterhölzl, Manuel Katzer, Andreas Knorr, Alexander Carmele

Published Sep 4, 2020• 14 citations
Publisher

Abstract

This paper presents an efficient algorithm for the time evolution of open quantum many-body systems using matrix-product states (MPS) proposing a convenient structure of the MPS-architecture, which exploits the initial state of system and reservoir. By doing so, numerically expensive re-ordering protocols are circumvented. It is applicable to systems with a Markovian type of interaction, where only the present state of the reservoir needs to be taken into account. Its adaption to a non-Markovian type of interaction between the many-body system and the reservoir is demonstrated, where the information backflow from the reservoir needs to be included in the computation. Also, the derivation of the basis in the quantum stochastic Schrödinger picture is shown. As a paradigmatic model, the Heisenberg spin chain with nearest-neighbor interaction is used. It is demonstrated that the algorithm allows for the access of large systems sizes. As an example for a non-Markovian type of interaction, the generation of highly unusual steady states in the many-body system with coherent feedback control is demonstrated for a chain length of N=30.

DOI
10.3390/e22090984

Related papers

Non-Markovian open quantum many-body system dynamics

Regina Finsterhölzl

The study of open quantum systems forms one of the main problems of modern physics. This thesis aims at contributing to this field, and focuses on the investigation of non-Markovian dynamics of open quantum systems. Here, as a special case, quantum many-body systems provide for a rich field of research. With the ability to control them being of essential importance for the development of quantum technologies, the characterization of their dynamical properties has also drawn a lot of interest lately. One of the main methodological tasks in the study of open quantum systems is to reduce the infinite degrees of freedom of the total system. In order to achieve this, the thesis makes use of the tensor network method matrix product states and combines it with the picture of the quantum stochastic Schrödinger equation. However, when computing the dynamics of many-body systems of system size N, even the reduced system dynamics scale with 2N, making larger system sizes difficult to access. Therefore, especially when computing non-Markovian time-dynamics where the information backflow from the reservoir has to be taken into account, efficiently constructed algorithms are required. This thesis proposes these type of algorithms and demonstrates their application. They enable an access to larger system sizes, thus making it possible to investigate non-equilibrium steady-states of quantum many-body systems in a non-Markovian interaction for system sizes up to N=30. This thesis consists of three parts. The first part provides an introduction to the theory of open quantum systems and to the methodology for modeling them. The second part examines many-body system-reservoir interaction in the non-Markovian regime by investigating the case of fully coherent self-feedback. It uses the one-dimensional Heisenberg spin-1/2 chain as a paradigmatic and recently intensely studied many-body system. It is demonstrated that this gives rise to new conditions for the feedback phase and to highly non-trivial states within the system as well as to the possibility of partially characterizing the system non-invasively and storing excitation within it. Next, a novel approach to the boundary-driven Heisenberg chain is presented, where the chain is exposed to a coherent driving field, thus enabling a comparison to the well-studied incoherently driven model. The third part demonstrates the use of an external pump to control the feedback phase of a Lambda-type three-level system, where the external pump pulse is shown to be a new control parameter. The thesis concludes with the proposition to enhance the indistinguishability of two photons by feedback control.

•Jan 1, 2021•0 citations

Non-Markovian Stochastic Schr\"odinger Equation: Matrix Product State Approach

Xing Bao Gao, Jiajun Ren, Alexander Eisfeld +1 more

We derive a hierarchy of matrix product states (HOMPS) method which is numerically exact and efficient for general non-Markovian dynamics in open quantum system. This HOMPS is trying to attack the exponential wall issue in the recently developed hierarchy of pure states (HOPS) scheme with two steps: a. finding an effective time-dependent Schrodinger equation which is equivalent to HOPS, b. propagating this equation within matrix product states/operators (MPS/MPO) representation. HOMPS works in linear form and takes into account finite temperature effect straightforwardly from the initial pure state. Applications of HOMPS to spin-boson model covering both high and low temperatures are provided, demonstrating the validity and efficiency of the new approach.

•Sep 14, 2021•0 citations

Efficient time-evolution of matrix product states using average Hamiltonians

Belal Abouraya, Jirawat Saiphet, Fedor Jelezko +1 more

Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.

quant-pharXiv:2602.04955•Feb 6, 2026•0 citations

Non-Markovian stochastic Schrödinger equation: Matrix-product-state approach to the hierarchy of pure states

Xing Bao Gao, Jiajun Ren, Alexander Eisfeld +1 more

We derive a stochastic hierarchy of matrix product states (HOMPS) for non-Markovian dynamics in an open quantum system at finite temperature, which is numerically exact and efficient. HOMPS is obtained from the recently developed stochastic hierarchy of pure states (HOPS) by expressing HOPS in terms of formal creation and annihilation operators. The resulting stochastic first-order differential equation is then formulated in terms of matrix product states and matrix product operators. In this way the exponential complexity of HOPS can be reduced to scale polynomial with the number of particles. The validity and efficiency of HOMPS is demonstrated for the spin-boson model and long chains where each site is coupled to a structured, strongly non-Markovian environment.

arXiv:2109.06393•Physical Review A•Mar 31, 2022•35 citations

Improving the efficiency of open-quantum-system simulations using matrix product states in the interaction picture

Hanggai Nuomin, David N. Beratan, Peng Zhang

Modeling open quantum systems---quantum systems coupled to a bath---is of value in condensed-matter theory, cavity quantum electrodynamics, nanosciences, and biophysics. The real-time simulation of open quantum systems was advanced significantly by the recent development of chain mapping techniques and the use of matrix product states that exploit the intrinsic entanglement structure in open quantum systems. The computational cost of simulating open quantum systems, however, remains high when the bath is excited to high-lying quantum states. We develop an approach to reduce the computational costs in such cases. The interaction representation for the open quantum system is used to distribute excitations among the bath degrees of freedom so that the occupation of each bath oscillator is ensured to be low. The interaction picture also causes the matrix dimensions to be much smaller in a matrix product state of a chain-mapped open quantum system than in the Schr\"odinger picture. Using the interaction representation accelerates the calculations by one to two orders of magnitude over the existing matrix-product-state method. In the regime of strong system-bath coupling and high temperatures, the speedup can be as large as three orders of magnitude. The approach developed here is especially promising to simulate the dynamics of open quantum systems in high-temperature and strong-coupling regimes.

arXiv:2111.14308•Physical Review A•Mar 3, 2022•15 citations

Matrix Product State Simulations of Non-Equilibrium Steady States and Transient Heat Flows in the Two-Bath Spin-Boson Model at Finite Temperatures

Angus Dunnett, Alex W. Chin

Simulating the non-perturbative and non-Markovian dynamics of open quantum systems is a very challenging many body problem, due to the need to evolve both the system and its environments on an equal footing. Tensor network and matrix product states (MPS) have emerged as powerful tools for open system models, but the numerical resources required to treat finite-temperature environments grow extremely rapidly and limit their applications. In this study we use time-dependent variational evolution of MPS to explore the striking theory of Tamascelli et al. (Phys. Rev. Lett. 2019, 123, 090402.) that shows how finite-temperature open dynamics can be obtained from zero temperature, i.e., pure wave function, simulations. Using this approach, we produce a benchmark dataset for the dynamics of the Ohmic spin-boson model across a wide range of coupling strengths and temperatures, and also present a detailed analysis of the numerical costs of simulating non-equilibrium steady states, such as those emerging from the non-perturbative coupling of a qubit to baths at different temperatures. Despite ever-growing resource requirements, we find that converged non-perturbative results can be obtained, and we discuss a number of recent ideas and numerical techniques that should allow wide application of MPS to complex open quantum systems.

Entropy•Jan 6, 2021•16 citations
This paper presents an efficient algorithm for the time evolution of open quantum many-body systems using matrix-product states (MPS) proposing a convenient structure of the MPS-architecture, which exploits the initial state of system and reservoir. By doing so, numerically expensive re-ordering protocols are circumvented. It is applicable to systems with a Markovian type of interaction, where only the present state of the reservoir needs to be taken into account. Its adaption to a non-Markovian type of interaction between the many-body system and the reservoir is demonstrated, where the information backflow from the reservoir needs to be included in the computation. Also, the derivation of the basis in the quantum stochastic Schrödinger picture is shown. As a paradigmatic model, the Heisenberg spin chain with nearest-neighbor interaction is used. It is demonstrated that the algorithm allows for the access of large systems sizes. As an example for a non-Markovian type of interaction, the generation of highly unusual steady states in the many-body system with coherent feedback control is demonstrated for a chain length of N=30.
The study of open quantum systems forms one of the main problems of modern physics. This thesis aims at contributing to this field, and focuses on the investigation of non-Markovian dynamics of open quantum systems. Here, as a special case, quantum many-body systems provide for a rich field of research. With the ability to control them being of essential importance for the development of quantum technologies, the characterization of their dynamical properties has also drawn a lot of interest lately. One of the main methodological tasks in the study of open quantum systems is to reduce the infinite degrees of freedom of the total system. In order to achieve this, the thesis makes use of the tensor network method matrix product states and combines it with the picture of the quantum stochastic Schrödinger equation. However, when computing the dynamics of many-body systems of system size N, even the reduced system dynamics scale with 2N, making larger system sizes difficult to access. Therefore, especially when computing non-Markovian time-dynamics where the information backflow from the reservoir has to be taken into account, efficiently constructed algorithms are required. This thesis proposes these type of algorithms and demonstrates their application. They enable an access to larger system sizes, thus making it possible to investigate non-equilibrium steady-states of quantum many-body systems in a non-Markovian interaction for system sizes up to N=30. This thesis consists of three parts. The first part provides an introduction to the theory of open quantum systems and to the methodology for modeling them. The second part examines many-body system-reservoir interaction in the non-Markovian regime by investigating the case of fully coherent self-feedback. It uses the one-dimensional Heisenberg spin-1/2 chain as a paradigmatic and recently intensely studied many-body system. It is demonstrated that this gives rise to new conditions for the feedback phase and to highly non-trivial states within the system as well as to the possibility of partially characterizing the system non-invasively and storing excitation within it. Next, a novel approach to the boundary-driven Heisenberg chain is presented, where the chain is exposed to a coherent driving field, thus enabling a comparison to the well-studied incoherently driven model. The third part demonstrates the use of an external pump to control the feedback phase of a Lambda-type three-level system, where the external pump pulse is shown to be a new control parameter. The thesis concludes with the proposition to enhance the indistinguishability of two photons by feedback control.
We derive a hierarchy of matrix product states (HOMPS) method which is numerically exact and efficient for general non-Markovian dynamics in open quantum system. This HOMPS is trying to attack the exponential wall issue in the recently developed hierarchy of pure states (HOPS) scheme with two steps: a. finding an effective time-dependent Schrodinger equation which is equivalent to HOPS, b. propagating this equation within matrix product states/operators (MPS/MPO) representation. HOMPS works in linear form and takes into account finite temperature effect straightforwardly from the initial pure state. Applications of HOMPS to spin-boson model covering both high and low temperatures are provided, demonstrating the validity and efficiency of the new approach.
Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.
We derive a stochastic hierarchy of matrix product states (HOMPS) for non-Markovian dynamics in an open quantum system at finite temperature, which is numerically exact and efficient. HOMPS is obtained from the recently developed stochastic hierarchy of pure states (HOPS) by expressing HOPS in terms of formal creation and annihilation operators. The resulting stochastic first-order differential equation is then formulated in terms of matrix product states and matrix product operators. In this way the exponential complexity of HOPS can be reduced to scale polynomial with the number of particles. The validity and efficiency of HOMPS is demonstrated for the spin-boson model and long chains where each site is coupled to a structured, strongly non-Markovian environment.
Modeling open quantum systems---quantum systems coupled to a bath---is of value in condensed-matter theory, cavity quantum electrodynamics, nanosciences, and biophysics. The real-time simulation of open quantum systems was advanced significantly by the recent development of chain mapping techniques and the use of matrix product states that exploit the intrinsic entanglement structure in open quantum systems. The computational cost of simulating open quantum systems, however, remains high when the bath is excited to high-lying quantum states. We develop an approach to reduce the computational costs in such cases. The interaction representation for the open quantum system is used to distribute excitations among the bath degrees of freedom so that the occupation of each bath oscillator is ensured to be low. The interaction picture also causes the matrix dimensions to be much smaller in a matrix product state of a chain-mapped open quantum system than in the Schr\"odinger picture. Using the interaction representation accelerates the calculations by one to two orders of magnitude over the existing matrix-product-state method. In the regime of strong system-bath coupling and high temperatures, the speedup can be as large as three orders of magnitude. The approach developed here is especially promising to simulate the dynamics of open quantum systems in high-temperature and strong-coupling regimes.
Simulating the non-perturbative and non-Markovian dynamics of open quantum systems is a very challenging many body problem, due to the need to evolve both the system and its environments on an equal footing. Tensor network and matrix product states (MPS) have emerged as powerful tools for open system models, but the numerical resources required to treat finite-temperature environments grow extremely rapidly and limit their applications. In this study we use time-dependent variational evolution of MPS to explore the striking theory of Tamascelli et al. (Phys. Rev. Lett. 2019, 123, 090402.) that shows how finite-temperature open dynamics can be obtained from zero temperature, i.e., pure wave function, simulations. Using this approach, we produce a benchmark dataset for the dynamics of the Ohmic spin-boson model across a wide range of coupling strengths and temperatures, and also present a detailed analysis of the numerical costs of simulating non-equilibrium steady states, such as those emerging from the non-perturbative coupling of a qubit to baths at different temperatures. Despite ever-growing resource requirements, we find that converged non-perturbative results can be obtained, and we discuss a number of recent ideas and numerical techniques that should allow wide application of MPS to complex open quantum systems.