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:1607.03920 · Physical Review A

Dynamical renormalization-group approach to the spin-boson model

Hassan Shapourian

Published Mar 14, 2016• 11 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.93.032119

Related papers

Spin-boson model under dephasing: Markovian versus non-Markovian dynamics

Naushad Ahmad Kamar, Daniel A. Paz, Mohammad F. Maghrebi

The spin-boson model, describing a two-level system coupled to a bosonic bath, is extensively studied as a paradigmatic dissipative quantum system, exhibiting rich dynamical behavior and even a localization transition in the strong coupling regime. Here, we additionally consider dephasing as a source of Markovian dissipation on top of the non-Markovian dynamics due to an Ohmic bath, and investigate the dynamics of the spin. We show that the characteristic frequency of the spin dynamics, while strongly renormalized by the bosonic bath, changes in a simple fashion (or does not change at all) with dephasing. To obtain these results, we develop an exact nonperturbative method known as the stochastic Schr\"odinger equation, mimicking the Ohmic bath via a stochastic magnetic field combined with the Lindblad quantum master equation due to dephasing, which allows us to numerically compute the dynamics. Furthermore, we derive weak-coupling analytic results utilizing the well-known noninteracting blip approximation. Our findings are relevant to quantum simulation of the spin-boson model in the regime of strong coupling in trapped ions and circuit QED architectures among others.

arXiv:2305.00110•Physical review. B./Physical review. B•Aug 12, 2024•6 citations

Dynamics of the two-spin spin-boson model with a common bath

Tianrui Deng, Yiying Yan, Lipeng Chen +1 more

Dynamics of the two-spin spin-boson model in the presence of Ohmic and sub-Ohmic baths is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Its accuracy in dynamics simulations of the two-spin SBM is improved significantly over the single D1 Ansatz, especially in the weak to moderately strong coupling regime. Validity of the multi-D1 Ansatz for various coupling strengths is also systematically examined by making use of the deviation vector which quantifies how faithfully the trial state obeys the Schrödinger equation. The time evolution of population difference and entanglement has been studied for various initial conditions and coupling strengths. Careful comparisons are carried out between our approach and three other methods, i.e., the time-dependent numerical renormalization group (TD-NRG) approach, the Bloch-Redfield theory, and a method based on a variational master equation. For strong coupling, the multi-D1 trial state yields consistent results as the TD-NRG approach in the Ohmic regime while the two disagree in the sub-Ohmic regime, where the multi-D1 trial state is shown to be more accurate. For weak coupling, the multi-D1 trial state agrees with the two master-equation methods in the presence of both Ohmic and sub-Ohmic baths, but shows considerable differences with the TD-NRG approach in the presence of a sub-Ohmic bath, calling into question the validity of the TD-NRG results at long times in the literature.

The Journal of Chemical Physics•Apr 8, 2016•42 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

Spin precession and real-time dynamics in the Kondo model: Time-dependent numerical renormalization-group study

Frithjof B. Anders, Avraham Schiller

A detailed derivation of the recently proposed time-dependent numerical renormalization-group (TD-NRG) approach to nonequilibrium dynamics in quantum-impurity systems is presented. We demonstrate that the method is suitable for fermionic as well as bosonic baths. Comparisons with exact analytical results for the charge relaxation in the resonant-level model and for dephasing in the spin-boson model establish the accuracy of the method. The real-time dynamics of a single spin coupled to each type of bath is investigated. We use the TD-NRG to calculate the spin relaxation and spin precession of a single Kondo impurity. The short- and long-time dynamics are studied as a function of temperature in the ferromagnetic and antiferromagnetic regimes. The short-time dynamics agrees very well with analytical results obtained at second order in the exchange coupling $J$. In the ferromagnetic regime, the transient spin decay is described by the scaling variable $x=2{\ensuremath{\rho}}_{F}\ensuremath{\mid}J(T)\ensuremath{\mid}Tt$. In the antiferromagnetic regime, the long-time decay is governed for $T<{T}_{K}$ by the Kondo time scale $1∕{T}_{K}$. Here ${\ensuremath{\rho}}_{F}$ is the conduction-electron density of states, ${T}_{K}$ is the Kondo temperature, and $J(T)$ is the effective exchange coupling at temperature $T$. Results for spin precession are obtained by rotating the external magnetic field from the $x$ axis to the $z$ axis.

