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:cond-mat/0405525 · Physical Review B

Rigorous Born approximation and beyond for the spin-boson model

David Peter Divincenzo, Daniel Loss

Published Jan 14, 2005• 117 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.71.035318

Related papers

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

Dissipative quantum mechanics beyond the Bloch-Redfield approximation: A consistent weak-coupling expansion of the Ohmic spin boson model at arbitrary bias

Carsten J. Lindner, Herbert Schoeller

We study the time evolution of the reduced density matrix for the Ohmic spin boson model out of an uncorrelated but otherwise arbitrary initial state. We consider arbitrary bias $\ensuremath{\epsilon}$ and tunneling $\mathrm{\ensuremath{\Delta}}$ at zero temperature for a weak coupling $\ensuremath{\alpha}$ to the bosonic bath. Using the real-time renormalization group method, we present a consistent weak-coupling expansion one order beyond the Bloch-Redfield approximation within a renormalized perturbation theory with analytical results covering the whole crossover regime from small times $\mathrm{\ensuremath{\Omega}}t\ensuremath{\ll}1$ to large times $\mathrm{\ensuremath{\Omega}}t\ensuremath{\gg}1$, where $\mathrm{\ensuremath{\Omega}}=\sqrt{{\ensuremath{\epsilon}}^{2}+{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\mathrm{\ensuremath{\Delta}}}}^{2}}$ denotes the Rabi frequency in terms of the renormalized tunneling $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\mathrm{\ensuremath{\Delta}}}$. In addition, for exponentially small or large times, we perform a nonperturbative resummation of all logarithmic terms. We show that standard Born approximation schemes calculating the effective Liouvillian of the kinetic equation up to first order in $\ensuremath{\alpha}$ are not sufficient to account for various important corrections one order beyond the Bloch-Redfield solution. (1) The resummation of all secular terms $\ensuremath{\sim}{(\mathrm{\ensuremath{\Gamma}}t)}^{n}$ is necessary to obtain the correct exponential decay of all terms of the time evolution with decay rate $\mathrm{\ensuremath{\Gamma}}$ or $\mathrm{\ensuremath{\Gamma}}/2$, together with the correct pre-exponential functions. (2) The resummation of all logarithmic terms at high and low energies leads to a renormalized tunneling $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\mathrm{\ensuremath{\Delta}}}$ and to pre-exponential functions of logarithmic and power-law form. (3) The fact that two eigenvalues of $L(E)$ are close to each other by $O(\mathrm{\ensuremath{\Gamma}})$ requires degenerate perturbation theory for times $\mathrm{\ensuremath{\Gamma}}t\ensuremath{\sim}O(1)$, where certain terms of the Liouvillian in $O({\ensuremath{\alpha}}^{2})$ are needed to calculate the stationary state and the time evolution of the nonoscillating purely decaying modes up to $O(\ensuremath{\alpha})$. In contrast to the zero-bias case, we find two further interesting results for the time dynamics of the oscillating modes. (4) The terms of the pre-exponential functions with a strong time dependence show a leading long-time tail $\ensuremath{\sim}\ensuremath{\alpha}/(\mathrm{\ensuremath{\Omega}}t)$, besides other subleading terms $\ensuremath{\sim}\ensuremath{\alpha}/{(\mathrm{\ensuremath{\Omega}}t)}^{2}$ well-known from the zero-bias case. (5) The terms of the pre-exponential functions with a weak (logarithmic) time dependence vary according to a power law ${(\frac{1}{\mathrm{\ensuremath{\Omega}}t})}^{2\ensuremath{\alpha}\frac{{\ensuremath{\epsilon}}^{2}}{{\mathrm{\ensuremath{\Omega}}}^{2}}}$ for exponentially large times. The power-law exponent depends on the bias and has to be contrasted to the one at exponentially small times where it crosses over to the bias-independent result $2\ensuremath{\alpha}$. We discuss that the complexity to calculate one order beyond Bloch-Redfield approximation is rather generic and applies also to other models of dissipative quantum mechanics.

arXiv:1802.09846•Physical review. B./Physical review. B•Sep 17, 2018•16 citations

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

Ilya Vilkoviskiy, Dmitry A. Abanin

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.

arXiv:2307.15592•Jul 28, 2023•0 citations

Non-Markovian effects in the spin-boson model at zero temperature

Sebastian Wenderoth, Heinz‐Peter Breuer, Michael Thoss

We investigate memory effects in the spin-boson model using a recently proposed measure for non-Markovian behavior based on the information exchange between an open system and its environment. Employing the numerical exact multilayer multiconfiguration time-dependent Hartree approach, we simulate the dynamics of the spin-boson model at zero temperature for a broad range of parameters. For a fast bath, i.e., in the scaling limit, we find non-Markovian dynamics for a coherently decaying spin at weak system-bath coupling, whereas memory effects are absent for stronger coupling in the regimes of incoherent decay and localization. If the timescales of system and bath are comparable, a complex, nonmonotonic dependence of non-Markovianity on the system-bath coupling strength is observed.

arXiv:2101.09463•Physical Review A•Jul 19, 2021•26 citations

Dynamical renormalization-group approach to the spin-boson model

Hassan Shapourian

We develop a semianalytical approach beyond the Born-Markov approximation to study the quench dynamics of the spin-boson model in the strong-coupling regime ($\ensuremath{\alpha}\ensuremath{\le}1/2$) for the Ohmic bath. The basic idea in our approach is to write an effective time-dependent model for the dynamics of the system coupled to the bosonic bath after integrating out high-frequency bath modes. By applying this procedure to the Heisenberg equations of motion, we derive a set of flow equations for the system parameters as a function of time. The final flow equations look similar to those of the equilibrium renormalization-group theory; however, in our derivation the scaling parameter is set by the real time. We solve the equations of motion with time-dependent renormalized parameters and show that the resulting dynamics is in decent agreement with the exact NRG calculations as well as the noninteracting blip approximation that is a well-known good solution in this limit.

