We analyze the static and dynamic properties of two Ising-coupled quantum spins embedded in a common bosonic bath as an archetype of dissipative quantum mechanics. First, we elucidate the ground-state phase diagram for an Ohmic and a sub-Ohmic bath using a combination of bosonic numerical renormalization group (NRG), analytical techniques, and intuitive arguments. Second, by employing the time-dependent NRG we investigate the system's rich dynamical behavior arising from the complex interplay between spin-spin and spin-bath interactions. Interestingly, spin oscillations can synchronize due to the proximity of the common non-Markovian bath and the system displays highly entangled steady states for certain nonequilibrium initial preparations. We complement our nonperturbative numerical results by exact analytical solutions when available and provide quantitative limits on the applicability of the perturbative Bloch-Redfield approach at weak coupling.
We analyze the static and dynamic properties of two Ising-coupled quantum spins embedded in a common bosonic bath as an archetype of dissipative quantum mechanics. First, we elucidate the ground-state phase diagram for an Ohmic and a sub-Ohmic bath using a combination of bosonic numerical renormalization group (NRG), analytical techniques, and intuitive arguments. Second, by employing the time-dependent NRG we investigate the system's rich dynamical behavior arising from the complex interplay between spin-spin and spin-bath interactions. Interestingly, spin oscillations can synchronize due to the proximity of the common non-Markovian bath and the system displays highly entangled steady states for certain nonequilibrium initial preparations. We complement our nonperturbative numerical results by exact analytical solutions when available and provide quantitative limits on the applicability of the perturbative Bloch-Redfield approach at weak coupling.
Employing the nonperturbative numerical renormalization group method, we study the dynamics of the spin-boson model, which describes a two-level system coupled to a bosonic bath with a spectral density J(omega) proportional to omega(s). We show that, in contrast with the case of Ohmic damping, the delocalized phase of the sub-Ohmic model cannot be characterized by a single energy scale only, due to the presence of a nontrivial quantum phase transition. In the strongly sub-Ohmic regime, s<<1, weakly damped coherent oscillations on short time scales are possible even in the localized phase--this is of crucial relevance, e.g., for qubits subject to electromagnetic noise.
The decoherence of a two-state system coupled with a sub-Ohmic bath is investigated theoretically by means of the perturbation approach based on a unitary transformation. It is shown that the decoherence depends strongly and sensitively on the structure of environment. Nonadiabatic effect is treated through the introduction of a function ${\ensuremath{\xi}}_{k}$ which depends on the boson frequency and renormalized tunneling. The results are as follows: (i) the nonequilibrium correlation function P(t), the dynamical susceptibility ${\ensuremath{\chi}}^{\ensuremath{''}}(\ensuremath{\omega})$, and the equilibrium correlation function C(t) are analytically obtained for $s\ensuremath{\leqslant}1$; (ii) the phase diagram of thermodynamic transition shows the delocalized-localized transition point ${\ensuremath{\alpha}}_{l}$ which agrees with exact results and numerical data from the numerical renormalization group; (iii) the dynamical transition point ${\ensuremath{\alpha}}_{c}$ between coherent and incoherent phase is explicitly given. A crossover from the coherent oscillation to incoherent relaxation appears with increasing coupling (for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{c}$, the coherent dynamics disappear); (iv) the Shiba's relation and sum rule are exactly satisfied when $\ensuremath{\alpha}\ensuremath{\leqslant}{\ensuremath{\alpha}}_{c}$; (v) an underdamping-overdamping transition point ${\ensuremath{\alpha}}_{c}^{*}$ exists in the function $S(\ensuremath{\omega})$. Consequently, the dynamical phase diagrams in both ohmic and sub-Ohmic case are mapped out. For $\ensuremath{\Delta}⪡{\ensuremath{\omega}}_{c}$, the critical couplings (${\ensuremath{\alpha}}_{l},{\ensuremath{\alpha}}_{c}$, and ${\ensuremath{\alpha}}_{c}^{*}$) are proportional to ${\ensuremath{\Delta}}^{1\ensuremath{-}s}$.
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.
Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.
By means of a unitary transformation, we study quantum phase transitions in the ground state of a two-qubit system interacting with a dissipative reservoir. First, the ground state phase diagram is analyzed in the presence of the Ohmic and sub-Ohmic bath using an analytic ground state wave function which takes into account the competition between intrasite tunneling and intersite correlation. The quantum critical point is determined as the transition point from non-degenerate to degenerate ground state and our calculated critical coupling strength αc agrees with that from the numerical renormalization group method. Moreover, by computing the entanglement entropy between the qubits and the bath in the ground state, we explore the nature of the quantum phase transition between the delocalized and localized states.
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.