Chain-mapping techniques in combination with the time-dependent density matrix renormalization group are a powerful tool for the simulation of open-system quantum dynamics. For finite-temperature environments, however, this approach suffers from an unfavorable algorithmic scaling with increasing temperature. We prove that the system dynamics under thermal environments can be nonperturbatively described by temperature-dependent system-environmental couplings with the initial environment state being in its pure vacuum state, instead of a mixed thermal state. As a consequence, as long as the initial system state is pure, the global system-environment state remains pure at all times. The resulting speed-up and relaxed memory requirements of this approach enable the efficient simulation of open quantum systems interacting with highly structured environments in any temperature range, with applications extending from quantum thermodynamics to quantum effects in mesoscopic systems.
Chain-mapping techniques in combination with the time-dependent density matrix renormalization group are a powerful tool for the simulation of open-system quantum dynamics. For finite-temperature environments, however, this approach suffers from an unfavorable algorithmic scaling with increasing temperature. We prove that the system dynamics under thermal environments can be nonperturbatively described by temperature-dependent system-environmental couplings with the initial environment state being in its pure vacuum state, instead of a mixed thermal state. As a consequence, as long as the initial system state is pure, the global system-environment state remains pure at all times. The resulting speed-up and relaxed memory requirements of this approach enable the efficient simulation of open quantum systems interacting with highly structured environments in any temperature range, with applications extending from quantum thermodynamics to quantum effects in mesoscopic systems.
Chain-mapping techniques combined with the time-dependent density matrix renormalization group are powerful tools for simulating the dynamics of open quantum systems interacting with structured bosonic environments. Most interestingly, they leave the degrees of freedom of the environment open to inspection. In this work, we fully exploit the access to environmental observables to illustrate how the evolution of the open quantum system can be related to the detailed evolution of the environment it interacts with. In particular, we give a precise description of the fundamental physics that enables the finite temperature chain-mapping formalism to express dynamical equilibrium states. Furthermore, we analyze a two-level system strongly interacting with a super-Ohmic environment, where we discover a change in the spin-boson ground state that can be traced to the formation of polaronic states.
The chain mapping of structured environments is a most powerful tool for the simulation of open quantum system dynamics. Once the environmental bosonic or fermionic degrees of freedom are unitarily rearranged into a one dimensional structure, the full power of Density Matrix Renormalization Group (DMRG) can be exploited. Beside resulting in efficient and numerically exact simulations of open quantum systems dynamics, chain mapping provides an unique perspective on the environment: the interaction between the system and the environment creates perturbations that travel along the one dimensional environment at a finite speed, thus providing a natural notion of light-, or causal-, cone. In this work we investigate the transport of excitations in a chain-mapped bosonic environment. In particular, we explore the relation between the environmental spectral density shape, parameters and temperature, and the dynamics of excitations along the corresponding linear chains of quantum harmonic oscillators. Our analysis unveils fundamental features of the environment evolution, such as localization, percolation and the onset of stationary currents
The chain mapping of structured environments is a most powerful tool for the simulation of open quantum system dynamics. Once the environmental bosonic or fermionic degrees of freedom are unitarily rearranged into a one dimensional structure, the full power of Density Matrix Renormalization Group (DMRG) can be exploited. Beside resulting in efficient and numerically exact simulations of open quantum systems dynamics, chain mapping provides an unique perspective on the environment: the interaction between the system and the environment creates perturbations that travel along the one dimensional environment at a finite speed, thus providing a natural notion of light-, or causal-, cone. In this work we investigate the transport of excitations in a chain-mapped bosonic environment. In particular, we explore the relation between the environmental spectral density shape, parameters and temperature, and the dynamics of excitations along the corresponding linear chains of quantum harmonic oscillators. Our analysis unveils fundamental features of the environment evolution, such as localization, percolation and the onset of stationary currents.
All systems are open to an essentially uncontrollable environment that acts as a source of decoherence and dissipation. In some cases the environment's only effect is to add a weak relaxation mechanism and thus can be ignored for short timescales. In others, however, the presence of the environment can fundamentally alter the behavior of the system. Such is the case in mesoscopic superconductors where the environment can stabilize superconductivity and in spin-boson systems where the environment induces a localization transition. Likewise, in technological applications we are often interested in systems operating far from equilibrium. Here the environment might act as a particle reservoir or strong driving force. In all these examples, we need accurate methods to describe the influence of the environment on the system and to solve for the resulting dynamics or equilibrium states. In this thesis, we develop computational and conceptual approaches to efficiently simulate quantum systems in contact with an environment. Our starting point is the use of numerical renormalization techniques. Thus, we restrict our attention to one-dimensional lattices or small quantum systems coupled to an environment. We have developed several complementary algorithms: a superoperator renormalization algorithm for simulating real-time Markovian dynamics and for calculating states in thermal equilibrium; a blocking algorithm for simulating integro-differential equations with long-time memory; and a tensor network algorithm for branched lattices, which can be used to simulate strongly dissipative systems. Further, we provide support for an idea that to generically and accurately simulate the real-time dynamics of strongly dissipative systems, one has to include all or part of the environment within the simulation. In addition, we discuss applications and open questions.
The combination of chain-mapping and tensor-network techniques provides a powerful tool for the numerically exact simulation of open quantum systems interacting with structured environments. However, these methods suffer from a quadratic scaling with the physical simulation time, and therefore, they become challenging in the presence of multiple environments. This is particularly true when fermionic environments, well-known to be highly correlated, are considered. In this work, we first illustrate how a thermo-chemical modulation of the spectral density allows replacing the original fermionic environments with equivalent, but simpler, ones. Moreover, we show how this procedure reduces the number of chains needed to model multiple environments. We then provide a derivation of the fermionic Markovian closure construction, consisting of a small collection of damped fermionic modes undergoing a Lindblad-type dynamics and mimicking a continuum of bath modes. We describe, in particular, how the use of the Markovian closure allows for a polynomial reduction of the time complexity of chain-mapping based algorithms when long-time dynamics are needed.
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.