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/0411443

Path integral derivation of Bloch-Redfield equations for a qubit weakly coupled to a heat bath: Application to nonadiabatic transitions

N. W. Watkins, David Waxman

Published Nov 17, 2004• 0 citations
arXivHTMLPDF

Abstract

Quantum information processing has greatly increased interest in the phenomenon of environmentally-induced decoherence. The spin boson model is widely used to study the interaction between a spin-modelling a quantum particle moving in a double well potential-and its environment-modelled by a heat bath of harmonic oscillators. This paper extends a previous analysis of the static spin boson study to the driven spin boson case, with the derivation of an exact integro-differential equation for the time evolution of the propagator of the reduced spin density matrix. This is the first main result. By specializing to weak damping we then obtain the next result, a set of Bloch-Redfield equations for the equilibrium fixed spin initial condition. Finally we show that these equations can be used to solve the classic dissipative Landau-Zener problem and illustrate these solutions for the weak damping case. The effect of dissipation is seen to be minimised as the speed of passage is increased, implying that qubits need to be switched as fast as possible.

DOI
10.48550/arxiv.cond-mat/0411443

Related papers

Nonadiabatic dynamics of a slowly driven dissipative two-level system

Canran Xu, Amrit Poudel, Maxim Vavilov

We study the dynamics of a two-level system described by a slowly varying Hamiltonian and weakly coupled to the Ohmic environment. We follow the Bloch-Redfield perturbative approach to include the effect of the environment on qubit evolution and take into account modification of the spectrum and matrix elements of qubit transitions due to the time dependence of the Hamiltonian. We apply this formalism to two problems. (1) We consider a qubit, or a spin-$\frac{1}{2}$, in a rotating magnetic field. We show that once the rotation starts, the spin has a component perpendicular to the rotation plane of the field that initially wiggles and eventually settles to the value proportional to the product of angular rotation velocity of the field and the Berry curvature. (2) We reexamine the Landau-Zener transition for a system coupled to the environment at arbitrary temperature. We show that as temperature increases, the thermal excitation and relaxation become leading processes responsible for transition between states of the system. We also apply the Lindblad master equations to these two problems and compare results with those obtained from the Bloch-Redfield equations.

arXiv:1402.4210•Physical Review A•May 6, 2014•41 citations

Qubit-oscillator relationships in the open quantum Rabi model: the role of dissipation

G. Di Bello, Loris Maria Cangemi, V. Cataudella +3 more

Abstract Using a dissipative quantum Rabi model, we study the dynamics of a slow qubit coupled to a fast quantum harmonic oscillator interacting with a bosonic bath from weak to strong and ultra-strong coupling regimes. Solving the quantum Heisenberg equations of motion, perturbative in the internal coupling between qubit and oscillator, we derive functional relationships directly linking the qubit coordinates in the Bloch sphere to oscillator observables. We then perform accurate time-dependent Matrix Product State simulations and compare our results both with the analytical solutions of the Heisenberg equations of motion, and with numerical solutions of a Lindblad master equation, perturbative in the external coupling between oscillator and environment. Indeed, we show that, up to the strong coupling regime, the qubit state accurately fulfils the derived functional relationships. We analyse in detail the case of a qubit starting with generic coordinates on the Bloch sphere of which we evaluate the three components of the Bloch vector through the averages of oscillator observables. Interestingly, a weak to intermediate oscillator coupling to the bath is able to simplify the Bloch vector evaluation since qubit-oscillator relationships are more immediate. Moreover, by monitoring the qubit fidelity with respect to free limit, we find the parameter regime where the combined effect of internal and external couplings is able to hinder the reliable evaluation of the qubit Bloch vector. Finally, in the ultra-strong coupling regime, non-Markovian effects become robust and the dynamics of qubit and oscillator are inextricably entangled making the qubit Bloch vector evaluation difficult. Graphical Abstract

arXiv:2208.00708•The European Physical Journal Plus•Feb 8, 2023•2 citations

