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 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 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.
The Lindblad theory of open quantum systems has been successfully applied in various fields ranging from quantum optics, condensed matter physics, and quantum thermodynamics to quantum chemistry. Nevertheless, there are situations where we can compare its predictions to those of exact methods or numerically exact results and find discrepancies between them. One example is the case where a single harmonic oscillator is coupled to a thermal bath of independent harmonic oscillators. In this case, for certain parameters, the exact master equation describes equilibration of the system to a thermal state with a temperature different from the one predicted by the Lindblad master equation. The Lindblad theory, based on a weak coupling between the system and environment as well as the Born-Markov approximation, predicts that the system will evolve to a thermal state whose temperature is equal to the initial temperature of the bath. This result is a consequence of the Born approximation, which supposes that throughout the evolution, the bath only fluctuates around its initial equilibrium state. In this master thesis we study the validity of the Born approximation as well as its consistency with other approximations. By starting from a Hamiltonian for the total system we show that, if the the Born approximation is not invoked, the dynamics of the system is determined by a hierarchy of equations for matrices which can be used to reconstruct the reduced density matrix of the system. This hierarchy approach is fundamentally different from the standard approach, which describes the system dynamics with a weak-coupling master equation, because the hierarchy does not restrict the amount of entanglement possible between the system and the environment. Furthermore, we show how invoking the Born approximation reduces to the usual result and how a generalised ansatz, which we refer to as a generalised Born approximation, reduces the hierarchy to a master equation where the bath temperature is an additional degree of freedom. Finally, we analyse the conservation of the energy of the total system and its relationship to the Markov and secular approximations, and discuss whether imposing it as a condition can fix the time-dependent temperature.
In various fields from quantum physics to biology, the open quantum dynamics of a system consisting of interacting subsystems emphasizes its fundamental functionality. The local approach, deriving a dissipator in a master equation by ignoring the inter-subsystem interaction, has been widely used to describe the reduced dynamics due to its robustness to keep the positivity of a density operator. However, one critique is that a stationary state obtained by the approach in the limit of weak system-environment coupling is written in the form of the Gibbs state for the partial Hamiltonian by excluding the inter-subsystem interaction from the total one of the relevant system. As an alternative, the global approach, deriving a dissipator with including the inter-subsystem interaction, under the Born--Markov and secular approximations has attracted much attention, and there is debate concerning its violation of positivity in the short-term region and/or limited parameter region for the Bohr frequencies of the subsystems. In this study, we present a formalism that leads to the time-convolutionless (time-local) master equation obtained by extending the global approach beyond the Born-Markov and secular approximations. We apply it to the excitation energy transfer between interacting sites in which only the terminal site weakly interacts with a bosonic environment of finite temperature in a manner beyond the rotating-wave approximation. We find that the formulation (1) gives the short-time behavior while preserving positivity, (2) shows the oscillatory features that the secular approximation would obscure, and (3) leads to a stationary state very near to the Gibbs state for the total Hamiltonian of the relevant system.
In various fields from quantum physics to biology, the open quantum dynamics of a system consisting of interacting subsystems emphasizes its fundamental functionality. The local approach, deriving a dissipator in a master equation by ignoring the intersubsystem interaction, has been widely used to describe the reduced dynamics due to its robustness to keep the positivity of a density operator. However, one critique is that a stationary state obtained by the approach in the limit of weak system-environment coupling is written in the form of the Gibbs state for the partial Hamiltonian by excluding the intersubsystem interaction from the total one of the relevant system. As an alternative, the global approach, deriving a dissipator with including the intersubsystem interaction, under the Born-Markov and secular approximations has attracted much attention, and there is debate concerning its violation of positivity in the short-time region and/or limited parameter region for the Bohr frequencies of the subsystems. In this paper, we present a formalism that leads to the time-convolutionless (time-local) master equation obtained by extending the global approach beyond the Born-Markov and secular approximations. We apply it to the excitation energy transfer between interacting sites in which only the terminal site weakly interacts with a bosonic environment of finite temperature in a manner beyond the rotating-wave approximation. We find that the formulation (1) gives the short-time behavior while preserving positivity, (2) shows the oscillatory features that the secular approximation would obscure, and (3) leads to a stationary state very near to the Gibbs state for the total Hamiltonian of the relevant system.
We present a numerical path-integral iteration scheme for the low-dimensional reduced density matrix of a time-dependent quantum dissipative system. Our approach simultaneously accounts for the combined action of a microscopically modeled pure-dephasing-type coupling to a continuum of harmonic oscillators representing, e.g., phonons, and further environmental interactions inducing non-Hamiltonian dynamics in the inner system represented, e.g., by Lindblad-type dissipation or relaxation. Our formulation of the path-integral method allows for a numerically exact treatment of the coupling to the oscillator modes and moreover is general enough to provide a natural way to include Markovian processes that are sufficiently described by rate equations. We apply this new formalism to a model of a single semiconductor quantum dot which includes the coupling to longitudinal acoustic phonons for two cases: (a) external laser excitation taking into account a phenomenological radiative decay of the excited dot state and (b) a coupling of the quantum dot to a single mode of an optical cavity taking into account cavity photon losses.
We introduce and study the adiabatic dynamics of free-fermion models subject to a local Lindblad bath and in the presence of a time-dependent Hamiltonian. The merit of these models is that they can be solved exactly, and will help us to study the interplay between nonadiabatic transitions and dissipation in many-body quantum systems. After the adiabatic evolution, we evaluate the excess energy (the average value of the Hamiltonian) as a measure of the deviation from reaching the final target ground state. We compute the excess energy in a variety of different situations, where the nature of the bath and the Hamiltonian is modified. We find robust evidence of the fact that an optimal working time for the quantum annealing protocol emerges as a result of the competition between the nonadiabatic effects and the dissipative processes. We compare these results with the matrix-product-operator simulations of an Ising system and show that the phenomenology we found also applies for this more realistic case.