In practical settings, the many-body quantum systems can not be truly isolated, underscoring the importance of open many-body quantum systems. A prominent example of an open quantum system is a quantum computer where qubits interact with each other and their surrounding environment. Studying open quantum systems brings its own set of challenges, especially in modeling the interaction between the system and environment. The study of open quantum systems is somewhat limited to the weak coupling between the system and the environment (Markov approximation). Recent developments utilizing process tensors have made handling the non-Markovian environment for studying open quantum systems accessible.I am currently working on the non-markovian dynamics of open quantum systems for my Master's thesis.
I started with the unbiased spin-boson model to investigate the effects of dissipation on a quantum system. This model describes a quantum particle in a dissipative bath of the harmonic oscillator. The Hamiltonian of this model is described below H=ΩSx+∑iSz(giai+g∗ia†i)+ωia†iai where Sx and Sz are spin operators, ai (a†i) are annihilation (creation) operators, wi is the ith frequency mode of the bath, and gi, the coupling strength between the system and the bath. The behavior of the bath is characterized by the spectral density function
For ohmic spectral density J(ω)=2αωexp(−ω/ωc), this model shows a quantum phase transition at critical value of the system-environment coupling α=αc. The time-evolving matrix product operator (TEMPO) algorithm has been employed for the calculation, which exploits the augmented density tensors (ADT) to capture the system's history over a finite bath memory time τc. Figure 1 shows the evolution of ⟨Sz(t)⟩, with an initial condition ⟨Sz(t=0)⟩=+12 and with no excitations in the environment. Before reaching the localized phase, there is a crossover at α≅0.5 from coherent decaying oscillation to incoherent decay. For α>0.5, ⟨Sz⟩ decays to zero asymptotically as ⟨Sz(t)⟩∝exp(−γt) as shown in figure 1.
However, at large α, ⟨Sz⟩ approaches a non-zero value asymptotically when the system localizes. The decay rate γ crosses zero at around αc≅1.25 with a 90% confidence interval as depicted in figure 2 where the inset shows the dependence of the decay rate γ with memory cutoff τc→KΔ. The system transitions from the delocalized phase to a localized phase at around αc≅1.25, consistent with the known result.