Dissipative dynamics of a driven qubit: Interplay between nonadiabatic dynamics and noise effects from the weak to strong coupling regime

Loris Maria Cangemi, V. Cataudella, Maura Sassetti +1 more

We study the exact solution of the Schr\"odinger equation for the dissipative dynamics of a qubit, achieved by means of short iterative Lanczos method (SIL), which allows us to describe the qubit and the bath dynamics from weak to strong coupling regimes. We focus on two different models of a qubit in contact with the external environment: the first is the spin boson model (SBM), which gives a description of the qubit in terms of static tunneling energy and a bias field. The second model describes an externally driven qubit, where both the bias field and the tunneling rate are controlled by a time-dependent magnetic field, following a finite time protocol. We show that in the SBM case, our solution correctly describes the crossover from coherent to incoherent behavior of the magnetization, occurring at the Toulouse point. Furthermore, we show that the bath response dramatically changes during the system dynamics, going from nonresonant at small times to resonant behavior at long times. When the external driving field is present, for fixed values of the drive duration our results show that the bath can provide beneficial effects to the success of the protocol. We find evidence for a complex interplay between nonadiabaticity of the protocol due to the external drive and dissipation effects, which strongly depends on the detailed form of the qubit-bath interaction.

arXiv:1905.01762•Physical review. B./Physical review. B•Jul 16, 2019•8 citations

Reduced dynamics in spin-boson models: A method for both slow and fast bath

Andrei A. Golosov, Richard A. Friesner, Philip Pechukas

We study a model for treating dissipative systems, a one dimensional quantum system coupled to a harmonic bath. The dynamics of such a system can be described by Feynman’s path integral expression for the reduced density matrix. In this formulation the interaction of the system with the environment is stored in the influence functional. Recently we showed that fast environmental modes that give rise to correlations in the influence functional which are short range in time can be treated efficiently by a memory equation algorithm, which is a discretized version of a master equation. In this work we extend this approach to treat slow environmental modes as well, thereby efficiently linking adiabatic and nonadiabatic regimes. In this extended method the long range correlations in the influence functional arising from slow bath modes are taken into account through Stock’s semiclassical self-consistent-field approach.

The Journal of Chemical Physics•Feb 1, 2000•46 citations

Nonperturbative stochastic method for driven spin-boson model

Peter P. Orth, Adilet Imambekov, Karyn Le Hur

We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems, and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\"odinger equation that allows us to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work [P. P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. A 82, 032118 (2010)] by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known noninteracting-blip approximation in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence ${\ensuremath{\sigma}}^{x}(t)$ and of ${\ensuremath{\sigma}}^{z}(t)$ under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of ${\ensuremath{\sigma}}^{z}(t)$ without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at nonzero bias and different temperatures.

