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

A bound on approximating non-Markovian dynamics by tensor networks in the time domain

Ilya Vilkoviskiy, Dmitry A. Abanin

Published Jul 28, 2023• 0 citations
arXivHTMLPDF

Abstract

Spin-boson (SB) model plays a central role in studies of dissipative quantum dynamics, both due its conceptual importance and relevance to a number of physical systems. Here we provide rigorous bounds of the computational complexity of the SB model for the physically relevant case of a zero temperature Ohmic bath. We start with the description of the bosonic bath via its Feynman-Vernon influence functional (IF), which is a tensor on the space of spin's trajectories. By expanding the kernel of the IF functional via a sum of decaying exponentials, we obtain an analytical approximation of the continuous bath by a finite number of damped bosonic modes. We bound the error induced by restricting bosonic Hilbert spaces to a finite-dimensional subspace with small boson numbers, which yields an analytical form of a matrix-product state (MPS) representation of the IF. We show that the MPS bond dimension $D$ scales polynomially in the error on physical observables $ε$, as well as in the evolution time $T$, $D\propto T^4/ε^2$. This bound indicates that the spin-boson model can be efficiently simulated using polynomial in time computational resources.

DOI
10.48550/arxiv.2307.15592

Related papers

Bound on approximating non-Markovian dynamics by tensor networks in the time domain

Ilya Vilkoviskiy, Dmitry A. Abanin

The spin-boson (SB) model plays a central role in studies of dissipative quantum dynamics, due to bothits conceptual importance and relevance to a number of physical systems. Here, we provide rigorous bounds of the computational complexity of the SB model for the physically relevant case of a zero temperature ohmic bath. We start with the description of the bosonic bath via its Feynman-Vernon influence functional (IF), which is a tensor on the space of the trajectory of an impurity spin. By expanding the kernel of the IF via a sum of decaying exponentials, we obtain an analytical approximation of the continuous bath by a finite number of damped bosonic modes. We bound the error induced by restricting bosonic Hilbert spaces to a finite-dimensional subspace with small boson numbers, which yields an analytical form of a matrix-product state (MPS) representation of the IF. We show that the MPS bond dimension D scales polynomially in the error on physical observables as well as in the evolution time T, $D\ensuremath{\propto}{T}^{4}/{\ensuremath{\epsilon}}^{2}$. This bound indicates that the SB model can be efficiently simulated using a polynomial in time-computational resources.

Physical review. B./Physical review. B•May 9, 2024•16 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

Simulating Bosonic Baths with Error Bars

Mischa P. Woods, M. Cramer, Martin Bodo Plenio

We derive rigorous truncation-error bounds for the spin-boson model and its generalizations to arbitrary quantum systems interacting with bosonic baths. For the numerical simulation of such baths, the truncation of both the number of modes and the local Hilbert-space dimensions is necessary. We derive superexponential Lieb-Robinson-type bounds on the error when restricting the bath to finitely many modes and show how the error introduced by truncating the local Hilbert spaces may be efficiently monitored numerically. In this way we give error bounds for approximating the infinite system by a finite-dimensional one. As a consequence, numerical simulations such as the time-evolving density with orthogonal polynomials algorithm (TEDOPA) now allow for the fully certified treatment of the system-environment interaction.

arXiv:1504.01531•Physical Review Letters•Sep 22, 2015•58 citations

Rigorous Born approximation and beyond for the spin-boson model

David Peter Divincenzo, Daniel Loss

Within the lowest-order Born approximation, we present an exact calculation of the time dynamics of the spin-boson model in the Ohmic regime. We observe non-Markovian effects at zero temperature that scale with the system-bath coupling strength and cause qualitative changes in the evolution of coherence at intermediate times of order of the oscillation period. These changes could significantly affect the performance of these systems as qubits. In the biased case, we find a prompt loss of coherence at these intermediate times, whose decay rate is set by $\sqrt{\ensuremath{\alpha}}$, where $\ensuremath{\alpha}$ is the coupling strength to the environment. We also explore the calculation of the next-order Born approximation: we show that, at the expense of very large computational complexity, interesting physical quantities can be rigorously computed at fourth order using computer algebra, presented completely in an accompanying MATHEMATICA file. We compute the $O(\ensuremath{\alpha})$ corrections to the long time behavior of the system density matrix; the result is identical to the reduced density matrix of the equilibrium state to the same order in $\ensuremath{\alpha}$. All these calculations indicate precision experimental tests that could confirm or refute the validity of the spin-boson model in a variety of systems.

arXiv:cond-mat/0405525•Physical Review B•Jan 14, 2005•117 citations

Tensor network simulations of open quantum systems

Dainius Kilda

The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.

•Mar 3, 2020•0 citations

Exact Born Approximation for the Spin-Boson Model

Daniel Loss, David Peter Divincenzo