arXiv:1607.03920•Physical Review A•Mar 14, 2016•11 citations

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
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 \ensuremathα\sqrt{\ensuremath{\alpha}}\ensuremathα​, where \ensuremathα\ensuremath{\alpha}\ensuremathα 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α)O(\ensuremath{\alpha})O(\ensuremathα) 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α\ensuremath{\alpha}\ensuremathα. All these calculations indicate precision experimental tests that could confirm or refute the validity of the spin-boson model in a variety of systems.
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.
We study the time evolution of the reduced density matrix for the Ohmic spin boson model out of an uncorrelated but otherwise arbitrary initial state. We consider arbitrary bias $\ensuremath{\epsilon}$ and tunneling $\mathrm{\ensuremath{\Delta}}$ at zero temperature for a weak coupling $\ensuremath{\alpha}$ to the bosonic bath. Using the real-time renormalization group method, we present a consistent weak-coupling expansion one order beyond the Bloch-Redfield approximation within a renormalized perturbation theory with analytical results covering the whole crossover regime from small times $\mathrm{\ensuremath{\Omega}}t\ensuremath{\ll}1$ to large times $\mathrm{\ensuremath{\Omega}}t\ensuremath{\gg}1$, where $\mathrm{\ensuremath{\Omega}}=\sqrt{{\ensuremath{\epsilon}}^{2}+{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\mathrm{\ensuremath{\Delta}}}}^{2}}$ denotes the Rabi frequency in terms of the renormalized tunneling $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\mathrm{\ensuremath{\Delta}}}$. In addition, for exponentially small or large times, we perform a nonperturbative resummation of all logarithmic terms. We show that standard Born approximation schemes calculating the effective Liouvillian of the kinetic equation up to first order in $\ensuremath{\alpha}$ are not sufficient to account for various important corrections one order beyond the Bloch-Redfield solution. (1) The resummation of all secular terms $\ensuremath{\sim}{(\mathrm{\ensuremath{\Gamma}}t)}^{n}$ is necessary to obtain the correct exponential decay of all terms of the time evolution with decay rate $\mathrm{\ensuremath{\Gamma}}$ or $\mathrm{\ensuremath{\Gamma}}/2$, together with the correct pre-exponential functions. (2) The resummation of all logarithmic terms at high and low energies leads to a renormalized tunneling $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\mathrm{\ensuremath{\Delta}}}$ and to pre-exponential functions of logarithmic and power-law form. (3) The fact that two eigenvalues of L(E)L(E)L(E) are close to each other by $O(\mathrm{\ensuremath{\Gamma}})$ requires degenerate perturbation theory for times $\mathrm{\ensuremath{\Gamma}}t\ensuremath{\sim}O(1)$, where certain terms of the Liouvillian in $O({\ensuremath{\alpha}}^{2})$ are needed to calculate the stationary state and the time evolution of the nonoscillating purely decaying modes up to $O(\ensuremath{\alpha})$. In contrast to the zero-bias case, we find two further interesting results for the time dynamics of the oscillating modes. (4) The terms of the pre-exponential functions with a strong time dependence show a leading long-time tail $\ensuremath{\sim}\ensuremath{\alpha}/(\mathrm{\ensuremath{\Omega}}t)$, besides other subleading terms $\ensuremath{\sim}\ensuremath{\alpha}/{(\mathrm{\ensuremath{\Omega}}t)}^{2}$ well-known from the zero-bias case. (5) The terms of the pre-exponential functions with a weak (logarithmic) time dependence vary according to a power law ${(\frac{1}{\mathrm{\ensuremath{\Omega}}t})}^{2\ensuremath{\alpha}\frac{{\ensuremath{\epsilon}}^{2}}{{\mathrm{\ensuremath{\Omega}}}^{2}}}$ for exponentially large times. The power-law exponent depends on the bias and has to be contrasted to the one at exponentially small times where it crosses over to the bias-independent result $2\ensuremath{\alpha}$. We discuss that the complexity to calculate one order beyond Bloch-Redfield approximation is rather generic and applies also to other models of dissipative quantum mechanics.
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.
We investigate memory effects in the spin-boson model using a recently proposed measure for non-Markovian behavior based on the information exchange between an open system and its environment. Employing the numerical exact multilayer multiconfiguration time-dependent Hartree approach, we simulate the dynamics of the spin-boson model at zero temperature for a broad range of parameters. For a fast bath, i.e., in the scaling limit, we find non-Markovian dynamics for a coherently decaying spin at weak system-bath coupling, whereas memory effects are absent for stronger coupling in the regimes of incoherent decay and localization. If the timescales of system and bath are comparable, a complex, nonmonotonic dependence of non-Markovianity on the system-bath coupling strength is observed.
We develop a semianalytical approach beyond the Born-Markov approximation to study the quench dynamics of the spin-boson model in the strong-coupling regime ($\ensuremath{\alpha}\ensuremath{\le}1/2$) for the Ohmic bath. The basic idea in our approach is to write an effective time-dependent model for the dynamics of the system coupled to the bosonic bath after integrating out high-frequency bath modes. By applying this procedure to the Heisenberg equations of motion, we derive a set of flow equations for the system parameters as a function of time. The final flow equations look similar to those of the equilibrium renormalization-group theory; however, in our derivation the scaling parameter is set by the real time. We solve the equations of motion with time-dependent renormalized parameters and show that the resulting dynamics is in decent agreement with the exact NRG calculations as well as the noninteracting blip approximation that is a well-known good solution in this limit.
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.