arXiv:cond-mat/0604517•Physical Review B•Dec 18, 2006•250 citations

Dynamical scaling in the Ohmic spin-boson model studied by extended hierarchical equations of motion

Qianlong Wang, Zhihao Gong, Chenru Duan +2 more

Through a decomposition of the bath correlation function, the hierarchical equations of motion are extended to the Ohmic spin-boson model at zero temperature. For two typical cutoff functions of the bath spectral density, the rate kernel of spin dynamics is numerically extracted by a time-convolution equation of the average magnetic moment. A characteristic time is defined accordingly as the inverse of the zeroth-order moment of the rate kernel. For a given Kondo parameter in the incoherent regime, the time evolution of average magnetic moments gradually collapses onto a master curve after rescaling the time variable with the characteristic time. The rescaled spin dynamics is nearly independent of the cutoff frequency and the form of cutoff functions. For a given cutoff frequency, the characteristic time with the change of the Kondo parameter is fitted excellently as a function of the renormalized tunneling amplitude. Despite a significant difference in definition, our result is in good agreement with the characteristic time of the noninteracting blip approximation.

The Journal of Chemical Physics•Feb 28, 2019•18 citations

Dynamics of open quantum systems

Herbert Schoeller

A tutorial introduction is presented for the calculation of the time dynamics for models of dissipative quantum mechanics where a small quantum system is coupled to noninteracting bosonic or fermionic reservoirs. We discuss the basics and applications of the real-time renormalization group method (RTRG), where the Fourier variable $E$ conjugate to the time variable $t$ is used as flow parameter. It is shown how a well-controlled perturbation theory can be set up in the coupling to the bath where secular and logarithmic terms are resummed to all order of perturbation theory. This is achieved by a set of perturbative RG equations for the effective Liouvillian and the effective vertices. Truncating the RG equations we show how the time dynamics can be obtained analytically for exponentially small, intermediate and exponentially large times. Basic applications are discussed for particle and potential energy exchange within the interacting resonant level model, spin exchange for the nonequilibrium Kondo model, and heat exchange within the ohmic spin boson model at zero bias.