Within the lowest-order Born approximation, we present an exact calculation of the time dynamics of the spin-boson model in the Ohmic regime. We observe non-Markovian effects at zero temperature that scale with the system-bath coupling strength and cause qualitative changes in the evolution of coherence at intermediate times of order of the oscillation period. These changes could significantly affect the performance of these systems as qubits. In the biased case, we find a prompt loss of coherence at these intermediate times, whose decay rate is set by $\sqrtα$, where $α$ is the coupling strength to the environment. These calculations indicate precision experimental tests that could confirm or refute the validity of the spin-boson model in a variety of systems.

arXiv:cond-mat/0304118•Apr 4, 2003•9 citations
Spin-boson (SB) model plays a central role in studies of dissipative quantum dynamics, both due its conceptual importance and relevance to a number of physical systems. Here we provide rigorous bounds of the computational complexity of the SB model for the physically relevant case of a zero temperature Ohmic bath. We start with the description of the bosonic bath via its Feynman-Vernon influence functional (IF), which is a tensor on the space of spin's trajectories. By expanding the kernel of the IF functional via a sum of decaying exponentials, we obtain an analytical approximation of the continuous bath by a finite number of damped bosonic modes. We bound the error induced by restricting bosonic Hilbert spaces to a finite-dimensional subspace with small boson numbers, which yields an analytical form of a matrix-product state (MPS) representation of the IF. We show that the MPS bond dimension DDD scales polynomially in the error on physical observables εεε, as well as in the evolution time TTT, D∝T4/ε2D\propto T^4/ε^2D∝T4/ε2. This bound indicates that the spin-boson model can be efficiently simulated using polynomial in time computational resources.
The spin-boson (SB) model plays a central role in studies of dissipative quantum dynamics, due to bothits conceptual importance and relevance to a number of physical systems. Here, we provide rigorous bounds of the computational complexity of the SB model for the physically relevant case of a zero temperature ohmic bath. We start with the description of the bosonic bath via its Feynman-Vernon influence functional (IF), which is a tensor on the space of the trajectory of an impurity spin. By expanding the kernel of the IF via a sum of decaying exponentials, we obtain an analytical approximation of the continuous bath by a finite number of damped bosonic modes. We bound the error induced by restricting bosonic Hilbert spaces to a finite-dimensional subspace with small boson numbers, which yields an analytical form of a matrix-product state (MPS) representation of the IF. We show that the MPS bond dimension D scales polynomially in the error on physical observables as well as in the evolution time T, $D\ensuremath{\propto}{T}^{4}/{\ensuremath{\epsilon}}^{2}$. This bound indicates that the SB model can be efficiently simulated using a polynomial in time-computational resources.
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.
We derive rigorous truncation-error bounds for the spin-boson model and its generalizations to arbitrary quantum systems interacting with bosonic baths. For the numerical simulation of such baths, the truncation of both the number of modes and the local Hilbert-space dimensions is necessary. We derive superexponential Lieb-Robinson-type bounds on the error when restricting the bath to finitely many modes and show how the error introduced by truncating the local Hilbert spaces may be efficiently monitored numerically. In this way we give error bounds for approximating the infinite system by a finite-dimensional one. As a consequence, numerical simulations such as the time-evolving density with orthogonal polynomials algorithm (TEDOPA) now allow for the fully certified treatment of the system-environment interaction.
Within the lowest-order Born approximation, we present an exact calculation of the time dynamics of the spin-boson model in the Ohmic regime. We observe non-Markovian effects at zero temperature that scale with the system-bath coupling strength and cause qualitative changes in the evolution of coherence at intermediate times of order of the oscillation period. These changes could significantly affect the performance of these systems as qubits. In the biased case, we find a prompt loss of coherence at these intermediate times, whose decay rate is set by $\sqrt{\ensuremath{\alpha}}$, where $\ensuremath{\alpha}$ is the coupling strength to the environment. We also explore the calculation of the next-order Born approximation: we show that, at the expense of very large computational complexity, interesting physical quantities can be rigorously computed at fourth order using computer algebra, presented completely in an accompanying MATHEMATICA file. We compute the $O(\ensuremath{\alpha})$ corrections to the long time behavior of the system density matrix; the result is identical to the reduced density matrix of the equilibrium state to the same order in $\ensuremath{\alpha}$. All these calculations indicate precision experimental tests that could confirm or refute the validity of the spin-boson model in a variety of systems.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
Within the lowest-order Born approximation, we present an exact calculation of the time dynamics of the spin-boson model in the Ohmic regime. We observe non-Markovian effects at zero temperature that scale with the system-bath coupling strength and cause qualitative changes in the evolution of coherence at intermediate times of order of the oscillation period. These changes could significantly affect the performance of these systems as qubits. In the biased case, we find a prompt loss of coherence at these intermediate times, whose decay rate is set by α\sqrtαα​, where ααα is the coupling strength to the environment. These calculations indicate precision experimental tests that could confirm or refute the validity of the spin-boson model in a variety of systems.