arXiv:1211.1201•Physical Review B•Jan 25, 2013•83 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
Quantum information processing has greatly increased interest in the phenomenon of environmentally-induced decoherence. The spin boson model is widely used to study the interaction between a spin-modelling a quantum particle moving in a double well potential-and its environment-modelled by a heat bath of harmonic oscillators. This paper extends a previous analysis of the static spin boson study to the driven spin boson case, with the derivation of an exact integro-differential equation for the time evolution of the propagator of the reduced spin density matrix. This is the first main result. By specializing to weak damping we then obtain the next result, a set of Bloch-Redfield equations for the equilibrium fixed spin initial condition. Finally we show that these equations can be used to solve the classic dissipative Landau-Zener problem and illustrate these solutions for the weak damping case. The effect of dissipation is seen to be minimised as the speed of passage is increased, implying that qubits need to be switched as fast as possible.
We study the dynamics of a two-level system described by a slowly varying Hamiltonian and weakly coupled to the Ohmic environment. We follow the Bloch-Redfield perturbative approach to include the effect of the environment on qubit evolution and take into account modification of the spectrum and matrix elements of qubit transitions due to the time dependence of the Hamiltonian. We apply this formalism to two problems. (1) We consider a qubit, or a spin-12\frac{1}{2}21​, in a rotating magnetic field. We show that once the rotation starts, the spin has a component perpendicular to the rotation plane of the field that initially wiggles and eventually settles to the value proportional to the product of angular rotation velocity of the field and the Berry curvature. (2) We reexamine the Landau-Zener transition for a system coupled to the environment at arbitrary temperature. We show that as temperature increases, the thermal excitation and relaxation become leading processes responsible for transition between states of the system. We also apply the Lindblad master equations to these two problems and compare results with those obtained from the Bloch-Redfield equations.
Abstract Using a dissipative quantum Rabi model, we study the dynamics of a slow qubit coupled to a fast quantum harmonic oscillator interacting with a bosonic bath from weak to strong and ultra-strong coupling regimes. Solving the quantum Heisenberg equations of motion, perturbative in the internal coupling between qubit and oscillator, we derive functional relationships directly linking the qubit coordinates in the Bloch sphere to oscillator observables. We then perform accurate time-dependent Matrix Product State simulations and compare our results both with the analytical solutions of the Heisenberg equations of motion, and with numerical solutions of a Lindblad master equation, perturbative in the external coupling between oscillator and environment. Indeed, we show that, up to the strong coupling regime, the qubit state accurately fulfils the derived functional relationships. We analyse in detail the case of a qubit starting with generic coordinates on the Bloch sphere of which we evaluate the three components of the Bloch vector through the averages of oscillator observables. Interestingly, a weak to intermediate oscillator coupling to the bath is able to simplify the Bloch vector evaluation since qubit-oscillator relationships are more immediate. Moreover, by monitoring the qubit fidelity with respect to free limit, we find the parameter regime where the combined effect of internal and external couplings is able to hinder the reliable evaluation of the qubit Bloch vector. Finally, in the ultra-strong coupling regime, non-Markovian effects become robust and the dynamics of qubit and oscillator are inextricably entangled making the qubit Bloch vector evaluation difficult. Graphical Abstract
We study the exact solution of the Schr\"odinger equation for the dissipative dynamics of a qubit, achieved by means of short iterative Lanczos method (SIL), which allows us to describe the qubit and the bath dynamics from weak to strong coupling regimes. We focus on two different models of a qubit in contact with the external environment: the first is the spin boson model (SBM), which gives a description of the qubit in terms of static tunneling energy and a bias field. The second model describes an externally driven qubit, where both the bias field and the tunneling rate are controlled by a time-dependent magnetic field, following a finite time protocol. We show that in the SBM case, our solution correctly describes the crossover from coherent to incoherent behavior of the magnetization, occurring at the Toulouse point. Furthermore, we show that the bath response dramatically changes during the system dynamics, going from nonresonant at small times to resonant behavior at long times. When the external driving field is present, for fixed values of the drive duration our results show that the bath can provide beneficial effects to the success of the protocol. We find evidence for a complex interplay between nonadiabaticity of the protocol due to the external drive and dissipation effects, which strongly depends on the detailed form of the qubit-bath interaction.
We study a model for treating dissipative systems, a one dimensional quantum system coupled to a harmonic bath. The dynamics of such a system can be described by Feynman’s path integral expression for the reduced density matrix. In this formulation the interaction of the system with the environment is stored in the influence functional. Recently we showed that fast environmental modes that give rise to correlations in the influence functional which are short range in time can be treated efficiently by a memory equation algorithm, which is a discretized version of a master equation. In this work we extend this approach to treat slow environmental modes as well, thereby efficiently linking adiabatic and nonadiabatic regimes. In this extended method the long range correlations in the influence functional arising from slow bath modes are taken into account through Stock’s semiclassical self-consistent-field approach.
We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems, and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\"odinger equation that allows us to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work [P. P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. A 82, 032118 (2010)] by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known noninteracting-blip approximation in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence ${\ensuremath{\sigma}}^{x}(t)$ and of ${\ensuremath{\sigma}}^{z}(t)$ under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of ${\ensuremath{\sigma}}^{z}(t)$ without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at nonzero bias and different temperatures.
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.