arXiv:1802.10014•Feb 27, 2018•10 citations
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α\ensuremath≤1/2\ensuremath{\alpha}\ensuremath{\le}1/2\ensuremathα\ensuremath≤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 model, describing a two-level system coupled to a bosonic bath, is extensively studied as a paradigmatic dissipative quantum system, exhibiting rich dynamical behavior and even a localization transition in the strong coupling regime. Here, we additionally consider dephasing as a source of Markovian dissipation on top of the non-Markovian dynamics due to an Ohmic bath, and investigate the dynamics of the spin. We show that the characteristic frequency of the spin dynamics, while strongly renormalized by the bosonic bath, changes in a simple fashion (or does not change at all) with dephasing. To obtain these results, we develop an exact nonperturbative method known as the stochastic Schr\"odinger equation, mimicking the Ohmic bath via a stochastic magnetic field combined with the Lindblad quantum master equation due to dephasing, which allows us to numerically compute the dynamics. Furthermore, we derive weak-coupling analytic results utilizing the well-known noninteracting blip approximation. Our findings are relevant to quantum simulation of the spin-boson model in the regime of strong coupling in trapped ions and circuit QED architectures among others.
Dynamics of the two-spin spin-boson model in the presence of Ohmic and sub-Ohmic baths is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Its accuracy in dynamics simulations of the two-spin SBM is improved significantly over the single D1 Ansatz, especially in the weak to moderately strong coupling regime. Validity of the multi-D1 Ansatz for various coupling strengths is also systematically examined by making use of the deviation vector which quantifies how faithfully the trial state obeys the Schrödinger equation. The time evolution of population difference and entanglement has been studied for various initial conditions and coupling strengths. Careful comparisons are carried out between our approach and three other methods, i.e., the time-dependent numerical renormalization group (TD-NRG) approach, the Bloch-Redfield theory, and a method based on a variational master equation. For strong coupling, the multi-D1 trial state yields consistent results as the TD-NRG approach in the Ohmic regime while the two disagree in the sub-Ohmic regime, where the multi-D1 trial state is shown to be more accurate. For weak coupling, the multi-D1 trial state agrees with the two master-equation methods in the presence of both Ohmic and sub-Ohmic baths, but shows considerable differences with the TD-NRG approach in the presence of a sub-Ohmic bath, calling into question the validity of the TD-NRG results at long times in the literature.
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.
A detailed derivation of the recently proposed time-dependent numerical renormalization-group (TD-NRG) approach to nonequilibrium dynamics in quantum-impurity systems is presented. We demonstrate that the method is suitable for fermionic as well as bosonic baths. Comparisons with exact analytical results for the charge relaxation in the resonant-level model and for dephasing in the spin-boson model establish the accuracy of the method. The real-time dynamics of a single spin coupled to each type of bath is investigated. We use the TD-NRG to calculate the spin relaxation and spin precession of a single Kondo impurity. The short- and long-time dynamics are studied as a function of temperature in the ferromagnetic and antiferromagnetic regimes. The short-time dynamics agrees very well with analytical results obtained at second order in the exchange coupling JJJ. In the ferromagnetic regime, the transient spin decay is described by the scaling variable $x=2{\ensuremath{\rho}}_{F}\ensuremath{\mid}J(T)\ensuremath{\mid}Tt$. In the antiferromagnetic regime, the long-time decay is governed for T<TKT<{T}_{K}T<TK​ by the Kondo time scale 1∕TK1∕{T}_{K}1∕TK​. Here ${\ensuremath{\rho}}_{F}$ is the conduction-electron density of states, TK{T}_{K}TK​ is the Kondo temperature, and J(T)J(T)J(T) is the effective exchange coupling at temperature TTT. Results for spin precession are obtained by rotating the external magnetic field from the xxx axis to the zzz axis.
Through a decomposition of the bath correlation function, the hierarchical equations of motion are extended to the Ohmic spin-boson model at zero temperature. For two typical cutoff functions of the bath spectral density, the rate kernel of spin dynamics is numerically extracted by a time-convolution equation of the average magnetic moment. A characteristic time is defined accordingly as the inverse of the zeroth-order moment of the rate kernel. For a given Kondo parameter in the incoherent regime, the time evolution of average magnetic moments gradually collapses onto a master curve after rescaling the time variable with the characteristic time. The rescaled spin dynamics is nearly independent of the cutoff frequency and the form of cutoff functions. For a given cutoff frequency, the characteristic time with the change of the Kondo parameter is fitted excellently as a function of the renormalized tunneling amplitude. Despite a significant difference in definition, our result is in good agreement with the characteristic time of the noninteracting blip approximation.
A tutorial introduction is presented for the calculation of the time dynamics for models of dissipative quantum mechanics where a small quantum system is coupled to noninteracting bosonic or fermionic reservoirs. We discuss the basics and applications of the real-time renormalization group method (RTRG), where the Fourier variable EEE conjugate to the time variable ttt is used as flow parameter. It is shown how a well-controlled perturbation theory can be set up in the coupling to the bath where secular and logarithmic terms are resummed to all order of perturbation theory. This is achieved by a set of perturbative RG equations for the effective Liouvillian and the effective vertices. Truncating the RG equations we show how the time dynamics can be obtained analytically for exponentially small, intermediate and exponentially large times. Basic applications are discussed for particle and potential energy exchange within the interacting resonant level model, spin exchange for the nonequilibrium Kondo model, and heat exchange within the ohmic spin boson model at zero bias.