Dynamics of Open Quantum Systems

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+giai)+ωiaiai where Sx and Sz are spin operators, ai (ai) 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

Fig. 1 : The dynamics of the spin-1/2 system captured at K = 200 for the different values of α indicated. The dotted lines show the exponential fits to the data above α=0.5 .
J(ω)=i|gi|2δ(ωωi)

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.

Fig. 2 : Shows the extrapolated decay rate γ which crosses zero at around α1.25.The inset shows the dependence of the decay rate γ of the exponential fit on 1/K \& is utilized to extract the value of γ when K.

However, at large α, Sz approaches a non-zero value asymptotically when the system localizes. The decay rate γ crosses zero at around αc1.25 with a 90% confidence interval as depicted in figure 2 where the inset shows the dependence of the decay rate γ with memory cutoff τcKΔ. The system transitions from the delocalized phase to a localized phase at around αc1.25, consistent with the known result.