We propose a numerical technique suitable for simulating the dynamics of reduced density matrix of a qubit interacting with its environment. Our approach, based on a combination of short-iterative Lanczos method (SIL) and a flexible truncation scheme, allows to include in the physical description multiple-excitation processes, beyond weak coupling and Markov approximations. We perform numerical simulations of two different model Hamiltonians, that are relevant in the field of adiabatic quantum computation (AQC), and we show that our technique is able to recover the correct thermodynamic behavior of the qubit-bath system, from weak to intermediate coupling regime.
We propose a numerical technique suitable for simulating the dynamics of reduced density matrix of a qubit interacting with its environment. Our approach, based on a combination of short-iterative Lanczos method (SIL) and a flexible truncation scheme, allows to include in the physical description multiple-excitation processes, beyond weak coupling and Markov approximations. We perform numerical simulations of two different model Hamiltonians, that are relevant in the field of adiabatic quantum computation (AQC), and we show that our technique is able to recover the correct thermodynamic behavior of the qubit-bath system, from weak to intermediate coupling regime.
We propose a numerical technique based on a combination of short-iterative Lanczos and exact diagonalization methods, suitable for simulating the time evolution of the reduced density matrix of a single qubit interacting with an environment. By choosing a mode discretization method and a flexible bath states truncation scheme, we are able to include in the physical description multiple-excitation processes, beyond weak coupling and Markov approximations. We apply our technique to the simulation of three different model Hamiltonians, which are relevant in the field of adiabatic quantum computation. We compare our results with those obtained on the basis of the widely used Lindblad master equation, as well as with well-known exact and approximated approaches. We show that our method is able to recover the thermodynamic behavior of the qubit-bath system, beyond the Born-Markov approximation. Finally, we show that even in the case of the adiabatic quantum annealing of a single qubit the bath can be beneficial in reaching the reduced system ground state.
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.
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.
We demonstrate an exact mapping of a class of models of two interacting qubits in thermal reservoirs to two separate problems of spin–boson-type systems. Based on this mapping, exact numerical simulations of the qubits dynamics can be performed, beyond the weak system–bath coupling limit and the Markovian approximation. Given the time evolution of the system population and coherences, we study as an application the dynamics of entanglement between the pair of qubits immersed in boson thermal baths, showing a rich phenomenology, including an intermediate oscillatory behavior, the entanglement sudden birth, sudden death and revival. We find that the occurrence of entanglement sudden death in this model depends on the portion of the zero and double excitation states in the subsystem initial state. In the long-time limit, analytic expressions are presented at weak system–bath coupling, for a range of relevant qubit parameters.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
In this PhD thesis we have been concerned with dissipative and decoherence processes which take place specifically for quantum systems in a dissipative bistable potential. This potential describes a Hamiltonian system whose spectrum is composed by a ladder of levels, grouped as doublets under the potential barrier and being almost equally spaced above it. The system interacts with a thermal environment, fully characterized by the so-called bath spectral density. The dynamical quantity we looked at is the so-called reduced density matrix (RDM), obtained by tracing out the full density matrix the environmental degrees of freedom. The matrix elements for the RDM are exactly evaluated by means of the Feynman-Vernon real-time path-integral technique. Depending on which regime of parameters we are interested in, i.e. temperature and coupling strength with the reservoir, a few or several energy levels contribute to the dynamics. In particular, a two-level system (TLS) is first analyzed, corresponding to the regime where only the lowest doublet is energetically accessible. We discussed in the beginning the case of a TLS interacting with Ohmic and super-Ohmic baths and then the situation in which a symmetric TLS is assumed to be coupled with an external detector, like a dc-SQUID, the latter being coupled to an Ohmic environment. This model provides a better description of a real qubit, as frequently discussed in various experiments. In the former case, we developed a theory called Weakly-Interacting Blip Approximation (WIBA), which is nonperturbative in the coupling strength as well as in the tunneling frequency and gives reliable results over a wide range of temperatures and coupling strength to the thermal environment. In the latter case, the dynamics was evaluated upon starting from the Non-Interacting Blip Approximation (NIBA) and approximating the NIBA kernels up to first order in the detector-bath coupling strength. Analytical results, valid also at finite detuning, were obtained and we showed that our theory correctly reproduces the onset of multiple oscillation frequencies in the dynamics due to the entanglement between TLS and detector. Finally, we discussed the more general case of an M-level system coupled with an Ohmic environment. This situation describes the case of a bistable potential where the temperature is high enough that higher-lying energy levels are being populated and the system cannot be reduced to a simple TLS. Many are the experimental situations which can be described by our theory. We developed a novel approximation called Weakly-Interacting VR-Blip Approximation (VR-WIBA) for the dynamics in a driven dissipative bistable potential. In the regimes of low to moderate temperatures and coupling strengths, we exploited a separation of time scales between tunneling and vibrational relaxation events and we could derive analytical kernels for the intra- and inter-well dynamics.